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1. Introduction 

Our present understanding of the evolution of the universe is based on 
the Standard Hot Big Bang Model. The expansion of the universe, the 
synthesis of the light elements and the Cosmic Microwave Background 
(CMB) radiation arc the pillars of this Model. Its definite observational 
confirmation was the CMB accidental detection in 1965 by Penzias & 
Wilson, who received the Nobel Prize for this discovery. In the 1990s 
the COBE^ (COsmic Background Explorer) satellite improved on the 
previous observations and measured an almost isotropic blackbody ra- 
diation with temperature Tcmb = 2.726 ±0.010 K (Mather et al. 1994), 
together with temperature anisotropics on angular scales 90° (Smoot 
et al. 1992). These discoveries warranted yet another Nobel Prize (in 
2006) to the principal investigators of the mission, Smoot and Mather. 
The isotropy of the microwave background indicates that on large scales 
(> 200 Mpc) the universe is indeed very smooth, as postulated by 
Einstein in his "cosmological principle". On small scales, in contrast, 
it presents inhomogeneities, from planets and stars, to galaxies, clus- 
ters and super-clusters of galaxies. These structures are not uniformly 
distributed, but rather show some spatial correlation, and regions of 
space almost totally devoid of galaxies are alternated to high den- 
sity regions. The commonly adopted theory for the formation of these 
structures is the gravitational instability scenario, in which primordial 
density perturbations grow through gravitational Jeans instability to 
form all the structures we observe today. The evidence for the exis- 
tence of such perturbations is based on combined CMB observations by 
BOOMERanG^ (Balloon Observations Of Millimetric Extragalactic Radi- 
ation and Geomagnetics), MAXIMA^ (Millimeter Anisotropy eXperiment 
IMaging Array), DASI^ (Degree Angular Scale Interferometer) and, 
more recently, WMAP^ (Wilkinson Microwave Anisotropy Probe) which 
provided the power spectrum of temperature anisotropics on a wide 
range of angular scales (Hu et al. 2001; Sievers et al. 2003; Spergel 
et al. 2003; Spergel et al. 2007). The anisotropics at scales < 90° 
are interpreted as the result of perturbations in the energy density at 
the decoupling epoch and hence are directly related to the primordial 
density fluctuations from which all present structures originate. 

While the observation of the very high-redshift universe, z k, 1000, 
is possible through the CMB, there is a lack of observational data in the 

^ http:/ /lambda.gsfc. nasa.gov/product/cobe 
^ http:/ /cmb.phys.cwrii.cdu/boomcrang 
^ http:/ /cosmology.berkeley.edu/group/cmb 

http: / / astro.uchicago.edu/dasi 
^ http://niap.gsfc.nasa.gov 
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redshift interval 1 < z < 1000, with only a few objects with z 6.5 — 8 
(Hu et al. 2002; Kodaira et al. 2003; Kncib ct al. 2004; Bouwens et 
al. 2004; Malhotra et al. 2005; Mobashcr ct al. 2005; lye et al. 2006; 
Schaerer & Pello 2005; Ota et al. 2007). A deep, narrow J-band search 
for Lya emitting galaxies at z 9 by Willis k, Courbin (2005) has 
given no result, but the NICMOS field at z « 10 has revealed 3 possible 
galaxies (Bouwens et al. 2005). Another tentative detection of a galaxy 
at z 10 has been reported (Pello et al. 2004a,b), but this has not 
been confirmed by an analysis of the same data from a diflFerent group 
(Weatherley, Warren & Babbcdgc 2004). Other two lensed candidates 
at z « 10 have been reported by Stark et al. (2007b). 

The state of the art of the telescopes observing in the NIR/optical/UV 
range is represented by HST^ (Hubble Space Telescope), that has al- 
lowed to take high resolution images ranging from supernova (SN) 
explosions, to nebulae with ongoing star formation, to the Hubble Deep 
Field and Ultra Deep Field. Additional important contributions are 
brought to the field by the Spitzer^ satellite, a space-based infrared 
telescope launched in 2003 working in the 3-180 jiia range. Thanks to 
its exquisite sensitivity in this band, it has allowed to better understand 
the evolution of galaxies up to z ~ 5 through a variety of mid- and far- 
infrared imaging surveys and to constrain the nature and amplitude 
of the Cosmic Infrared Background (CIB). In order to observe radia- 
tion from objects at even higher redshift, telescopes with exceptional 
sensitivity in the IR and radio bands are needed. JWST^ (James Webb 
Space Telescope), for example, with its nJy sensitivity in the 1 — 10 /xm 
infrared regime, is ideally suited for probing optical-UV emission from 
sources sX z > 10. Similarly, the planned generation of radio telescopes 
as SKA^ (Square Kilometer Array), LOFAR^*^ (LOw Frequency ARray), 
21cniA^^ (21cm Array; Pen, Wu & Peterson 2004; Peterson, Pen & 
Wu 2005) and MWA^^ (Murchison Widefield Array) will open a new 
observational window on the high redshift universe. 

Thus, in the near future, we should be able to image the first sources 
of light that had formed in the universe, and unveil their nature. Nowa- 
days, observational evidence suggests that they are stellar type objects, 
although the existence of primordial mini-quasars is not completely 
ruled out. 



® http://www.stsci.edu/hst 

http://www.spitzer.caltech.edu/spitzer/index.shtml 

* http://ngst.gsfc.nasa.gov 

® http://www.skatclcscope.org 
^'^ http://www.lofar.org 

http://astrophysics.phys.cmu.edu/ jbp 
http://web.haystack.mit.edu/arrays/MWA/ 
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In the standard cosmological hierarchical scenario for structure for- 
mation, the objects which form first are predicted to have masses 
corresponding to virial temperatures Tyir < 10*^ K. Once the gas has 
virialized in the potential wells of pre-existing dark matter halos, ad- 
ditional cooling is required to further collapse the gas and form stars. 
For a gas of primordial composition at such low temperatures the main 
coolant is molecular hydrogen. The typical primordial H2 fraction is 
usually lower than the one required for the formation of such objects, 
but during the collapse phase the H2 content can reach suitable values, 
depending on the mass of the object itself. Once the gas has collapsed 
and cooled, star formation (SF) is ignited. The study of the star for- 
mation process in primordial objects is still in its early stages, and 
although large improvements have been done in both analytical and 
numerical approaches, a solid physical scenario requires substantial 
further investigation. 

Once the first sources have formed, their mass deposition, energy 
injection and emitted radiation can deeply affect the subsequent galaxy 
formation process and influence the evolution of the Intergalactic Medium 
(IGM) via a number of so-called "feedback" effects. Although a rigorous 
classification of the various effects is not feasible, they can be divided 
into three broad classes: radiative, mechanical and chemical feedback. 
Into the first class fall all those effects associated, in particular, with 
ionization/dissociation of hydrogen atoms/molecules; the second class 
is produced by the mechanical energy injection of massive stars in form 
of winds or supernova explosions; the chemical feedback is, in contrast, 
related to the postulated existence of a critical metallicity governing 
the cosmic transition from very massive stars to "normal" stars. 

In this review we will summarize the recent achievements in under- 
standing the first structure formation process. 



2. First Structure Formation: Overview 

In this Section we will describe the main processes leading from the 
evolution of primordial density perturbations to the formation of the 
first structures in the universe. 

2.1. The Cosmological Model 

Throughout this review we will usually refer to the popular ACDM 
model, in which dark matter is composed of cold, weakly interacting, 
massive particles. The cosmological model is completely defined once 
the value of the following parameters is specified: the adimensional 
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density of the universe, = + ^^A (where and are the con- 
tribution from matter and vacuum, respectively), the Hubble constant, 
Hq = lOO/i km s""*^ Mpc~^, the adimensional baryon density, ^^5, the 
rms mass fluctuations on 8h~^ Mpc scale, as, and the spectral index 
of the primordial density fluctuation, n. Recent observations have set 
stringent constraints on the values allowed for these parameters. Studies 
of the temperature anisotropics in the CMB strongly suggest that the 
universe is fiat, f^o ^ 1 (Balbi et al. 2000, 2001; O'Meara et al. 2001; 
Pryke et al. 2002), with a matter contribution of flm ~ 0.10 — 0.36. The 
latter values are derived from the local abundance of galaxy clusters 
(Pen 1998; Viana k Liddle 1999; Schuecker et al. 2003), their redshift 
evolution (Bahcall & Fan 1998; Eke et al. 1998; Borgani et al. 2001), 
mass-to-light ratio of groups of galaxies (Hoekstra et al. 2001), X-ray 
gas mass fraction of clusters of galaxies (Allen, Schmidt & Fabian 2002), 
clustering measurements from large galaxy surveys (e.g. Blake et al. 
2007) and high-2; supernovae (Schmidt et al. 1998; Perlmutter et al. 
1999; but see also Rowan- Robinson 2002). The baryon contribution 
to matter, Q^/i^ ^ 0.020 — 0.030, is derived through measurements of 
CMB anisotropics (Balbi et al. 2000, 2001; Pryke et al. 2002), clustering 
measurements from large galaxy surveys (e.g. Blake et al. 2007) and 
D/H abundance in quasars spectra (Tytler et al. 2000; Pettini & Bowen 
2001; Levshakov et al. 2002; Crighton et al. 2004). AH the estimates 
of the Hubble constant are based on secondary distance indicators, 
such as the TuUy-Fisher relation, the fundamental plane of elliptical 
galaxies. Type la supernovae and surface brightness fluctuations, and 
give h 0.5 - 0.8 (Theureau et al. 1997; Tonry et al. 1997; Jha et 
al. 1999; Mould et al. 2000; Freedman et al. 2001; Gibson & Stetson 
2001; Liu k Graham 2001; Tutui et al. 2001; Willick & Batra 2001). 
Local and high-z evolution of galaxy clusters is used to set a limit 
to (78 ^ 0.7 - 1.2 (Bahcah & Fan 1998; Eke et al. 1998; Pen 1998; 
Viana & Liddle 1999; Borgani et al. 2001; Hoekstra, Yee & Gladders 
2002). Finally, studies of the temperature anisotropics in the CMB 
give n 1.01 - 1.08 (Balbi et al. 2000, 2001; Pryke et al. 2002). 
These parameters have been confirmed by the CMB data from WMAP 
(e.g. Spergel et al. 2003, 2006; Dunklcy et al. 2008). But, although the 
excess power at I « 2000 — 6000 reported by CBI^"^ (Cosmic Background 
Imager) is roughly consistent with a secondary contribution resulting 
from the Sunyaev-Zeldovich effect, this requires a normalization for the 
matter power spectrum higher than measured by other means (Sievers 
et al. 2003). 



http://www.astro.caltech.edu/ tjp/CBI 
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The most popular, or concordance, model is a ACDM with = 
1.02 ± 0.02, Qa = 0.73 ± 0.04, = 0.27 ± 0.04, = 0.044 ± 0.004, 
h = 0.7ltg;|j^, 0-8 = 0.84 ± 0.04 and n = 0.93 ± 0.03 (Bennett et al. 
2003). The above parameter set, determined from CMB experiments, 
is in impressively good agreement with that obtained independently by 
experiments based on high-redshift Type la supernovae (Knop et al. 
2003) and galaxy cluster evolution data. Throughout the review we will 
adopt these values. Slight changes to the above values are constantly 
being made as WMAP is collecting additional data. The results after 5 
years of mission have been released very recently (see e.g. Hinshaw et 
al. 2008). The reader is then referred to the WMAP web page for further 
updates. In Fig. 1 expected values for the cosmological parameters 
as measured by WMAP are plotted, while Fig. 2 shows the most recent 
comparison (i.e. after the 5-th year of WMAP observations) between the 
CMB temperature fluctuations predicted by such ACDM model and 
those measured by CMB experiments. 
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Figure 2. CMB temperature fluctuations as measured by difi'erent missions (sym- 
bols). The curve is the best-fit ACDM model to the WMAP data (see Nolta et al. 2008 
for details). 

2.2. Linear Growth of Fluctuations 

The commonly adopted theory for structure formation is the gravita- 
tional instability scenario, in which primordial density perturbations 
grow through gravitational Jeans instability to form all the structures 
we observe today. The most favored model for the origin of these pertur- 
bations is inflation. In this scenario, the universe expands exponentially 
for a brief period of time at early epochs {t ~ 10~^^ — 10^'^^ s), and, at 
the end of this inflationary period, it is highly homogeneous on large 
scales, but is locally perturbed as a consequence of quantistic fluctua- 
tions. To derive the evolution of these primordial density fluctuations 
into bound objects, we can proceed as follows. 

Let us describe the universe in terms of a fluid made of collisionless 
dark matter and baryons, with a mean mass density p. At any time and 
location, the mass density can be written as p{x,t) = p{t)[l + 6{x,t)], 
where x indicates the comoving spatial coordinates and S{x,t) is the 
mass density contrast. The time evolution equation for 5 during the 
linear regime {6 <^ 1) reads (Peebles 1993): 

5(x, t) + 2H{t)S{x, t) = 47rG/)(t)d(x, t) + —^V^6{x, t). (1) 

a{ty 

Here, Cs is the sound speed, a = {1 -\- z)^^ is the scale factor which 
describes the expansion and H{t) = Ho[n,n{l+zf+^AV^^- The second 
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term on the left hand side of the above equation represents the effect 
of cosmological expansion. This, together with tlic pressure support 
(second term on the right hand side) acts against the growth of the 
perturbation due to the gravitational collapse (first term on the right 
hand side). The pressure in the baryonic gas is essentially provided by 
collisions, while in the collisionless dark matter component the pressure 
support arises from the readjustment of the particle orbits. The above 
equation can be used also to follow separately the evolution of the 
different components of a multi-component medium. In this case, Cg 
would be the velocity of the perturbed component (which provides 
the pressure support) and p would be the density of the component 
which is most dominant gravitationally (as it drives the collapse of 
the perturbation). The equation has two independent solutions, one of 
which grows in time and governs the formation of structTircs. The total 
density contrast at any spatial location can be described in the Fourier 
space as a superposition of modes with different wavelengths: 

j^-^Sk{t)exp(ik-x.), (2) 

where k is the comoving wave number of the Fourier series. The evo- 
lution of the single Fourier components is then given by: 

4 + 2il4 = ( 47rGp - ^ ) 5k. (3) 



a2 



This sets a critical wavelength, the Jeans length, at which the compet- 
ing pressure and gravitational forces cancel (Jeans 1928): 

For A S> A J the pressure term is negligible because the response time for 
the pressure wave is long compared to the growth time for the density 
contrast, and the zero pressure solutions apply. On the contrary, at 
A < Aj the pressure force is able to counteract gravity and the density 
contrast oscillates as a sound wave. It is conventional to introduce also 
the Jeans mass as the mass within a sphere of radius Aj/2: 

M. = ^p(^)'. (5) 



3 V 2 , 

In a perturbation with mass greater than Mj the pressure force is not 
counteracted by gravity and the structure collapses. This sets a limit 
on the scales that are able to collapse at each epoch and has a different 
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value according to the component under consideration, reflecting the 
difi^erences in the velocity of the perturbed component. 

Given the initial power spectrum of the perturbations, P{k) = |(5fcp, 
the evolution of each mode can be followed through eq. 3 and then 
integrated to recover the global spectrum at any time. The inflationary 
model predicts that P{k) oc fe", with n = 1. This value corresponds 
to a scale invariant spectrum, in which neither small nor large scales 
dominate. Although the initial power spectrum is a pure power-law, 
perturbation growth results in a modified final power spectrum. In fact, 
while on large scales the power spectrum follows a simple linear evolu- 
tion, on small scales it changes shape due to the additional non-linear 
gravitational growth of perturbations and it results in a bended spec- 
trum, P{k) oc k"'^^. The amplitude of the power spectrum, however, 
is not specified by current models of inflation and must be determined 
observationally. Note that most of the power of the fluctuation spec- 
trum of the standard CDM model is on small scales; therefore these are 
the first to become non-linear. In the following Section we will discuss 
the non-linear evolution of the dark matter halos. 

2.3. NoN Linear Evolution: Dark Matter Halos 

Because dark matter is made of coUisionless particles that interact very 
weakly with the rest of matter and with the radiation field, density 
contrast in this component can start to grow at early times. Once a 
scale has become non-linear, the linear perturbation theory described 
in the previous Section does not apply anymore. In overdense regions 
{S > 0) the self-gravity of the local mass concentration will work against 
the expansion of the universe. This region will expand at a progressively 
slower rate compared to the background universe and, eventually, will 
collapse and form a bound object. The details of this process depend 
on the initial density profile. A simple and elegant approximation to 
describe the non-linear stage of gravitational evolution has been de- 
veloped by Zel'dovich (1970). In this approach, sheetlike structures 
("pancakes") are the first non- linear structures to form from collapse 
along one of the principal axes, when gravitational instability amplifies 
density perturbations. Other structures, like filaments and knots, would 
result from simultaneous contractions along two and three axes, re- 
spectively. As the probability distribution for simultaneous contraction 
along more axes is small, in this scenario pancakes arc the dominant 
features arising from the first stages of non-linear evolution. Numerical 
simulations have been employed to test the Zeldovich approximation, 
finding that it works remarkably well at the beginning of the non-linear 
evolution. At later times, however, its predictions are not as accurate. 
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Hence, although the perturbations in Gaussian density fields are in- 
herently triaxial (e.g. Bardeen et al. 1986) and an ellipsoidal collapse 
would be more suitable to follow their evolution (e.g. Eisenstein &: 
Loeb 1995; Sheth, Mo & Tormen 2001; Sheth & Tormen 2002), the 
simplest model is the one of a spherically symmetric, constant density 
region, for which the collapse can be followed analytically. At a certain 
point the region reaches the maximum radius of expansion, then it 
turns around and starts to contract. In the absence of any symmetry 
violation, the mass would collapse into a point. However, long before 
this happens, the dark matter experiences a violent relaxation process 
and quickly reaches virial equilibrium. If we indicate with z the redshift 
at which such a condition is reached, the halo can be described in terms 
of its virial radius, r^r, circular velocity, Vc = \/ GMj r^ir , and virial 
temperature, Tyir = fj,mpV^/2kB, whose expressions are (Barkana & 
Loeb 2001): 



= 0.784 



M 







1/3 



-1/3 



l + Z 

10 



h-\-pc, (6) 



23.4 



M 



(lOS/i-^Mo 



in^ A, 11/6 /i + z\V' 



T,„> = 2 X lOM 



M 



0.6/ VlOS/i-iM« 



2 

2/3 



10 



ni, 187r2 



1/3 



km s 



1 + z 
10 



(7) 



K. (8) 



Here, f^ is the mean molecular weight, nip is the proton mass, and 
(Bryan & Norman 1998): 



IStt^ + 82(17^ - 1) - 39(a 

Q„(l + Z)3 



1) 



(9) 

Although spherical collapse captures some of the physics governing 
the formation of halos, their inner structure should be investigated 
through numerical simulations. Navarro, Prenk & White (1996, 1997 
NFW) have simulated the formation of dark matter halos of masses 
ranging from dwarfs to rich clusters, finding that their density profile 
has a universal shape, independent of the halo mass, the initial density 
fluctuation spectrum and the cosmological parameters: 



p{r) 



(r/rs)(l + r/r< 



^2' 



(11) 



where ps and are a characteristic density and radius. Usually, the 
quantity c = r^ir/rg, the concentration parameter, is introduced. As 
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Figure 3. Density profiles of halos in a simulation by Navarro at al. (2004), scaled 
to the radius, r_2, where the local logarithmic slope of the density profile takes the 
isothermal value of 2. Densities are scaled to p_2 ~ p{r-2)- 



Ps can be written in terms of c, the above equation is a one-parameter 
form. An improved fitting formula (which reproduces the more gradual 
shallowing of the density profile towards the centre, the apparent lack 
of evidence for convergence to a well-defined central power law and the 
scatter in profile for different halos) has been developed by Navarro et 
al. (2004; eqs. 4-5). Fig. 3 shows that, with proper scaling, there is little 
difference in the shape of density profile of halos of different mass. The 
halo profiles are approximately isothermal over a large range of radii, 
but they are shallower {p cx r~'^) near the center and steeper (p oc r~^) 
near the virial radius. The halo mass and the characteristic density 
are strongly correlated: low-mass halos are denser than more massive 
systems, and this reflects the higher collapse redshift of small halos. 
The existence of a universal shape is explained as the result of the 
violent relaxation process producing an equilibrium largely independent 
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of initial conditions. On the other hand, simulations which accurately 
follow the mass accretion histories of dark matter halos (Wechsler et 
al. 2002; Zhao et al. 2003), find an early phase of fast accretion, during 
which the inner structure of halos is mainly built, followed by a late 
phase of gentle accretion due to the secondary infall, which gives rise to 
an outer profile matching the universal form. Violent relaxation would 
be an active process only during the first phase. 

A universal shape of the density profile with an internal cusp, has 
been confirmed by other authors (Moore et al. 1999; Del Popolo et al. 
2000; Ghigna et al. 2000; Jing 2000; Jing & Suto 2000; Power et al. 
2003; Fukushige, Kawai & Makino 2004), although it has been noted 
that the halos simulated by NFW are intentionally selected to be in an 
equilibrium state, and thus the shape would not be, strictly speaking, 
universal. Moreover, compared with the original NFW simulations, 
some of these studies find a steeper central profile (p oc r~^, 1.2 < 
< 1.5). This result could be related to an increased resolution of 
the numerical simulations (Ghigna et al. 2000), which nowadays can 
resolve down to less than 1% of the virial radius. The simulations by 
Reed at al. (2005a), which resolve 0.5% of the virial radius, confirm the 
existence of an internal cusp with 1 < (3 < 1.5, depending on the mass 
of the halo (smaller halos have steeper cusps). 

Not all the studies converge towards a universal shape of the halo 
density profile. Rather, some authors (e.g. Kravtsov et al. 1998; Subra- 
manian, Cen & Ostriker 2000; Ricotti 2003; Cen et al. 2004) find that 
the scatter in the halo profile is large, depending on a different number 
of items. For example, the outer slope of the density profile strongly 
depends on the environment (Avila-Reese et al. 1999); the slope at 1% 
of the virial radius increases with decreasing halo mass (Jing & Suto 
2000); the slope depends on the initial fiuctuation field and cosmology 
(Nusser & Sheth 1999); the concentration parameter is not constant but 
decreases with increasing rcdshift and mass and halos in dense envi- 
ronment are more concentrated than isolated ones (Bullock et al. 2001; 
Maccio et al. 2007) ; clusters (dwarf galaxies) have steeper (shallower) 
cusps than Milky Way type galaxies with a concentration parameter 
that is a universal constant at the virialization rcdshift (Ricotti 2003). 
Some of the apparent tensions with works that find a universal shape 
of the density profile might be due to the different method used to 
analyse the simulations (Ricotti, Pontzen & Viel 2007). 

Also on the presence of a central cusp there is no general agreement 
(e.g. Stoehr 2006). For example, if the effect of gas is included in the cal- 
culation a flat core is more plausible than a central cusp. In fact, if the 
gas is not initially smoothly distributed, but is concentrated in clumps, 
dynamical friction acting on these clumps dissipates their orbital energy 
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and deposits it in the dark matter, resulting in a heating of the halo 
and the formation of a finite, non-divergent core (El-Zant, Shlosman & 
Hoffman 2001; Mashchenko, Couchman &: Wadsley 2006). A core would 
be also expected if the standard description of the post-collapse object 
were replaced by a more realistic truncated, non-singular, isothermal 
sphere (TIS) in virial and hydrostatic equilibrium (Shapiro, Iliev &; 
Raga 1999). An alternative mechanism that could prevent a cusp form- 
ing is the following. The build up of dark matter halos by merging 
satellites inevitably leads to an inner cusp with /? > 1. A flatter core 
with P < 1 exerts on each satellite tidal compression which prevents 
deposit of stripped satellite material in this region. An inner cusp is 
expected as long as enough satellite material is captured by the inner 
halo. If, instead, satellites are disrupted in the outer halo, a flatter core 
is maintained (Dckcl, Dcvor & Hetzroni 2003). A flat core could be 
also explained with the existence of a dark spheroid of baryons of mass 
comparable with that of the cold dark matter component (Burkert Sz 
Silk 1997). 

The density of the dark matter halos and their spatial distribution 
are better known. Two main methods have emerged to evaluate them: 
numerical computations that solve the equations of gravitational col- 
lapse, and analytical techniques that approximate these results with 
simple one-dimensional functions. While only numerical simulations 
give the spatial distribution of halos, analytical techniques are ex- 
tremely useful as they are much faster and allow the analysis of a wide 
range of parameters. The most commonly applied method of this type 
was first developed by Press Sz Schechter (1974). In their approach, the 
abundance of halos at a redshift z is determined from the linear density 
field by applying a model of spherical collapse to associate peaks in the 
field with virialized objects in a full non-linear treatment. This simple 
model, later refined by Bond et al. (1991) and Lacey Sz Cole (1994), 
has had great success in describing the formation of structures and 
reproducing the numerical results. The method provides the comoving 
number density of halos, dn, with mass between M and M + dM as 
(Peebles 1993): 



where po is the present mean mass density, a is the standard deviation 
of the density contrast smoothed through a certain window and fc is 
the minimum number of standard deviations of a collapsed fluctuation; 
Fig. 4 shows an application of the above equation. 

While the model accurately reproduces the number density of the 
dark matter halos, it does not predict their spatial distribution, which 




(12) 
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Figure 4- Each curve indicates the variation with redshift of the comoving number 
density of dark matter halos with masses exceeding a specific value M in the stan- 
dard ACDM model. The label on each curve indicates the corresponding value of 
log(M/M©) (see Mo & White 2002 for details). 



is important to study structure formation. For this reason, Porciani et 
al. (1998) and Scannapieco & Barkana (2002) have developed analytical 
techniques to derive halo correlation functions and non-linear biasing, 
by studying the joint statistics of dark matter halos forming at two 
points. An alternative approach to the analytical determination of the 
abundance and spatial distribution of dark matter halos has been de- 
veloped by Sheth &; Tormen (2002). This is known as the excursion set 
approach and can be applied in the case of both a spherical (to repro- 
duce the standard Press-Schechter results) and an elliptical collapse, 
which seems to be in better agreement with numerical simulations. 
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2.4. Formation of protogalaxies 

In contrast to dark matter, as long as the gas is fully ionized, the ra- 
diation drag on free electrons prevents the formation of gravitationally 
bound systems. When recombination (defined as the time at which 
the electron fraction has dropped to 0.1) occurs Qit Zf^Q. ~ 

1250 the 

primeval plasma combines into neutral atomic hydrogen. The resid- 
ual ionization of the cosmic gas keeps its temperature locked to the 
CMB temperature through Compton scattering, down to a rcdshift 
1 + zt 1000(^^6/1^)2/5 (Peebles 1993). Following recombination, the 
nearly neutral matter decouples from the radiation and perturbations in 
this component are finally able to grow in the pre-existing dark matter 
halo potential wells and eventually form the first bound objects. The 
process leading to the virialization of the gaseous component of matter 
is similar to the dark matter one. In this case, during the contraction 
following the turnaround, the gas develops shocks and gets reheated to 
a temperature at which pressure support can prevent further collapse. 

The mass of these first bound objects can be derived from eqs. 4- 
5, where Cg is the velocity of the baryonic gas. In particular, at 2; > 
zt the Jeans mass is time-independent, while at z < zt, when the 
gas temperature declines adiabatically, Mj decreases with decreasing 
redshift: 

— ..o3(^)-"^(iM!)-"'(i±.)-M.. as) 

As the determination of the Jeans mass is based on a perturbative 
approach, it can only describe the initial phase of the collapse. Mj is 
simply a limit given by the linear theory to the minimum mass that is 
able to collapse. It is worth noticing that the Jeans mass represents only 
a necessary but not sufficient condition for collapse: in addition, one has 
to require that the atomic/molecular cooling time is shorter than the 
Hubble time. The dynamical and chemical evolution of collapse must 
then be derived in detail to assess if the object can actually become 
a protogalaxy. In the following we will discuss the minimum mass of 
protogalaxies. 

In the standard cosmological hierarchical scenario for structure for- 
mation, the objects which collapse first are predicted to have masses 
corresponding to virial temperatures T^i^ < 10^ K. Once the gas has 
virialized in the potential wells of dark matter halos, additional cooling 
is required to further collapse the gas and form stars. For a gas of 
primordial composition at such low temperatures the main coolant is 
molecular hydrogen. We define Pop III objects those for which H2 cool- 
ing is required for collapse (note that in the literature Pop III objects 
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are sometimes also referred to as minihalos) . After an H2 molecule gets 
rotationally or vibrationally excited through a collision with an H atom 
or another H2 molecule, a radiative de-excitation leads to cooling of the 
gas. This mechanism has been investigated in great detail by several au- 
thors (Lepp & Shull 1984; Hollenbach & McKee 1989; Martin, Schwarz 
& Mandy 1996; Galli & Palla 1998) and has been applied to the study 
of the formation of primordial small mass objects (e.g. Haiman, Thoul 
& Loeb 1996; Haiman, Rees k Loeb 1996; Abel et al. 1997; Tegmark 
et al. 1997; Abel et al. 1998; Omukai & Nishi 1999; Abel, Bryan & 
Norman 2000; Yoshida et al. 2003; Maio et al. 2007). Primordial H2 
forms with a fractional abundance of 10"^ at redshifts ^ 400 via the 
II2 formation channel. At redshifts ^ 110, when the CMB radiation 
intensity becomes weak enough to allow for significant formation of H~ 
ions, even more H2 molecules can be formed: 

H e~ Il~ + hi' 
H- + H ^ H2 + e- 

Due to the lack of molecular data, it has unfortunately not been possible 
to follow the details of the H2' chemistry as its level distribution de- 
couples from the CMB. Assuming the rotational and vibrational states 
of 112" to be in equilibrium with the CMB, the H2' photo-dissociation 
rate is much larger than the one obtained by considering only photo- 
dissociations out of the ground state. Conservatively, one concludes that 
these two limits constrain the H2 fraction to be in the range 10-^-10"''. 
In fact, both limits have been used in the literature (e.g. Lepp &: Shull 
1984; Palla, Galh & Silk 1995; Haiman, Thoul & Loeb 1996; Tegmark 
et al. 1997). If we assume that the H~ channel for H2 formation is the 
dominant mechanism, i.e. that the H2' photo-dissociation rate at high 
redshift is close to its equilibrium value, this leads to a typical primor- 
dial H2 fraction of ///j ~ 2 x 10~^ (Anninos &; Norman 1996) which 
is found for model universes that satisfy the constraint Qi,h'^ = 0.019. 
Hirata & Padmanabhan (2006) have recently revised the formation of 
primordial H2 including (i) the effect of spectral distortions due to 
radiation emitted during hydrogen and helium recombination (which 
destroys H~ at z > 70) and (ii) a detailed treatment of the H^ level 
population (which suppresses the efficiency of the H2' formation chan- 
nel). As a result, they find the slightly lower primordial value of 6x 10~^. 
It should be noted that the above effects do not influence the H2 
formation and cooling at lower redshift and in overdense conditions. 
In any case, this primordial fraction is usually lower than the one 
required for the formation of Pop III objects, but during the collapse 
the molecular hydrogen content can reach high enough values to trigger 
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Figure 5. Minimum mass able to cool and collapse as a function of redshift in 
the formulation of Tegmark et al. (1997) (solid line), Fuller & Couchman (2000) 
(dashed) and Shchekinov & Vasiliev (2004) (dashed-dotted) . The latter curves are 
derived for two different values of the rate of ionizing photon production by UHECR, 
e — (0.1, 3)/(l + z) for the upper and lower line, respectively. See original papers for 
details. 



star formation. Thus, the fate of a virialized lump depends crucially on 
its ability to rapidly increase its H2 content during the collapse phase. 

Haiman, Thoul &: Loeb (1996), have followed the growth of spher- 
ical perturbations into the non-linear regime, using a one-dimensional 
spherical Lagrangian hydrodynamics code to follow the dynamical, 
thermal and non-equilibrium chemical evolution of the gas. They note 
that H2 cooling becomes important only after virialization. Interest- 
ingly, though, they find that shell crossing by dark matter allows bary- 
onic objects with masses well below the Jeans mass to form. A similar 
result is obtained by Haiman &: Loeb (1997). Tegmark et al. (1997) 
have addressed the same question by calculating the evolution of the 
H2 abundance for different halo masses and initial conditions for a 
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standard CDM cosmology. They use an analytical approach to the 
problem and find that the minimum baryonic mass M5 is strongly 
redshift dependent, dropping from 10^ Mq at z ~ 15 to 5 x 10^ Mq at 
z « 100, as molecular hydrogen cooling becomes effective (see Fig. 5). 
They conclude that if the prevailing conditions are such that a molecu- 
lar hydrogen fraction ~ 5 x 10~^ is produced, then the lump will cool, 
fragment and eventually form stars. This criterion is met only by larger 
halos, implying that for each virialization redshift there is a critical 
mass. Merit, such that protogalaxies with total mass M > Merit will 
be able to collapse and form stars and those with M < Merit will fail. 
In a subsequent paper Nishi &; Susa (1999) claim that the primordial 
number fraction of H2 as calculated by Tegmark et al. (1997) is overes- 
timated by about two orders of magnitude, because the destruction rate 
of by CMB radiation at high redshift is underestimated (Galli & 
Palla 1998). Thus, since their primordial value ~ 10~^ is comparable to 
the one necessary to cool, their cooling criterion is not generally reliable. 
Fuller & Couchman (2000) have included the chemical rate network 
responsible for the formation of H2 in an N-body hydrodynamical code, 
finding a value for Merit considerably lower than the one derived by 
Tegmark et al. (1997) (see Fig. 5), since they have used a different 
molecular hydrogen cooling function and have included different chem- 
ical reactions, but they nevertheless find the same critical molecular 
hydrogen fraction for the collapse of fH2 ~ 5 x 10"^ . Most studies (e.g. 
Abel, Bryan & Norman 2000; Machacek, Bryan k Abel 2001; Reed 
et al. 2005b; O'Shea & Norman 2007) agree that the minimum mass 
allowed to collapse is as low as 10^ Mq. 

The results of recent quantum mechanical calculations of cross- 
sections for rotational transitions within the vibrational ground state of 
HD enabled the contribution of HD to the primordial cooling function 
to be reliably calculated for the first time. In spite of its low abundance 
(10"^) relative to H2, HD can be as important in the thermal balance 
of the primordial gas. In particular, the permanent dipole moment and 
lower rotational constant of HD favor cooling by HD at high densities 
and low kinetic temperatures (Flower 2000; Flower et al. 2000; Galli 
& Palla 2002), possibly allowing even smaller objects to collapse. If 
also vibrational transitions are taking into account, the cooling rate 
is increased (Lipovka, Nuiiez-Lopez Sz Avila-Reese 2005). In Fig. 6 a 
comparison between the H2 and HD cooling rates is shown. 

A considerably smaller characteristic mass (a factor ~ 10 below 
the one predicted by Tegmark et al. 1997) of first luminous objects 
could also result from the presence of ultra-high energy cosmic rays 
(UHECRs), if they formed in the so-called top-to-bottom scenario from 
decaying superheavy dark matter particles with masses > 10^^ GeV. 
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Figure 6. Left panel: cooling rate per HD molecule calculated for an ortho/para H2 
density ratio of 1 and an H/H2 abundance ratio of 1. Right panel: cooling rate per 
H2 moloculo. The parameters are the same as in the left panel. For details refer to 
Flower et al. (2000) and Le Bourlot, Pineau des Forets & Flower (1999). 



In fact, the small additional electron fraction produced by UHECR 
photons rcsuhs in more efficient H2 and HD molecule formation and 
cooling (Shchekinov & Vasiliev 2004, 2006; see Fig. 5). Also the presence 
of an early population of cosmic rays (that would produce free electrons; 
Jasche, Ciardi & Enfilin 2007; Stacy & Bronim 2007) or of sterile neutri- 
nos with masses of several keV (that would decay in X-rays and produce 
free electrons; Biermann &; Kusenko 2006) could boost the formation of 
the above molecules. On the other hand, Ripamonti, Mapelli k, Ferrara 
(2007b) find that moderately massive dark matter particle (as sterile 
neutrinos and light dark matter) decays and annihilations favors the 
formation of H2 and HD, but at the same time induce heating of the 
gas. The net effect in this case is a slight increase of yicHt- 

It should be noted though that all the above estimates do not take 
into account the large uncertainties on the reaction rates that govern 
the H2 formation and that can have a non negligible impact on the 
prediction for early structure formation (e.g. Glover, Savin &; Jappsen 
2006). 
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Figure 7. Matter power spectra (from Sugiyama, Zaroubi & Silk 2003). The WMAP 
best fitted power law result is indicated by PLACDM, and it is compared with 
the WMAP Running Spectral Index model (RSI) with errors combined with different 
additional data (hatched regions); the curves are isocurvature spectra with different 
prescriptions. See original paper for details. 

2.5. Key Observations 
2.5.1. Cosmological Model 

The main observational evidence in support of the inflationary scenario 
comes from observations of the CMB, showing temperature anisotropics 
on large angular scales (e.g. Hu et al. 2001), which are interpreted as the 
results of the primordial density fluctuations at the decoupling epoch. 
The same measurements confirm that the fluctuations are adiabatic 
and the value of the spectral index is the one predicted by the infla- 
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tionary model (Balbi et al. 2000, 2001; Pryke et al. 2002; Sievers et al. 
2003), although there is an indication for a possible running spectral 
index (e.g. Spergel et al. 2003, 2007), with n > 1 on large scales and 
n < 1 on small scales (see Fig. 7 for possible power spectra). This is 
not confirmed though when the constraints from the recent Lya forest 
and galaxy bias analysis of the SDSS^^ (Sloan Digital Sky Survey) are 
combined with previous constraints from SDSS galaxy clustering, the 
latest supernovae and the WMAP results (Seljak et al. 2005; Seljak, Slosar 
&; McDonald 2006). In fact, while CMB anisotropies, galaxy clustering 
from 2dF and SDSS, supernovae and weak Icnsing trace large scale 
structures (scales larger than 10 Mpc), tracers on smaller scale are 
needed to get correct information on the running spectral index. For 
the same reason, a tension is present also in the cjg value measured with 
or without the Lyo: forest data, giving a higher value in the latter case 
(Seljak, Slosar & McDonald 2006; Viel, Haehnelt & Lewis 2006). More 
specifically, if only the 3-rd year WMAP data are used, ag = 0.76 ± 0.05, 
but the value increases with additional constraints (Spergel et al. 2007), 
in particular those from the Lya forest (see Fig. 8; Seljak, Slosar &; 
McDonald 2006; Viel, Haehnelt Sz Lewis 2006). In addition, analyses of 
the WMAP observations have reported evidence for a positive detection 
of non-Gaussian features (e.g. Chiang et al. 2003; Rath, Schuecker &; 
Banday 2007; Eriksen et al. 2007; McEwen et al. 2008). 

As discussed in the previous Sections, the amplitude of the power 
spectrum is not predicted by the model and has to be set obscrva- 
tionally. Historically, this has been derived by comparing the observed 
CMB quadrupole anisotropies with theoretical ones (e.g. Wright et 
al. 1992; Efstathiou, Bond & White 1992) or from the abundance of 
galaxy clusters (see Sec. 2.1). More recently, this method has been 
complemented by measuring the amplitude of the power spectrum 
from available large galaxy redshift surveys (e.g. Lahav et al. 2002). 
The main problem related to the latter method is the attempt to 
estimate the distribution of matter using galaxies as tracers. In fact, 
as the galaxies form preferentially in high density regions of the mass 
distribution, they are more strongly correlated than the underlying 
distribution and an additional factor, the bias, should be specified. It 
has been found from the 2dF Galaxy Redshift Survey that, at least on 
large scales (5 — 30^~^ Mpc), optically selected galaxies do indeed trace 
the underlying mass distribution (Verde et al. 2002). 

An alternative method to measure P{k) is based on observations of 
the high-redshift Lya forest. The method is motivated by the physical 



http:/ /www.sdss.org/ 
http://magnum.anu.edu.au/ TDFgg/ 
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Figure 8. la and 2a likelihoods for erg and the spectral index (in this figure indicated 
with Tis), marginalized over all other parameters. Left panel: constraints for 1-yr WMAP 
only (green), the LUQAS (Kim et al. 2004) + Croft et al. (2002) data set as analyzed 
by Viel, Haehnelt & Springel (2004, VHS) (cyan) and the SDSS Lya forest data of 
McDonald et al. (2005) (blue). The thick dashed white empty contours refer to 1-yr 
WMAP + VHS, while the sohd blue contours are for 1-yr WMAP + SDSS. Right panel: 
as in the left panel but for the 3-yr WMAP. 



picture that has emerged from hydro dynamical cosmological simula- 
tions and related semi-analytical models, in which typical Lya forest 
lines arise in a diffuse IGM. The thermal state of this low-density 
gas is governed by simple physical processes, which lead to a tight 
correlation between the Lya optical depth, r, and the underlying mat- 
ter density (e.g. Bi & Davidsen 1997; Hui, Gnedin & Zhang 1997). 
Thus, as the transmitted flux in a QSO spectrum \s F = e~'^ , one can 
extract information about the underlying mass density field from the 
observed flux distribution (Croft et al. 1998, 1999, 2002; Viel, Haehnelt 
& Springel 2004; McDonald et al. 2005). These works confirm a basic 
prediction of the inflationary CDM scenario, an approximately scale 
invariant spectrum of primeval fluctuations modulated by a transfer 
function that bends it toward A;"~^ on small scales. A possible caveat 
on these applications is that radiative transfer effects (Abel &: Haen- 
helt 1999; Bolton, Meiksin & White 2004) are not included in current 
numerical simulations of the Lya forest. Such effects might (i) blur 
the polytropic temperature-density relation derived from the simula- 
tions and usually adopted, hence introducing a non-negligible error on 
the predicted cosmological spectrum, and (ii) produce an inversion of 
the temperature-density relation in regions of low density (voids), as 
pointed out by Bolton et al. (2008). 
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Figure 9. Left panel: Forecasts of the la uncertainties in estimates of the convergence 
power spectrum for a LOFAR-lyr (blue) and L0FAR-3yr (red) type experiment of 21 cm 
weak lensing. The sohd straight lines give the noise per mode, while the dotted curve 
is the underlying model power spectrum. Right panel: As the left panel but for a 
SKA type experiment. For details see Metcalf & White (2008). 



It has been proposed (Loeb & Zaldarriaga 2004; Ali, Bharadwaj & 
Panday 2005; Barkana & Loeb 2005a; Bharadwaj & Ah 2005; Pihepich, 
Porciani & Matarrese 2007; Lewis & Chalhnor 2007; Cooray, Li &; 
Melchiorri 2008; Mao et al. 2008) that a wealth of information about 
the initial density fluctuations (e.g. the presence of a running spectral 
index or deviations from gaussianity) could come from the 21 cm line in 
absorption against the CMB, from z ~ 200 down to lower redshift (as 
long as the effects of reionization are not dominant; Bowman, Morales 
& Hewitt 2007) or from the application of the Alcock-Paczyhski test to 
the 21 cm line (Nusser 2005b). In these cases, measurements from up- 
coming radio telescopes in combination with CMB data, would produce 
constraints on cosmological parameters tighter that those obtained by 
CMB data alone (McQuinn et al. 2006). To make precision cosmology 
though, particular care should be used in the modeling of the spin 
temperature (Hirata & Sigurdson 2007). 

Another 21 cm related experiment that could be used to reconstruct 
matter density fluctuations and to obtain constraints on cosmological 
parameter is weak lensing, using 21 cm emission as background source 
(see Fig. 9). Zahn & Zaldarriaga (2006) show that the lensing effect is 
lower compared to the case in which CMB fluctuations are used, but it 
has the advantage of retrieving information from a range of redshifts. 
In any case, 21 cm lensing experiments could provide constraints on 
cosmological parameters better than those obtained from comparably 
sized surveys of galaxy lensing (Metcalf & White 2007; Hilbert, Metcalf 
& White 2007; Lu & Pen 2008; Metcalf & White 2008). 
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2.5.2. Dark Matter Halo Density Profile 

To prove or disprove the existence of a universal density profile for 
dark matter halos, accurate and extensive observations are needed. A 
method largely applied in the past to the study of the density profiles 
relies on HI rotation curves, although the spatial resolution of these 
observations is rarely good enough to set meaningful constraints. This 
is true both for the rotation curves of Low Surface Brightness Galaxies 
(LSBGs) and late-type dwarf galaxies (e.g. van den Bosch et al. 2000; 
Borriello &; Salucci 2001; de Blok et al. 2001; van den Bosch & Swaters 
2001; de Blok, Bosma & McGaugh 2003; Kleyna et al. 2003), which are 
sometimes consistent with both a constant density core and an inner 
cusp, indicating that HI rotation curves do not have enough resolution 
to discriminate between the two models. Only in a few cases it has 
been possible to derive a meaningful value of /? (/3 is the power-index 
of the density profile in the central regions), giving 0.55 < (3 < 1.26 
at a 99.73% confidence level for an LSBG, which would be consistent 
with the NFW slope, and a (3 < 0.5 for two dwarf galaxies at the same 
confidence level (van den Bosch et al. 2000) , more consistent with a flat 
core, as the nearby Ursa Minor dwarf spheroidal measured by Kleyna 
et al. (2003) or the sample of LBGs analysed by de Blok, Bosma &; 
McGaugh (2003), who flnd an average (3 = 0.2. Hq rotation curves 
have a higher spatial resolution and rise more steeply in the inner parts 
than the HI rotation curves (Swaters, Madore Sz Trewhella 2000). Also 
optical rotation curves have been used, indicating that the NFW profile 
provides a good fit to 66% of the 400 galaxies in the sample analyzed 
by Jimenez, Verde Sz Oh (2003), while 68% of the galaxies are well 
fitted by an isothermal profile with a core. Based on optical and radio 
rotation curves, Persic, Salucci & Stel (1996) have confirmed that spiral 
galaxies have a universal rotation curve, characterized by one single free 
parameter, the (I-band) luminosity: low-luminosity spiral galaxies show 
rising rotation curves out to the optical radius, while high-luminosity 
ones are fiat or even decreasing. These findings are reiterated in Salucci 
et al. (2007). On the other hand, recent HST observations have revealed 
that elliptical galaxies have cusps which continue toward the center 
until the resolution limit. Bright elliptical galaxies have a shallow cuspy 
core with 0.5 < /3 < 1, while faint ones have /3 ~ 2 (Merritt &i Fridman 
1996). 

Alternatively, the density profile of dark matter halos can be inves- 
tigated through gravitational lensing experiments. It has been argued 
that radially distorted, gravitationally lensed images of background 
sources in galaxy clusters, the so-called radial arcs, require a flat core 
in the cluster density proflle. This would be consistent with the high 
resolution map of a cluster in which a very smooth, symmetric and 
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non-singular core is observed (Tyson, Kochanski & Dell' Antonio 1998). 
Nevertheless, Bartelmann (1996) shows that the NFW profile can pro- 
duce radial arcs despite its central singularity, as indicated by more 
recent observational data (Oguri, Taruya &: Suto 2001). Observations 
of radial arcs though are so scarce that larger samples are needed for 
a more significant study. 

Finally, as integral measures of weak gravitational lensing by dark 
matter halos, like the aperture mass, are sensitive to the density profile, 
these can be used to discriminate between an isothermal and a NFW 
profile. In particular, as the halo mass range probed by the aperture 
mass is much wider for a NFW profile, counts of halos with significant 
weak lensing signal are powerful discriminators (Bartelmann, King &; 
Schneider 2001). 

In summary, available observations are not sufficient to significantly 
constrain the profile of dark matter halos. 

2.5.3. Formation of Protogalaxies 

As primordial gas condenses within dark matter potential wells, form- 
ing luminous galaxies, it emits copious radiation. The energy radiated 
in this process is comparable to the gravitational binding energy of 
the baryons. Most of this cooling radiation is emitted by gas with 
T < 20,000 K. As a consequence, roughly 50% of it emerges in the 
Lya line. As such emission is predominantly emitted from the outskirts 
of the collapsing system, it is less likely to be extinguished by dust (if 
present) than the more deeply embedded stellar light. This radiation 
has been advocated (Haiman, Spaans & Quataert 1999; Fardal et al. 
2001) to explain the large (si 100 kpc), luminous {L !=a 10^ erg s"-*^) 
"blobs" of Lya emission found in narrow-band surveys of z ~ 3 proto- 
clusters (Steidel et al. 2000) and similar arguments could be applied to 
higher redshift objects. 

In addition to Lya radiation, molecules could also produce a de- 
tectable signal. Flower & Pineau des Forets (2001) calculated the spec- 
trum of H2 that is produced by collisional excitation in the shock waves 
generated when the speed of collapse becomes locally supersonic, find- 
ing that the rotational transitions within the vibrational ground state 
might be observable as inhomogeneities in the CMB. The H2 cooling 
emission of forming galaxies has been studied by Omukai & Kitayama 
(2003), who discuss its observability (in particular the rotational lines 
0-0S(3) at 9.7 fim and O-OS(l) at 17 fiui) through SAFIR^*^ (Single 
Aperture Far-Infrared observatory) and by Mizusawa, Omukai & Nishi 
(2005), who find that both rovibrational and pure rotational H2 lines 



http://safir.jpl.nasa.gov/index.asp 
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can be detected by SPICA^^ (SPace Infrared telescope for Cosmology 
and Astrophysics), if emitted by metal- free clouds aX z < 10. 

3. First Stars 

As described in the previous Section, hierarchical models of structure 
formation predict that the first collapsed objects have small masses and 
arc of primordial composition and thus the formation and cooling of 
these objects is mainly governed by the molecular hydrogen chemistry. 
As the collapse proceeds, the gas density increases and the first stars 
are likely to form. However, the primordial star formation process and 
its final products are presently quite unknown. This largely depends on 
our persisting ignorance of the fragmentation process and on its rela- 
tionship with the thermodynamical conditions of the gas. Despite these 
uncertainties, it is presently accepted that the first stars, being formed 
out of a gas of primordial composition, are metal-free {Pop III stars). 
In addition to this theoretical justification, the existence of Pop III 
stars have been historically invoked for many different reasons. Among 
other issues (see e.g. Larson 1998 and Chiosi 2000 for more complete 
reviews), Pop III stars could help explaining: 

— the gap between the Big Bang nucleoshynthesis (BBN) metal abun- 
dances {Z 10"^*^ — 10~^^) and those observed in the lowest 
metallicity Pop II stars {Z ^ 10"^ - 10"^); 

— the G-dwarf problem, i.e. the paucity of metal-poor stars in the 
solar neighborhood relative to predictions from simple models of 
chemical evolution (van den Bergh 1962; Norris, Peterson & Beers 
1993; Primas, Molaro & Castelh 1994; Sneden et al. 1994); 

— the oxygen anomaly or enhancement of a-elements in galactic 
metal-poor stars (Sneden, Lambert & Whitaker 1979) and the 
existence of extremely metal-poor stars showing s-process elements 
in their envelopes (Truran 1980); 

— the missing mass in clusters of galaxies and galactic halos (White 
& Rees 1978) formed by their dark remnants; 

— the reionization of the universe and the starting engine for the 
formation of the first galaxies; 

— the constant lithium abundance of metal-poor halo stars (the Spite 
Plateau) which is lower that the predictions from the BBN (Piau 
et al. 2006; Asplund et al. 2006); 

^"^ http://www.ir.isas.ac.jp/SPICA/ 
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— the contaminants of the intergalactic medium as inferred from 
metalhc absorption hnes in the Lya forest seen in quasar hght; 

— the cosmological hehum abundance (Talbot & Arnett 1971; Marigo 
et al. 2003; Salvaterra & Ferrara 2003b); 

— the formation of massive black holes. 

In the following Section we will discuss the physical processes governing 
the formation and evolution of the first stars. 

3.1. Formation Process 

To study the formation of the first stars it is crucial to understand the 
physics of the cooling in the early universe, primarily the number of 
channels available for the gas to cool, their efficiency (Rees & Ostriker 

1977) and the physical processes that set the scales of fragment masses 
and hence the stellar mass spectrum. For a recent review see McKee Sz 
Tan (2008). 

In general (e.g. Schneider et al. 2002), cooling is efficient when 

^cool ^ %i where tcool = 3nkT/2A{n,T) is the cooling time, iff = 
(37r/32G/9)^/^ is the free-fall time, n (p) is the gas number (mass) 
density and A(n,T) is the net radiative cooling rate (in units of erg 
cm~^s~^). This efficiency criterium implies that the energy deposited 
by gravitational contraction cannot balance the radiative losses; as a 
consequence, temperature decreases with increasing density. When the 
above condition is satisfied, the cloud cools and possibly fragments. At 
any given time, fragments form on a scale, Rp, that is small enough 
to ensure pressure equilibrium at the corresponding temperature, i.e. 
Rp « Aj oc Cgtff oc n'l'^~^ where the sound speed Cs = (JIT / iiY^'^ , 
T oc rC~^, IZ is the universal gas constant and 7 is the polytropic 
index. Since Cg varies on the cooling timescale, the corresponding Rp 
becomes smaller as T decreases. Similarly, the corresponding fragment 
mass is the Jeans mass, 

Mp oc nRl oc n'1"'/^+(^-'^) , (14) 

with T] = 2 for filaments and r] = 3 for spherical fragments (Spitzer 

1978) . This hierarchical fragmentation process comes to an end when 
cooling becomes inefficient because (i) the critical density for Local 
Thermodynamical Equilibrium (LTE) is reached or (ii) the gas becomes 
optically thick to cooling radiation. In both cases, ^cooi becomes larger 
than tff. At this stage, the temperature cannot decrease any further 
and it cither remains constant (if energy deposition by gravitational 
contraction is exactly balanced by radiative losses) or increases. The 
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necessary condition to stop fragmentation and start gravitational con- 
traction within cacli fragment is that the Jeans mass does not decrease 
any further, thus favoring fragmentation into sub-clumps. From eq. 14, 
this implies the condition 7 ^ 2(r/ — l)/ry, which translates into 7 4/3 
for a spherical fragment and 7 ^ 1 for a filament. Thus, a filament is 
marginally stable and contracts quasi-statically when tcooi ~ is, and 
the gas becomes isothermal. Finally, when tcooi ^ iff or the fragments 
become optically thick to cooling radiation, the temperature increases 
as the contraction proceeds adiabatically. 

The very first generation of stars must have formed out of prob- 
ably unmagnetized, pure H/He gas, since heavy elements can only 
be produced in the interior of stars. These characteristics render the 
primordial star formation problem very different from the prcscnt-day 
case, and lead to a significant simplification of the relevant physics. 
Nevertheless, the complexity and interactions of the hydrodynamical, 
chemical and radiative processes, have forced many early studies to use 
the steady state shock assumption, a spherical or highly idealized col- 
lapse model. Only recently have more realistic models been proposed. 
The first self-consistent three-dimensional cosmological hydrodynami- 
cal simulations were presented by Abel et al. (1998). Because of the 
limited spatial resolution, it was not possible to study the collapse to 
stellar densities and address the nature of the first objects. A higher 
resolution simulation has been performed by Bromm, Coppi &z Larson 
(1999). They investigate, by means of a Smooth Particle Hydrodynam- 
ics (SPH) simulation, the evolution of an isolated 3-0" peak of mass 
2 X 10^ Mq collapsing at z 30. They find that the gas dissipatively 
settles into a rotationally supported disk, which eventually fragments 
in high density clumps (n ~ 10^ cm^'^) of mass M^ump ~ 10^ — 10^ 
M0. Mass accretion and merging could subsequently raise the clump 
masses up to !=a 10^ Mq. Changing initial conditions does not affect the 
disk formation and fragmentation, although the accretion and merg- 
ing history can change (Bromm, Coppi &; Larson 2002). Bromm &; 
Loeb (2004) extended the results of the simulations by Bromm and 
collaborators down to a scale of ?alOO AU, following the collapse of 
the cloud to a fully-molecular protostellar core and the subsequent 
accretion phase. Although improving on previous simulations, radiative 
and mechanical feedback from the protostar on the accretion flow are 
not treated self-consistently. 

Similar results have been found by Abel, Bryan & Norman (2000). 
Their simulation is based on a 3-D Adaptive Mesh Refinement (AMR), 
covering comoving scales from 128 kpc down to 1 pc and follows the 
collapse of primordial molecular clouds and their subsequent fragmen- 
tation within a cosmologically representative volume. The authors find 
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Figure 10. Projected distribution for a first star cosmological simulation at z « 19. 
The physical side length is indicated in each panel. See Yoshida et al. (2006) for 
details. 



that these molecular clouds, of mass Mdoud ~ 10^ Mq, will eventually 
form a cold (~ 200 K) region at the center of the halo. Within this cold 
region a quasi-hydrostatic contracting clump with mass ~ 200 Mq 
and density n > 10^ cm~^ (and as high as ~ 10*^ cm~^) is found. 
The authors speculate that more than one clump could be found in 
the same object and that their number is likely to be proportional to 
the total amount of cooled gas. Clearly, the ability of these clumps 
to form stars is closely related to the strength of feedback processes 
(see Sec. 4). In addition, they suggest that the formation of a disk and 
its fragmentation, found by Bromm, Coppi & Larson (1999, 2002), is 
related to their highly idealized initial conditions (top-hat spheres that 
initially rotate as solid bodies on which smaller density fluctuations 
are imposed). But, although the process that leads to the formation of 
dense clumps in the two simulations is different, the final result is quite 
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similar, with very dense, massive protostellar clumps out of which the 
first star will eventually form. See Fig. 10 for a typical simulation of 
Pop III star formation. 

It should be noted that the fragmentation and the thermal evolution 
depend strongly on the geometry of the clouds. For filamentary clouds 
the gravitational contraction is slower and the compressional heating 
rate is lower than for spherical ones. Hence, for the filamentary cloud 
the isothermality breaks down at higher density, leading to a smaller 
mass of the final clumps. Nakamura & Umemura (2001) have studied 
the collapse and fragmentation of a filamentary cloud using a two- 
dimensional hydrodynamic simulation, finding that the mass of the 
resulting clumps depends on the initial density of the collapsing cloud: 
filaments with low initial density (n < 10^ cm""') tend to fragment into 
dense clumps before the cloud becomes optically thick to the H2 lines 
and the H2 cooling with the three-body reaction is effective. This would 
result in clumps much more massive than some tens Mq. However, 
relatively dense filaments result in clumps of a few Mq (the result is in 
practice somewhat ambiguous as proper "initial conditions" are rarely 
known). The more effective H2 cooling with the three-body reaction 
allows the filament to contract up to n ~ 10^^ cm^'', when it becomes 
optically thick to H2 lines and the radial contraction stops. 

As the first stars form in regions of the universe corresponding to 
peaks of the large-density field, a significant contribution to the local 
overdcnsity comes from large wavelength fluctuations. Thus, simula- 
tions that take into account this effect are needed to properly study the 
formation of the first stars. Gao et al. (2007) have run a variety of SPH 
simulations (from cosmological scales down to densities at which the gas 
is no more optically thin to molecular hydrogen lines, i.e. ~ lO^'^ cm~^) 
to investigate whether there is a significant scatter in the properties of 
primordial star-forming clouds and whether they depend on redshift. 
They find that the formation path of primordial stars is similar in 
all simulations, independently on redshift and environment. The only 
difference is that the timescale of the formation is much shorter at 
high redshift. Moreover, the star-forming clouds exhibit a variety of 
morphologies and different accretion rates. In general though, at the 
end of the simulation, none of the clouds have fragmented and stars of 
more than a few tens of solar masses can be expected, although nothing 
more precise can be said on the mass of the first stars. 

It should be noted that SPH simulations that do not resolve the 
Jeans mass at any time might suffer from the problem called artificial 
fragmentation, which, if not properly addressed, can induce errors in 
the determination of the mass of the stars and their Initial Mass Func- 
tion (IMF). A relatively new technique, the particle splitting (Kitsionas 
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& Whitworth 2002), has been used by Martel, Evans & Shapiro (2006) 
to investigate the artificial fragmentation of a molecular cloud. 

O'Shea & Norman (2007) have run AMR simulations with the code 
ENZO to study systematic effects in the formation of primordial proto- 
stellar cores. As previous studies, they find no evidence of fragmentation 
in any of their simulations, suggesting that Pop III stars forming in 
halos of ~ 10^ — 10^ Mq form in isolation. The minimum mass of halos 
hosting Pop III stars (~ 10^ M©), as well as other bulk halo properties, 
is independent on redshift. 

Small mass clumps can be obtained when the collapse of a cloud is 
driven by HD or atomic cooling rather than molecular cooling. For ex- 
ample, Uehara Sz Inutsuka (2000) find that, including HD cooling, frag- 
mentation of primordial gas can produce clumps as small as 0.1 Mq, but 
more typically of few x Mq (Johnson & Bromm 2006: Yoshida, Omukai 
& Hernquist 2007). Thus, any physical process promoting HD forma- 
tion might induce the formation of metal-free, small mass stars. Possible 
locations for such stars are relic H II regions (Johnson & Bromm 2006; 
Yoshida, Omukai & Hernquist 2007), shocked gas during merger of 
halos (Vasiliev Sz Shchekinov 2007), regions shocked by protostellar 
jets (Machida ct al. 2006), SN explosions (Johnson Sz Bromm 2006) 
or during structure formation (Vasiliev & Shchekinov 2005; Johnson 
&; Bromm 2006), in presence of cosmic rays (Jasche, Ciardi & EnBlin 
2007; Stacy Sz Bromm 2007). Omukai (2001) argues that, in the absence 
of molecular hydrogen (which can be easily dissociated by feedback 
effects, Sec. 4), sufficiently massive clouds can start dynamical collapse 
by Lya emission, two-photon emission and H~ free-bound emission and 
continue to collapse almost isothermally at several thousand K along 
the "atomic cooling track", until the clump mass is possibly reduced 
to about 0.03 Mq. In conclusion, depending on the cooling mecha- 
nism, a clump of Mciump ~ 10^ ~ 10^ (H2 cooling) or Mdump ~ 
0.1 — few X 10 Mq (H or HD cooling) forms out of the collapsing 
molecular cloud. 

Once a clump has reached densities of the order of n 10^ cm~^, 
the three-body formation of molecular hydrogen becomes dominant, the 
assumption of optically thin cooling begins to break down and radiative 
transfer effects become important. Tracking the subsequent evolution 
of these clumps is a very challenging problem, as it requires the simul- 
taneous solution of the hydrodynamic equations and of line radiative 
transfer. In their pioneering work, Omukai Sz Nishi (1998) studied the 
contraction of protostellar clumps into stars, for various masses and 
initial conditions, by means of 1-D hydrodynamical calculations (see 
Fig. 11 for a summary of their results). Coincidcntally, some of their 
initial conditions are very close to the final state of Abel, Bryan &; 
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Figure 11. Evolutionary sequence for a contracting protostollar clump. Gas number 
density (upper left panel), temperature (upper right), velocity (lower left) and H2 
concentration (lower right) as a function of radial distance. (0) Initial conditions. 
(1) After 5.7 X 10^ yr thrcc-body processes are active in the central region and 
temperature inversion occurs: (2) 8.5 x 10"' yr after (1). The cloud becomes optically 
thick to some lines: (3) 2.8 x 10^ yr after (2). The central region becomes fully 
molecular: (4) 12 yr after (3). The central region becomes optically thick to H2 CIA 
continuum: (5) 0.32 yr after (4). In the midst of dissociation: (6) 2.4 x 10~^ yr after 
(5). Shortly after the core formation: (7) 4.1 x 10~^ yr after (6). Final state. See 
Omukai & Nishi (1998) for details. 



Norman (2000). The authors find that the evolution of a gravitationally 
unstable clump proceeds in a highly non-homologous fashion, with the 
central parts collapsing first. This runaway phase is induced by H2 line 
radiation cooling up to densities of n « 10^^ cm~^, and by H2 collision- 
induced emission at higher densities. The resulting gas temperature is 
nearly constant at several 10^ K and the innermost region of ~ 1 Mq 
becomes fully molecular due to the three-body reaction. At densities 
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Figure 12. Radial profiles for density, temperature, molecular fraction, and infall ve- 
locity of a protostellar clump of primordial composition collapsing at rcdshift z ~ 19. 
The density is plotted as a function of distance from the center, whereas the other 
three quantities are plotted as a function of enclosed gas mass. The density profile 
is close to the power law oc R~^'^. The characteristic features in the temperature 
profile (marked AG) are described in Yoshida et al. (2006). 



n ~ 10^^ cm~^ the clump becomes optically thick to collision-induced 
absorption and H2 dissociation works as an effective cooling agent. 
Finally, lacking further cooling mechanisms, at n ~ 10^^ cm~^, a small 
hydrostatic core of 10 ^ Mq is formed. The core mass is highly 

independent of initial conditions. 

The same results have been confirmed by other authors (e.g. Ripa- 
monti et al. 2002; Machida et al. 2006). Among them, Yoshida et al. 
(2006) and Yoshida, Omukai Sz. Hernquist (2007) pushed 3-D Simula^ 
tions up to number densities of ~ 10^^ cm~^, by improving in both 
physics implementation (a treatment of the regime in which the gas is 
optically thick to H2 lines which allows to reach such densities) and 
resolution (which allows to reach a particle mass of 60 times the Earth 
mass and scale of about 1 AU) . The run-away collapse of the internal 
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core is followed down to the scales of the central 1 M© fully molecular 
core and the 0.01 which is completely opaque to H2 lines. No 
indication of fragmentation is found. Their results are quantitatively 
summarized in Fig. 12 in terms of the gas density, H2 fraction, tem- 
perature and velocity profiles, which can be readily compared with 
earlier semi-analytical results obtained by Omukai & Nishi (1998) and 
discussed in Fig. 11. In spite of some minor discrepancies, the degree 
of agreement between the two widely different methods is outstanding. 

As the core is surrounded by a large amount of reservoir gas (cor- 
responding to the mass of the initial clump), in the absence of any 
effect quenching accretion, the protostar can grow by several orders 
of magnitude in mass by accreting the envelope matter. It has yet to 
be firmly established whether standard mechanisms proposed to halt 
the infall continue to work under primordial conditions: the radiation 
force could be opacity limited; bipolar flows need some magnetohydro- 
dynamic acceleration process, and therefore seem to be excluded by the 
weak primordial magnetic field (e.g. Gnedin, Ferrara & Zweibcl 2000; 
Langer, Puget &: Aghanim 2003; Banerjee &: Jedamzik 2005; Langer, 
Aghanim & Puget 2005). 

The mass accretion rate onto the protostar can be written as Mace ~ 
c^g/G (Stabler, Shu & Taam 1980), where is the isothermal sound 
speed of the protostellar clump. Thus, the mass accretion rate is higher 
for protostars formed by atomic cooling because of the higher tem- 
perature of the protostellar clump, although the mass available for 
accretion is smaller. The evolution of zero-metal protostars in the main 
accretion phase has been recently studied by Omukai & Palla (2001) 
under a constant mass accretion rate of 4.4 x 10~^ Mq yr~^. The high 
accretion rate used, combined with a low opacity of the infalling gas 
due to the lack of dust grains, implies a milder effect of radiation 
pressure by protostellar photons and a higher inflow momentum. In 
fact, fast accreting protostars enter a phase of rapid expansion at a 
mass of K, 300 when the luminosity becomes close to the Eddington 
limit. This event may determine the onset of a powerful stellar wind 
driven by radiation pressure that can effectively quench further accre- 
tion. The effect of a timc-dcpcndcnt accretion rate, initially as high as 

10~^ M0 yr~^ and rapidly decreasing once the stellar mass exceeds 
f5i 90 M0, has been studied by Omukai & Palla (2003). They find 
that if Mace < Merit ~ 4 X 10~^ Mq yr~^, stars with mass > 100 M0 
can form, provided there is sufficient matter in the parent clumps. For 
Mace > Merit instead, the maximum mass limit decreases with Maec- In 
the context of cosmological simulations, a large scatter in the mass ac- 
cretion rates (10"'^ — 10~^ Mq/yt) is found, with a general trend of lower 
accretion rates at higher redshift (O'Shea & Norman 2007). Although 
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Figure 13. Summary of the relevant physical processes and possible evolution paths 
for Pop III star formation. 



this points toward a minimum mass of Pop III stars ^ 1 Mq, the 
upper limit remains uncertain. Similar numbers (10~^ — 10""*^ M0/yr) 
are obtained by Yoshida et al. (2006), who also estimate a final stellar 
mass of 60-100 M©. 

An ingredient which has so far received little attention but which 
could be quite important to determine the final mass of the first stars, 
is rotation of the infalling envelope. If protostellar clumps are rotating 
above a critical value of the rotational energy, which is found to be 
about 10~^ times the gravitational energy (Machida et al. 2007), frag- 
mentation of the clump might occur. Under the conditions prevailing 
at the time of first stars formation, such value is not implausible and as 
a results a non-neglible fraction of this ancient stellar population might 
be constituted by binary or multiple systems. In addition, as pointed 
out by Tan &: McKee (2004), rotation might also lead to dramatic 
differences in the way the H ii region breaks out of the infalling enve- 
lope and in radiation pressure effects. Radiation pressure, in particular, 
appears to be the dominant mechanism for suppressing infall, becoming 
dynamically important around 20 M©. 
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Table I. Classification of stars. The columns indicate: the 
name of the stellar type, their typical metallicity Z and mass 
M*, their formation mechanism or site. 



Name 


Z 


M. [M0] 


Formation 


Pop Ilia 





Z 100 


minihalos, H2 cooling 


Pop Illb 





^ 10 


HD, H cooling 


Pop II.5 


i$ Zcrit 


>io 


pre-enriched SF sites 


Pop II 




local IMF 


metal enriched SF sites 



As previously mentioned, primordial magnetic fields seem to be 
extremely weak. But during the collapse magnetic pressure increases 
more rapidly than gas pressure. For this reason, although the initial 
magnetic field is small, it can be amplified by 10 orders of magnitude 
and thus affect the formation of a star. Machida et al. (2006) use 3-D 
magneto-hydrodynamical (MHD) simulations to follow the collapse of 
a slowly rotating spherical cloud of 5 x 10^ with central number 
density ric, permeated by a uniform magnetic field. They find that if 
the initial field is > 10~^(nc/10^cm~^)^/^ G, a protostellar jet forms 
blowing away 3-10% of the accreting matter and reducing the final mass 
of the star. Studies in this area are still in their infancy and definite 
conclusions have to await for more detailed calculations. To summarize, 
if Mace < Merit protostellar cores can accrete all the matter of the 
parent clumps and possibly become very massive stars of M^, ~ Mdumpi 
otherwise accretion is less efficient and M*. <^ -^dump- Fig- 13 
a summary of the relevant physical processes and possible paths for 
Pop III star formation is sketched. A possible classification scheme of 
Pop III and Pop II stars has been recently agreed upon at the First 
Stars III Conference held in Santa Fe (O'Shea et al. 2008). A slightly 
modified version of such classification is reported in Tab. I. Here III, 
II. 5 and II indicate a metallicity sequence, while (differently from the 
original classification), a and b indicate a temporal sequence, i.e. first 
and second generation of Pop III stars. 

3.1.1. Initial Mass Function 

As we have seen above, recent studies seem to indicate that the first 
stars were massive. Nevertheless, as the primordial fragmentation pro- 
cess is very poorly understood, the distribution of masses with which 
stars are formed, the so-called IMF, is still very uncertain. The deter- 
mination of the IMF is also of basic importance because, although the 
luminosities of galaxies depend primarily on stars of masses fsi 1 Mq, 
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metal enrichment and feedback effects on galactic scales depend on the 
number of stars with masses above ^ 10 Mq. 

For many years, beginning with Schwarzschild &: Spitzcr (1953), 
there have been speculations in the literature that the IMF was domi- 
nated by massive stars at early times. More recently, Larson (1998) has 
proposed a modified Salpeter IMF for the first stars. If we define A'^ as 
the number of stars formed per logarithmic mass interval, the standard 
parameterization of the Larson IMF reads: 

dN/dlogM^ DC (1 + M*/Mc)-^■^^ (15) 

where Mc is a characteristic mass scale, related to the Jeans mass, and 
more in general to the scale at which there is a transition from a chaotic 
regime dominated by non-thermal pressure on larger scales to a regular 
regime dominated by thermal pressure on smaller scales. As the Jeans 
mass depends on temperature and pressure of a star-forming cloud as 
variable mass scale Mc naturally arises with a higher value 
at earlier times. The essential efi^ect of this type of variability would 
be to alter the relative number of low-mass to large-mass stars formed, 
with an IMF biased toward massive stars at high redshift. 

If indeed the mass of the clumps originating the protostars depends 
on the density of the collapsing cloud, as discussed earlier, this could 
produce a bimodal IMF with peaks of 10^ and 1 Mq (Nakamura & 
Umemura 2001). A similar IMF would arise in a scenario in which the 
first very massive stars create intense radiation fields which deplete the 
H2 content inside their parent clouds, decreasing, as mentioned earlier, 
the subsequent fragmentation mass scale. The resulting IMF would be 
bimodal, with a high-mass peak at ~ 40 Mq and a low-mass peak at 
a! 0.3 M0 (Omukai k Yoshii 2003). The high-mass portion of the IMF 
is a very steep function of mass, with a power-law index of -4; thus, 
although some of the first stars are very massive, their typical mass 
scale is smaller. A Larson IMF with an upper cut off of ^ 100 Mq 
seems to be favored also by theoretical models of the Pop III to Pop II 
transition (see below) which include constraints from the NICMOS UDFs 
field (Bouwens et al. 2005) and the WMAP satellite (Spergel et al. 2003). 

On the contrary, it has been known since the work of Salpeter 
(1955) that the present-day IMF of stars in the solar neighborhood 
can be approximated by a declining power-law for masses above 1 Mq, 
flattening below 0.5 Mq, and possibly even declining below 0.25 Mq 
(Scalo 1986, 1998; Kroupa 2001). While the behavior of the IMF at the 
lowest masses remains uncertain because of the poorly known mass- 
luminosity relation for the faintest stars, above ~ 1 Mq there is a 
consensus on the universality of the power- law part of the IMF, with 
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Figure 14- Temperature evolution of prestellar clouds with different metallicities. 
Those with metallicities [Z/H]=-oo (Z=0), -5, -3, and -1 (-6, -4, -2, and 0) are shown 
by solid (dashed) lines. Only the present-day CMB is considered as an external 
radiation field. The lines for constant Jeans mass are indicated by thin dotted 
lines. The positions where the central part of the clouds becomes optically thick 
to continuum self-absorption is indicated by the thin solid line. The intersection of 
the thin solid line with each evolutionary trajectory corresponds to the epoch when 
the cloud becomes optically thick to the continuum. To the right of this line, the 
clouds are optically thick and there is little radiative cooling. See Omukai et al. 
(2005) for details. 

dN/dlogM^ oc (von Hippel et al. 1996; Hunter et al. 1997; 

Massey 1998). 

Following the early study by Yoshii & Sabano (1980), it has been 
proposed (Bromm et al. 2001; Schneider et al. 2002; Bromm & Loeb 
3003; Santoro & Shull 2006; Smith & Sigurdsson 2007) that the mech- 
anism inducing a transition from a top-heavy to a more conventional 
IMF is metal enrichment. In fact, it is only when metals change the 
composition of the gas that further fragmentation occurs producing 
stars with significantly lower masses. Bromm et al. (2001) have nu- 
merically studied the effect of metals on the evolution of the gas in 
a collapsing dark matter minihalo, using an SPH simulation, for two 
different values of metallicity. They have shown that fragmentation of 
the gas cloud proceeds very differently in the two cases and concluded 
that exists a critical metallicity, Zcr, which marks a transition from 
a high-mass to a low-mass fragmentation mode. The value of Zcr is 
estimated to be 10~^ — 10~^ Zq. These values should be taken as a 
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first order approximation since the model does not include some key 
cooling agents, such as H2, other molecules and dust grains. Santoro &; 
Shull (2006), though, found similar values including cooling from CII, 
01, Sill and Fell lines and from H2. The same values are found by Smith 
k. Sigurdsson (2007) using the code ENZO to follow the formation and 
fragmentation of clouds starting from cosmological initial conditions. 
Schneider et al. (2002) though have investigated the problem including 
a much larger number (> 500) of chemical reactions. Using the 1-D 
code described in Omukai (2000), they have studied the fragmentation 
properties of a gravitationally unstable gas cloud for different initial 
metallicities in the range — 1 Zq. Including the effect of H2 and 
molecular cooling as well as the presence of dust grains in the thermal 
evolution, they concluded that the transition between the two fragmen- 
tation modes takes place for metallicities in the range 10^^ — 10^^ Zq. 
This critical threshold is almost two orders of magnitude lower than 
that found by Bromm and collaborators because the fragmentation 
properties of a gas cloud with Z = Z^x depend on the fraction of metals 
depleted onto dust grains rather than on the metals present in the gas 
phase (Schneider et al. 2003). A more recent study by Omukai et al. 
(2005), which includes also the chemistry of D previously neglected, 
confirms the above results (sec Fig. 14). The equation of state provided 
by the authors at different metallicities has been used in SPH simula- 
tions of primordial star formation by Tsuribe h Omukai (2006) and 
Clark, Glover & Klessen (2008). Both simulations confirm the results 
of the 1-D calculations (see also Tsuribe & Omukai 2008) but show 
also that, if substantial rotation or elongation of the collapsing core 
is present, fragmentation can happen at metallicities Z > 10~^ Zq 
(Tsuribe & Omukai 2006) or even lower. 

Finally, Omukai & Palla (2003) have considered the effect of heavy 
elements with abundances in the range Z = 5 x 10~^ — 5 x 10~^ Zq 
on the accretion process. The main evolutionary features of protostars 
are similar to those of metal-free objects, except that the value of 
Merit increases for metal-enriched protostars. Since the accretion rate 
is lower in a slightly polluted environment, the condition Mace < Merit 
is expected to be more easily met and the formation of massive stars 
is favored, provided the clump mass is large enough. 

3.2. Emission Spectrum 

The metal-free composition of the first stars restricts the stellar energy 
source to proton-proton burning rather than to the more efficient CNO 
cycle. Consequently, they are hotter and have harder spectra then their 
present-day counterparts of finite metallicity. These unique physical 
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Figure 15. Synthetic spectra of Pop II and Pop III clusters of lO'' M,;) with a Salpeter 
IMF. The numbers in the lower left corner, near each continuum mark, represent 
the rate (in units of photons s~^) of ionizing photons production (log) for that 
continuum (Tumlinson & Shull 2000). 



characteristics enhance the ionizing photon production of Pop III stars, 
particularly in the He ii continuum, in which they produce up to 10^ 
times more photons than Pop II (see Fig. 15 for a comparison of Pop II 
and Pop III emission spectra). If, in addition, the first stars were also 
massive, they would be even hotter and have harder spectra. Bromm, 
Kudritzki Sz Loeb (2001) find that metal-free stars with mass above 
300 Mq resemble a blackbody with an effective temperature of f» 10^ K, 
with a production rate of ionizing radiation per stellar mass larger by 
~ 1 order of magnitude for H and He i and by ~ 2 orders of magnitude 
for He II than the emission from Pop II stars. In the less extreme case of 
metal-free stars with masses < 100 Mq, Tumlinson, Shull & Venkatesan 
(2003) predict that the H-ionizing photon production takes twice as 
long as that of Pop II to decline to 1/10 of its peak value. Nevertheless, 
due to the redward stellar evolution and short lifetimes of the most 
massive stars, the hardness of the ionizing spectrum decreases rapidly. 
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leading to the disappearance of the characteristic He ll recombination 
hnes after ~ 3 Myr in instantaneous bursts (Schacrcr 2002, 2003). 

Models of emission spectra from primordial stars have been cal- 
culated by different authors for stars with masses up to « 50 Mq 
(Brocato et al. 2000 and Cassisi, Castellani Sz Tornambe 1996 in the 
implementation of Ciardi et al. 2001), ~ 100 Mq (Tumlinson & Shull 
2000; Tumlinson, Shull & Venkatesan 2003) and ^ 1000 Mq (Bromm, 
Kudritzki & Loeb 2001; Marigo et al. 2003). In particular, Marigo et al. 
(2003) have carried out extensive evolutionary calculations over a large 
range of stellar masses (0.7 Mq < < 1000 Mq), covering the H- 
and He-burning phases and allowing for a moderate overshooting from 
convective cores, but neglecting rotation. Additionally, for very massive 
stars {M^ > 120 Mq), they apply recent mass-loss rate prescriptions for 
the radiation-driven winds at very low metallicities and the amplifica- 
tion of the loss rate caused by stellar rotation. Nevertheless, the above 
studies are based on some strong simplifying assumptions; e.g. all the 
stars are assumed to be on the Zero Age Main Sequence (ZAMS) (i.e. 
stellar evolution is neglected) and nebular continuum emission is not 
included. In particular, this last process cannot be neglected for metal- 
poor stars with strong ionizing fluxes, as it increases significantly the 
total continuum flux at wavelengths redward of Lya and leads in turn to 
reduced emission line equivalent widths (Schaerer 2002, 2003). Nebular 
emission has been included in a more complete and extended study by 
Schaerer (2002, 2003), who presents realistic models for massive Pop III 
stars and stellar populations based on non-LTE model atmospheres, 
recent stellar evolution tracks and up-to-date evolutionary synthesis 
models, including also different IMFs. For more recent tabulated emis- 
sion rates refer also to Heger &: Woosley (2008). In Table II we have 
summarized the emission properties of Pop HI stars. The numbers have 
been derived integrating the ionizing photon rate in the absence of 
stellar winds given in Schaerer (2002) over different IMFs, i.e. Salpeter, 
Larson and Gaussian. 

3.3. Final Fate 

As the pp-chain is never sufficiently efficient to power massive stars, 
stars of initial zero metallicity contract until central temperatures > 
10^ K are reached and CNO seed isotopes are produced by the triple- 
a process. As a consequence of their peculiar behavior in hydrogen 
burning, the post-main sequence entropy structure of stars of zero 
initial metallicity is different from that of stars having a metallicity 
above 10~^Zq. A major uncertainty in the evolution of these stars is 
mass loss. Although radiative mass loss is probably negligible for stars 
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Table II. ZAMS Properties of Pop III stars. The columns show: the 
IMF type; the IMF mass interval [Mg]: the characteristic (central) 
mass scale for the Larson (Gaussian) IMF [M©]; rms deviation for 
the Gaussian IMF [M©]; logarithm of H i -, Ho i He ii -ionizing 
and H2-dissociating photons per baryon in Pop HI stars. 
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of such low metallicity, they stih might lose an appreciable fraction 
of their mass because of nuclear-driven pulsations. Smith & Owocki 
(2006) started from the observational evidence that nebulae around 
luminous blue variables (LBV) and LBV candidates (e.g. rj Car) show 
high ejecta masses, to discuss the possibility that the mass loss during 
the evolution of very massive stars may be dominated by optically 
thick, continuum-driven outbursts or explosions, rather than by steady 
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line-driven winds. Unlike the latter, the former may be independent of 
metallicity and play an important role in the evolution of Pop III stars. 
In any case, recent theoretical analysis on the evolution of metal-free 
stars predict that their fate can be classified as follows (see also the 
comprehensive paper by Heger k. Woosley 2008): 

— Stars with masses 10 Mq < <40 Mq proceed through the 
entire series of nuclear burnings accompanied by strong neutrino 
cooling: hydrogen to helium, helium to carbon and oxygen, then 
carbon, neon, oxygen and silicon burning, until finally iron is pro- 
duced. When the star has built up a large enough iron core, exceed- 
ing its Chandrasekhar mass, it collapses, followed by a supernova 
explosion (Woosley k. Weaver 1995). In particular, stars with > 
30 would eventually collapse into a Black Hole (BH) (Woosley 
& Weaver 1995; Fryer 1999). 

— For stars of 40 M© < My, <100 M© the neutrino-driven explosion 
is probably too weak to form an outgoing shock. A BH forms 
and either swallows the whole star or, if there is adequate angular 
momentum, produces a jet which could result in a Gamma Ray 
Burst (GRB; Fryer 1999). 

— Stars with > 100 (see e.g. Portinari, Chiosi & Bressan 
1998) form large He cores that reach carbon ignition with masses in 
excess of about 45 Mq. It is known that after helium burning, cores 
of this mass will encounter the electron-positron pair instability, 
collapse and ignite oxygen and silicon burning explosively. If ex- 
plosive oxygen burning provides enough energy, it can reverse the 
collapse in a giant nuclear-powered explosion (the so-called Pair 
Instability Supernova, PISN) by which the star would be partly 
or completely (if > 140 Mq) disrupted (Fryer, Woosley k. 
Heger 2001). The onset of a PISN though depends whether the 
star can retain enough of its initial mass to have a helium core 
bigger than about 64 Mq at the end of the core He-burning phase, 
or rather lose the mass during the previous stages of its evolution 
(Meynet, Ekstrom &: Maeder 2006). For even more massive stars 
(M* > 260 Mq) a new phenomenon occurs as a sufficiently large 
fraction of the center of the star becomes so hot that the photo- 
disintegration instability is encountered before explosive burning 
reverses the implosion. This uses up all the energy released by 
previous burning stages and accelerates the collapse leading to the 
prompt formation of a massive BH, and, again, either a complete 
collapse or a jet-powered explosion (Bond, Arnett and Carr 1984; 
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Fryer, Woosley & Heger 2001; Ohkubo et al. 2006; Suwa et al. 
2007a). 

— At even higher masses (M* > 10^ M©) the evolution depends on 
the metaUicity (Fuller, Woosley & Weaver 1986): if Z < 0.005 the 
star collapses to a BH as a result of post-Newtonian instabilities 
without ignition of the hydrogen burning; for higher metallicities 
it explodes, as it could generate nuclear energy more rapidly from 
/^-limited cycle. 

A summary of the final fate of Pop III stars is sketched in Fig. 16. 

3.4. Black Holes and Gamma Ray Bursts from the First 

Stars 

As discussed in the previous Section, metal-free (and also moderately 
metal-poor) stars more massive than approximately 260 Mq collapse 
completely to BHs. Similar arguments apply to stars in a lower mass 
window (30 Mq — 140 Mq), which are also expected to end their evolu- 
tion as BHs. If this is the case, a numerous population of Intermediate 
Mass Black Holes (IMBHs) - with masses M, « 10^ - 10^ Mq, between 
those of stellar and SuperMassive Black Holes (SMBHs) - may be 
the end-product of those episodes of early star formation and offer 
a possible observable of primordial stars. As these pre-galactic BHs 
become incorporated through a series of mergers into larger and larger 
halos, they sink to the center because of dynamical friction, accrete 
a fraction of the gas in the merger remnant to become supermassive, 
form a binary system, and eventually coalesce. 

Volonteri, Haardt &i Madau (2003) thoroughly studied this scenario 
including a number of physical ingredients such as gas accretion, hard- 
ening of the IMBH binary and triple interactions. Their results show 
that this hierarchical growth can reproduce the observed luminosity 
function of optically selected quasars in the redshift range 1 < z < 5. 
A prediction of the model is that a population of "wandering" black 
holes in galactic halos and in the IGM should exist (and possibly ob- 
served), contributing around 10% of the present day total black hole 
mass density, 4 x 10^ Mq Mpc~^ (a similar study has been performed 
by Islam, Taylor & Silk 2003). IMBHs that have not yet ended up in 
SMBHs could also be either (i) en route toward galactic nuclei, thereby 
accounting for the X-ray-bright off-center sources detected locally by 
ROSAT-*^^, or (ii) constituting the dark matter candidates composing the 
entire baryonic halos of galaxies (Schneider et al. 2002). A possible 

^® http://wave.xray.mpe.mpg.de/rosat 
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Figure 16. Initial-final mass function of Pop III stars (Hcgor & Woosley 2002). The 
horizontal axis gives the initial stellar mass. The j/-axis gives the final mass of the 
collapsed remnant (thick black curve) and the mass of the star at the beginning of 
the event producing that remnant (thick gray curve) . Since no mass loss is expected 
from Pop III stars before the final stage, the gray curve is approximately the same 
as the lino of no mass loss (dotted). Exceptions are « 100 — 140 M0 where the 
pulsatioiial pair-instability ejects the outer layers of the star before it collapses, 
and above 500 M© where pulsational instabilities in red supergiants may lead to 
significant mass loss (dashed) . The hydrogen- rich envelope and parts of the helium 
core (dash-double-dotted curve) are ejected in a supernova explosion. 



method to constrain the existence and nature of BHs at z > 6 would 
be with the next generation of X-ray and NIR space telescopes, that 
could detect BHs down to masses of 10^ - 10^ Mq (Salvaterra, Haardt 
& Volonteri 2007). 

Once a proto-BH has formed into the stellar core, accretion con- 
tinues through a disk. It is widely accepted, though not confirmed, 
that magnetic fields drive an energetic jet (Fig. 17) which produces 
a burst of TeV neutrinos (Linke et al. 2001; Meszaros, & Waxman 
2001; Schneider, Guetta k, Ferrara 2002) by photon-meson interaction. 
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and eventually breaks out of the stellar envelope appearing as a GRB. 
Based on recent numerical simulations and neutrino emission models, 
the expected neutrino diffuse flux from these Pop III GRBs could be 
within the capabilities of present and planned detectors as AMANDA^^ 
(Antartic Muon And Neutrino Detector Array) and IceCube^'^. High- 
energy neutrinos from Pop III GRBs could dominate the overall flux 
in two energy bands, 10^ — 10^ GeV and 10^ — 10® GeV, of neutrino 
telescopes. The enhanced sensitivities of forthcoming detectors in the 
high-energy band (AMANDA-II, IceCube) will provide a fundamental 
insight into the characteristic explosion energies of Pop III GRBs, and 
will constitute a unique probe of the IMF of the first stars. Based 
on such results, Pop III GRBs could be associated (e.g. Schneider, 
Guetta & Ferrara 2002; Arefiev, Priedhorsky & Borozdin 2003) with a 
new class of events detected by BeppoSax^^, the Fast X-ray Transients 
(FXTs), bright X-ray sources with peak energies in the 2 — 10 keV band 
and durations between 10 s and 200 s (e.g. Heise et al. 2001). locco 
et al. (2007) though, using updated models of Pop III star formation 
history and neutrino emission yields show that the planned generation 
of neutrino telescope will not be able to detect the contribution from 
Pop III stars and that, in any case, this would be highly contaminated 
by the contribution from Pop II stars. 



3.5. Key Observations 



Despite the efforts made in the past decades, the search for zero- 
metallicity stars has been proven unfruitful as no Pop III star has been 
detected yet. The first very low metallicity star, with a logarithmic 
abundance relative to solar of [Fe/H]=— 4.5^^, was identified more than 
15 years ago (Bessel &: Norris 1984). Since then, more stars with —4.5 < 
[Fe/H] < —2 have been detected (McWilliam et al. 1995; Beers et al. 
1998; Norris, Ryan & Beers 2001; Carretta et al. 2002). The record was 
established by Christlieb et al. (2002), who measured an abundance of 
[Fe/H] =-5.3 in HE0107-5240, a star with a mass oi^ 0.8 Mq. Does the 
existence of this star suggest that Pop III also contained low-mass and 
long-lived objects? To answer this question, we should be able to infer 
the metallicity of the gas cloud out of which HE0107-5240 formed from 
the abundance of the elements observed on its surface. In particular, 
the observed abundance pattern in elements heavier than Mg (which 



http : // amanda . berkeley . edu/www/ amemda . html 

http : / / icecube . wise . edu/ 

http : www . asdc . asi . it/bepposax/ 

In the usual notation, [X/H] = log (X/H) - log (X/H) 
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Figure 17. Collapse of a 300 Mq Pop III star (Fryer, Woosley & Heger 2001). The 
proto-BH 0.5 s before BH formation is shown. Color denotes temperature in 10^ K, 
and the vectors represent the direction and magnitude of the particle velocity. At 
this time, the proto-BH has a mass of roughly 78 M© and size of 1100 km. The inner 
core first forms a black hole, but as soon as it collapses, the entire proto-BH (which 
is 90 M0 at collapse) quickly accretes onto the BH. 

cannot be synthesized in the interior of a 0.8 Mq star) is a distinct 
signature of the previous generation of stars. 

HE0107-5240 is very iron-deficient but overabundant in C, N and 
O (Bessel et al. 2004; Christheb et al. 2004). The origin of these light 
elements can be twofold: (i) either they were already present in the 
parent gas cloud out of which the star was born, or (ii) their origin 
is due to a post-formation mechanism, such as mass transfer from a 
companion star, self-enrichment from the star itself or accretion due to 
repeated passage through the Galactic disk (for a thorough discussion 
of these issues we refer to Christlieb et al. 2004). In the first scenario, 
the gas cloud out of which the star formed would have had to be pre- 
enriched to a metallicity as high as Z = IO^^Zq. This metallicity is 
significantly higher than Zcr and low-mass stars can form. Moreover, 
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Umeda & Nomoto (2005) have shown that if the polluter is a zero- 
metallicity star with mass in the range ~ 25 — 130 Mq exploding as a 
sub-luminous Type-II SN, a good fit of the observed abundance pattern 
can be obtained (although some problems remain with the updated 
oxygen abundance; Bessel et al. 2004; Christlieb et al. 2004). In the 
second scenario, only the elements heavier than Mg, which cannot 
be synthesized in the interior of such a low-mass star, were already 
present before its formation. The composition of these elements can 
be equally well reproduced by the predicted elemental yields of either 
an intermediate mass star that exploded as SN (Christlieb et al. 2002; 
Limongi, Chieffi &: Bonifacio 2003; Weiss et al. 2004; Iwamoto et al. 
2005) or of a 200 - 220 Mq star that exploded as a PISN (Schneider 
et al. 2003). Thus, HE0107-5240 may have formed from a gas cloud 
which had been pre-enrichcd by such supernovae up to an estimated 
metallicity of 10^^'^ < Z/Zq < 10^^'^. These values fall within the 
proposed range for Z^r- 

A second very iron deficient star (HE1327-2326, with [Fe/H]= —5.4) 
with an overabundance of C and N has been detected by Frebel et al. 
(2005). Differently from HE0107-5240, which is a giant and has evolved 
off the main-sequence, HE1327-2326 is relatively unevolved. 

The most massive star is the Pistol star, with an estimated initial 
mass of 200 — 250 Mq (Figer et al. 1998). Although there are no in- 
dications that the star is a binary, it might have recently experienced 
a rejuvenation through a merger with another star. Figer (2005) has 
investigated the presence of very massive stars in the Arches cluster, 
which has an estimated age of 2-2.5 Myr. Thus, stars with masses above 

150 Mq should still be visible if they were formed. The typical mass 
function would predict 18 stars with masses greater than 130 Mq, but 
none were found, suggesting a firm upper limit for the star masses of 
150 M0. Also Oey & Clarke (2005), based on a statistical analysis of 
OB associations, find that, if the IMF has a slope close to a Salpeter, 
it would be extremely unlikely to have stars with masses greater than 
120-200 M0. 

As a novel strategy to search for Pop III stars, Scannapieco, Schnei- 
der &: Ferrara (2003) have suggested that at least some intermediate 
redshift Lya-emitting galaxies might have their flux powered by Pop III 
stellar clusters. In Fig. 18 the isocontours in the Lya luminosity-redshift 
plane indicate the probability to find Pop III star-hosting galaxies in 
a given sample of Lya emitters for various feedback efficiencies. Eg. 
Such objects populate a well defined region of the Lya-redshift plane, 
whose extent is governed by the feedback strength. Furthermore, their 
fraction increases with redshift, independently of the assumed feedback 
strength, suggesting that Lya emission from already observed high-z 
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Figure 18. Fraction of Pop III star-hosting galaxies as a function of Lya luminosity 
and redshift. Isocontours of fractions > 10~^, 10"^'^, 10~^ and lO"'^ are shown for a 
burst-mode star formation. Each panel is labeled by the assumed feedback parameter 
Eg value. For reference, the dashed line gives the luminosity corresponding to an 
observed flux of 1.5 x 10~^^ ergs cm~^ s~^, and the various points correspond to 
observed galaxies (see details in Scannapieco, Schneider & Ferrara 2003). 



sources can be due to Pop III stars. Hence, collecting large data samples 
to increase the statistical leverage may be crucial for detecting the 
elusive first stars. 

Confirmation of the presence of Pop III stars in the Lya emitters can 
come from follow-up spectroscopy. Tumlinson and collaborators (Tum- 
linson, Giroux &: Shull 2001; Tumlinson, Shull & Venkatesan 2003) 
have suggested that the harder ionizing spectrum expected from metal- 
free stars would result in stronger He ii recombination lines. Their 
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Figure 19. Flux of He ii A1640 line as a function of the source redshift. The different 
lines refer to different models for the time evolution of the stellar ionizing radiation. 

The horizontal lines show the limits of the current Lya emission-line surveys. At 
the top, the ranges of redshift probed by the He ii A1640 (above) and A4686 (below) 
lines for optical and infrared searches are marked. See Tumlinson, Giroux & ShuU 
(2001) for details. 

calculations indicate that line fluxes of the He II A1640 and A4686 
are sufficiently large to be detected by narrow-band and spectroscopic 
searches for high-redshift emission-line sources (see Fig. 19). Together 
with measurements of the Ha emission line, the He ii recombination 
lines can be used also to infer the ratio of the He ii to Hi ionizing 
photons (Oh, Haiman & Rees 2001). 

Nagao et al. (2008) performed a photometric search of Pop HI stars 
by looking at galaxies with strong emission both in Lya and He ii A1640 
in the redshift intervals 3.93 < z < 4.01 and 4.57 < z < 4.65. As none 
of the observed dual emitters could be identified with a galaxy hosting 
a Pop HI star, an upper limit of 5 x 10~^ yr~^ Mpc~^ to the Pop HI 
SFR has been set. 

In addition to He lines, also H2 lines could be used to observe 
the first stars. For example, H2 lines (principally J = 5 — 3(?; = 0), 
J = 2 — 0(f = 0) and Aw = 1 lines) emitted by collapsing primordial 
clouds, although not observable from an individual object, could pro- 
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duce cumulative fluxes in the sub-mJy level, marginally observable with 
ASTRO-F^''^ and ALMA^^ (Atacama Large Millimeter Array). If HD is an 
important coolant as discussed in the previous Section, then additional 
line emission might be observed. In this case, the lines emitted during 
the collapse would show a double-peak feature, with the higher peak 
associated with the H2 J = 2 — 0(v = 0) transition and the lower peak 
with the HD J = 4 — 3{v = 0) transition (Kamaya & Silk 2002; Kamaya 
Sz. Silk 2003). A similar strategy has been proposed by Ripamonti et 
al. (2002), who suggest that the sum of the H2 infrared line and con- 
tinuum emission radiation produced during the collapse and accretion 
on a small hydrostatic protostellar core could be detected by JWST. In 
contrast with the above studies, Mizusawa, Nishi Sz Omukai (2004), 
who include in their calculation the emission of H2 lines both in the 
collapse and in the accretion phase, find that observation of such lines 
with forthcoming facilities will be highly improbable. 

While direct observations of the early stars are very difficult (if not 
impossible) to obtain, one can search for evidence of their past existence 
with other observable. Such evidence can in particular be seen in binary 
systems, where one of the components is a sub-solar-mass star and 
the other was a larger mass star. The presently observed spectrum of 
the lower-mass companion would exhibit the chemical signature of the 
larger-mass companion. The high frequency of carbon-enhanced stars 
found among very metal-poor stars in the halo of our Galaxy provides 
a constrain on the past IMF and indicates that it was different from the 
present and shifted toward higher masses (e.g. Lucatello et al. 2005). 
These findings are consistent with an early prediction by Hernandez 
& Ferrara (2001) of an upward shift of the characteristic stellar mass 
with rcdshift based on an analysis of the observational data of Galactic 
metal poor stars. Such mass increases above the Jeans mass fixed by 
the increase of the CMB temperature with redshift. 

Ballero, Matteucci & Chiappini (2006) though reach a different con- 
clusion. They have studied the effect of Pop III star metal pollution 
on the chemical evolution of the Milky Way, focusing their analysis 
on a-elements, carbon, nitrogen and iron, finding that the effect of 
Pop III stars is negligible and it is not possible to prove or disprove their 
existence based on chemical abundance in the Milky Way. Moreover, 
their analysis favors a constant or slightly varying IMF to a strongly 
variable one. Along the same lines, Damped Lyman Alpha (DLA) 
systems have been searched for the nucleosynthetic ashes of Pop III 
stars. Erni et al. (2006) compare the chemical enrichment pattern of 



http:/ /www.ir.isas.ac.jp/ASTRO-F/index-e.html 
http://www.alma.nrcio.edu 
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the most metal-poor DLAs with predictions from different explosive 
nucleosynthesis models and find that the observed abundance pattern 
favors massive (10-50 Mq) mctal-free supernovac/hypernovae. Pair- 
instability SNe seem to be excluded; however, given the metallicity 
level of DLAs (p^ 10~^ Zq), this is not surprising as it is likely that the 
gas in these systems has already been enriched by a large number of 
stellar populations. 

Signatures of Pop III stars might be seen in puzzling observations of 
the abundance of lithium isotopes (e.g. Asplund et al. 2006). While ^Li 
is produced both in BBN and through cosmic ray (CR) spallation, ®Li 
is produced only via the second process. On this basis a correlation be- 
tween ^Li abundance and metallicity is expected. Two problems arise: 
(i) ^Li/H abundance observed in metal-poor halo stars is '-^ 10^^'^, 
i.e. much lower than expected from BBN alone ('^ 4.5 x 10^^*^); (ii) 
^Li abundance does not correlate with metallicity and rather shows a 
plateau that extends to low metallicity values. Different scenarios have 
been proposed to explain the plateau (CR in shocks during the forma- 
tion of the Galaxy, decay of relic particles during BBN, production in 
connection with gamma rays, initial burst of cosmological CRs) without 
overproducing ^Li. A tentative solution has been found by Rollinde, 
Vangioni &: Olive (2006) who have calculated the production of ^Li 
from CR produced by Pop III and Pop II stars using a prescription for 
the cosmic SN rate. However, the models makes extreme and somewhat 
unphysical assumptions concerning the cosmic star formation history 
and it is not clear to what extent it can reproduced the metallicity 
distribution function of halo stars. Useful reviews for readers interested 
in the subject can be found in Coc et al. (2004) and Steigman (2005). 

Jimenez & Haiman (2006) explain with the presence of 10%-30% of 
metal-free stars at 2; ~ 3 — 4, four puzzling observations that cannot 
be accounted for with a more normal stellar population: (i) the signif- 
icant UV emission from LBGs at A < 912 A; (ii) the strong Lyman-a 
emission from extended blobs at z = 3.1 with little or no associated 
apparent ionizing continuum; (Hi) a population of 2; ~ 4.5 galaxies 
with a very wide range of Lyman-a equivalent width emission lines; 
(iv) strong Hell emission line in a composite of LBGs. 

GRBs have been proposed as possible tools to study the high-z 
universe and, in principle, a fraction of them could be triggered by 
Pop HI progenitors. To produce a GRB though, a star needs to lose its 
outer H envelope. As the mass loss from Pop HI stars is not expected 
to be significant, mass transfer to a companion star in binaries could 
help. If Pop III binaries are common. Pop III star formation at high-z 
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Figure 20. Extragalactic background light data compilation from various IR and 
UV experiments. Also shown are (lower) limits coming from galaxy counts (open 
symbols). Lines denote different theoretical models (for details see Aharonian et al. 
2006). 

could be probed by SWIFT^^ (Bromm & Loeb 2006). Belczynski et al. 
(2007) find that, although, about 10% of Pop III binaries lose the H 
envelope, only 1% acquire sufficient angular momentum to produce a 
GRB. Thus, it will be extremely difficult for SWIFT to detect GRBs 
from Pop III binaries. 

In the last years the alleged discovery of an excess in the near IR 
background (an updated compilation on the extragalactic background 
light measurements is shown in Fig. 20) and of its fluctuations above 
the "normal" galaxy contribution claimed by Matsumoto et al. (2005; 
but see, e.g., Aharonian et al. 2006 for a different view), has produced 
considerable excitement given the fact that the excess could be inter- 
preted as an imprint of redshifted light of the first (massive) stars, 

http://www.swift.psu.edu/xrt/ 
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(Santos, Bromm & Kamionkowski 2002; Salvaterra &; Ferrara 2003a; 
Madau &: Silk 2005; Kashlinsky 2005). In addition, the emission from 
metal-free stars is needed to explain the level of fluctuations in the 
same background, at least at the shortest wavelengths. In fact, clus- 
tering of unresolved Pop III stars can account for the entire signal at 
almost all the ~ 1 — 30 arcscc scales probed by observations in the J 
band, while their contribution fades away at shorter frequencies and 
becomes negligible in the K band, where "normal" galaxies are dom- 
inant (Magliocchetti, Salvaterra k. Ferrara 2003; Cooray et al. 2004; 
Kashlinsky ct al. 2004; Kashlinsky et al. 2005; Cooray et al. 2007). The 
new measurements of fluctuations with the IRAC instrument on board 
of Spitzer and the associated analysis confirm that the anisotropics 
on scales ^ 0.5' — 10' must be of extragalactic origin (Kashlinsky et al. 
2007a). Salvaterra ct al. (2006) find that Pop III stars can contribute 
only less than 40% to the NIRB and that the fluctuations can be repro- 
duced by the clustering of Pop II hosting galaxies at z > 5. That the 
fluctuations must arise from a population of sources at those redshifts 
is confirmed also by Kashlinsky et al. (2007b), although they find that 
such sources should be much brighter per unit mass than the present 
stellar population. 

Dwek, Arendt & Krennrich (2005) also noticed that explaining the 
excess with Pop III contribution only, would require uncomfortably 
high formation efficiencies. In fact, more than 5% of all baryons must 
be processed into Pop III stars and, to avoid excessive metal pollution, 
they should end up in IMBHs with steep UV/X-ray spectra (not to 
overproduce the present-day soft X-ray background). In addition, if 
the excess in the NIRB were of extragalactic origin (rather than due 
to zodiacal light), its 7 — 7 interactions with 7-ray photons emitted 
by blazars would produce a distinct absorption feature in the spectra 
of these sources. Dwek, Krennrich k. Arendt (2005) show that, if the 
NIRB excess were due to Pop III stars it would be difficult to rcprodTicc 
the absorption-corrected spectra of some observed distant TeV blazars. 

A substantial IR background could arise from Pop III stars not only 
due to the direct emission associated with these stars, but also due 
to indirect processes that lead to free-free and Lya emission from the 
ionized nebulae or H 11 regions surrounding the stars. In view of more 
recent observations though, a substantial contribution to the near IR 
background for Pop III stars seems to be excluded. Salvaterra &: Ferrara 
(2006) find that if this is the case, a large fraction (in some case all) 
of Spitzer counts (at 3.5 /xm and 4.5 /Ltm; Fazio at al. 2004) at the 
faintest fluxes is due to galaxies at 2; > 8, but worse, the number of 
J-dropout and Lya emitters identified in NICMOS (Bouwens et al. 2005) 



review_08.tex; 28/04/2008; 12:20; p. 54 



55 




P t i X 14 I I I I ,i t I .« I t t. ) J V. I 

1^ 16 1^ 1© ^ m 23^ m M 



Figure 21. Predicted surface density of PISNe on the sky per unit redshift as a 

function of redshift for four different observing wavelengths (labeled in /im), assum- 
ing that survey images reach 1 mag below maximum light, after Weinmann & Lilly 
(2005). 



and ZEN (Willis &; Courbin 2005) deep surveys would be overpredicted 
by a factor of thousands. 

An additional NIRB contribution could come from supcrnovae (al- 
though this is usually minor) and high redshift quasars. The latter 
could give a substantial contribution if Pop III were indeed very massive 
and ended up as IMBH (Cooray & Yoshida 2004). In this case though 
they would probably overestimate the measured soft X-ray background 
(Salvaterra, Haardt & Ferrara 2005). 

Regardless of the precise contribution from primordial stars to the 
IR background, they would leave enough photons to provide a large 
optical depth for high-energy photons from high-z GRBs. Observations 
of their spectra could thus provide information on the emission from 
the first stars (Kashlinsky 2005). 

JWST could observe the upper mass range of pair instability super- 
novae from the first stars, which would provide, in addition to an 
observational proof of the existence of very massive metal- free stars, 
stringent constraints on the understanding of primordial star formation 
and IMF (Wise k Abel 2005). Ahhough PISN are bright enough to be 
detected by JWST, Weinmann &; Lilly (2005) claim that their rate is 
so small (0.2-4 deg^^ yi"^^) that either a dedicated wide-angle search 
or a serendipitous search is needed. The numbers quoted by the latter 
authors are less optimistic than previous estimates (0-2200 deg~^ yr~^; 
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e.g. Cen 2003a; Mackey, Bromm & Hernquist 2003; Wise & Abel 2005) 

because very approximate liglit curves were adopted. Although, as 
noted by Scannapicco ct al. (2005), the light curves used by Wcinmann 
&: Lilly (2005) contain a normalization error, the numbers derived are 
in agreement with those derived adopting the correct light curves (see 
Fig. 21). In any case. Smith at al. (2007) reported the detection of the 
peculiar Type II SN 2006gy which could be the first SN powered by a 
pair-instability explosion. 

If indeed the first stars were very massive and left BH remnants, then 
gravitational waves emitted during each collapse (Schneider et al. 2000; 
Fryer, Holz & Hughes 2002; Buonanno et al. 2005; Sandick et al. 2006; 
Suwa et al. 2007b) would be seen as a stochastic background, whose pre- 
dicted spectral strain peaks in the frequency range 5x10-^-5x10-3 Hz 
(where it has a value 10^^'^ — 10^"*^^ Hz^^/^) and might be detected by 
LISA'^^ (Laser Interferometer Space Antenna). Similarly, the coalescence 
of massive black hole binaries due to galaxy mergers has also been 
considered as a source of low-frequency gravitational radiation (Wyithe 
& Loeb 2003c; Sesana et al. 2004; Belczynski, Bulik k Rudak 2004). 
Several discrete mHz massive black hole sources detectable by LISA 
are likely to originate at redshifts z > 7, while the advanced LIGO^''' 
should detect at least several events per year with high signal to noise 
ratio (> 10). Core collapse SN from primordial stars would produce 
a gravitational wave background that peaks at 30-100 Hz (compared 
to 300-500 Hz from a normal mode of star formation) and is within 
the sensitivity of planned interferometers as LIGO, BBO (Big Bang Ob- 
servatory) and DECIGO (DECi-hertz Interferometer Gravitational wave 
Observatory) (Sandick et al. 2006). Using an improved spectrum, Suwa 
et al. (2007b) find that the range of observability is 0.1-1 Hz. 



4. Feedback Effects 

Once the first sources have formed, their mass deposition, energy in- 
jection and emitted radiation can deeply affect the subsequent galaxy 
formation process and influence the evolution of the IGM via a number 
of so-called "feedback" effects. The word "feedback" is by far one of the 
most used in modern cosmology, where it is applied to a vast range of 
situations and astrophysical objects. However, for the same reason, its 
meaning in the context is often unclear or fuzzy. Generally speaking, 
the concept of feedback invokes a back reaction of a process on itself or 
on the causes that have produced it. Secondly, the character of feedback 

^® http://lisa.jpl.nasa.gov 

^"^ http://www.ligo.caltecli.edu 



review_08.tex; 28/04/2008; 12:20; p. 56 



57 

Table III. Classification of different feedback effects. 



Radiative 


NEGATIVE 

Mechanical 


Chemical 


1. Photoionization/evaporation 

2. H2 Photodissociation 

3. Photoheating filtering 


1. Blowout /blowaway 

2. Impinging shocks 

3. Preheating 


1. Fragmentation 


Radiative 


POSITIVE 

Mechanical 


Chemical 


1. In front of H 11 regions 

2. Inside relic H 11 regions 

3. X-ray background 


1. Behind shocks 

2. Shell fragmentation 





can be either negative or positive. Finally, and most importantly, the 
idea of feedback is intimately linked to the possibility that a system 
can become self-regulated. Although some types of feedback processes 
are disruptive, the most important ones in astrophysics are probably 
those that are able to drive the systems towards a steady state of some 
sort. To exemplify, think of a galaxy which is witnessing a burst of star 
formation. The occurrence of the first supernovae will evacuate/heat 
the gas thus suppressing the subsequent SF process. This feedback 
is then acting back on the energy source (star formation); it is of a 
negative type, and it could drive the SF activity in such a way that 
only a sustainable amount of stars is formed. However, feedback can 
fail to produce such regulation either in small galaxies, where the gas 
can be ejected by a handful of SNc, or in cases when the star formation 
timescale is too short compared to the feedback timescale. In the spirit 
of the present review we will only discuss feedback occurring at high 
redshifts and hence shaping the first structures. 

Although a rigorous classification of the various effects is not feasi- 
ble, they can be divided into three broad classes: radiative, mechanical 
and chemical feedback. In the first class fall all those effects associated, 
in particular, with ionization/dissociation of hydrogen atoms /molecules; 
the second class is produced by the mechanical energy injection of 
massive stars in form of winds or SN explosions; and chemical feedback 
is instead related to the postulated existence of a critical metallicity 
governing the cosmic transition from very massive stars to "normal" 
stars. 
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Attempting a classification of all proposed feedback effects is almost 
desperate, due to the large number of applications and definitions, 
often discordant, present in the literature. Nevertheless, we offer in 
Table III a working classification aimed essentially at organizing the 
material discussed in this Section. Whereas the radiative feedback is 
intimately connected with cosmic reionization, chemical feedback is 
instead strongly dependent on the history of metal enrichment of the 
universe. Although reionization and metal enrichment could then be 
analyzed as aspects of feedback, given their importance and the large 
amount of literature on these subjects, we have reserved a separate 
discussion for them in Sec. 5. 

4.1. Radiative Feedback 

Radiative feedback is related to the ionizing/dissociating radiation pro- 
duced by massive stars or quasars. This radiation can have local effects 
(i.e. on the same galaxy that produces it) or long-range effects, either 
affecting the formation and evolution of nearby objects or joining the 
radiation produced by other galaxies to form a background. In spite of 
the different scenarios implied, the physical processes are very similar. 

Photoionization/ evaporation 

The collapse and formation of primordial objects exposed to a UV 
radiation field can be inhibited or halted for two main reasons: (i) 
cooling is considerably suppressed by the decreased fraction of neu- 
tral hydrogen, and (ii) gas can be photoevaporated out of the host 
halo. In fact, the gas incorporated into small mass objects that were 
unable to cool efficiently, can be boiled out of the gravitational po- 
tential well of the host halo if it is heated by UV radiation above the 
virial temperature. Such effects are produced by the same radiation 
field and act simultaneously. For this reason, it is hard to separate 
their individual impact on the final outcome. The problem has been 
extensively studied by several authors (e.g. Thoul &: Weinberg 1996; 
Nagashima, Gouda Sz Sugiura 1999; Barkana & Loeb 1999; Susa Sz 
Kitayama 2000; Chiu, Gnedin & Ostriker 2001; Kitayama et al. 2000, 
2001; Machacck, Bryan & Abel 2001; Haiman, Abel k Madau 2001; 
Susa & Umemura 2004b; Ahn & Shapiro 2007; Whalen et al. 2007). 
In particular, Susa Sz Kitayama (2000) and Kitayama et al. (2000, 
2001), assuming spherical symmetry, solve self-consistently radiative 
transfer of photons, non-equilibrium H2 chemistry and gas hydrody- 
namics of a collapsing halo. They find that at weak UV intensities 
(J < 10~^^ erg s"-*" cm~^ sr~^ Hz~^), objects as small as Vc ~ 15 km s~^ 
are able to collapse, owing to both self-shielding of the gas and H2 cool- 
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ing. At stronger intensities though, objects as large as Vc ~ 40 km 
can be photoevaporatcd and prohibited from collapsing, in agreement 
with previous investigations based on the optically thin approximation 
(Thoul &: Weinberg 1996). Similar results have been found by Susa &: 
Umemura (2004b) by means of a 3-D hydrodynamic simulation which 
incorporates the radiative transfer of ionizing photons. Also the 3-D 
hydrodynamic simulation of Machacek, Bryan k, Abel (2001) confirms 
that the presence of UV radiation delays or suppresses the formation 
of low mass objects. In contrast with these studies, Dijkstra et al. 
(2004), applying the 1-D code by Thoul k Weinberg (1996) to z > 10 
objects, find that objects as small as Vc ~ 10 km s~^ can self-shield and 
collapse because the collisional cooling processes at high redshift are 
more efficient and the amplitude of the ionizing background is lower. 
The second condition might not always apply, as the ionizing flux at 
high redshift is dominated by the direct radiation from neighboring 
halos rather than the background (Ciardi et al. 2000). 

The photoevaporation effect might be particularly important for 
Pop III objects, as their virial temperatures are below the typical 
temperatures achieved by a primordial photoionized gas (T ^ 10^ K). 
For such an object, the ionization front gradually burns its way through 
the collapsed gas, prodTicing a wind that blows backwards into the IGM 
and that eventually evaporates all the gas content. According to recent 
studies (Haiman, Abel k Madau 2001; Shapiro, Iliev k Raga 2004; 
Whalen et al. 2007) these sub-kpc galactic units were so common to 
dominate the absorption of ionizing photons. This means that estimates 
of the number of ionizing photons per H atom required to complete 
reionization should not neglect their contribution to absorption. As 
Fig. 22 (Shapiro, Iliev k Raga 2004) shows, the number of ionizing 
photons absorbed per initial minihalo atom, ^, increases gradually with 
time; in addition, it depends on the ionizing spectrum assumed. For 
the hard QSO spectrum, the ionization front is thicker and penetrates 
deeper into the denser and colder parts of the halo, increasing the rate 
of recombinations per atom, compared to stellar type sources. However, 
this same pre-heating effect shortens the evaporation time, ultimately 
leading to a rough cancellation of the two effects and the same total 
^ as for a black-body spectrum with T = 5 x 10^ K (mimicking a low 
metallicity Pop II stellar emission). An even lower ^ is needed for a 
black-body spectrum with T = 10^ K (more typical of a Pop III stellar 
emission), because of an increased evaporation vs. penetration ability. 
Thus, overall. Pop III stellar sources appear significantly more efficient 
than Pop II or QSO sources in terms of the total number of ionizing 
photons needed to complete the photoevaporation process. It should 
be noted though that the typical timescales for photoevaporation are 
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Figure 22. Evolution of the cumulative number of ionizing photons absorbed per 
initial minihalo atom for a QSO and two black-body spectra with T — 5 x 10* K 
and T = 10^ K, as labelled (see Shapiro, Iliev & Raga 2004 for detailed definitions). 



larger than the lifetime of massive Pop III stars, and thus minihalos 
generally survive photoevaporation (Alvarez, Bromm &; Shapiro 2006). 
The degree of photoevaporation also depends sensively on the evolution 
stage of halos, with diffuse halos with central densities below few cm~^ 
being completely photoevaporated, while more evolved halos generally 
retaining their core (Whalcn ct al. 2007). 

H2 Photodissociation 

As intergalactic H2 is easily photodissociated, a soft-UV background 
in the Lyman- Werner bands could quickly build up and have a nega- 
tive feedback on the gas cooling and star formation inside small halos 
(Haiman, Rees & Loeb 1997; Ciardi, Ferrara & Abel 2000; Ciardi et 
al. 2000; Haiman, Abel & Rees 2000; Ricotti, Gnedin & Shull 2002; 
Mackey, Bromm & Hernquist 2003; Yoshida et al. 2003a; Wise & Abel 
2008; Johnson, Greif Sz Bromm 2008). In addition to an external back- 
ground, the evolution of structures can also be affected by internal 
dissociating radiation. In fact, once the first generation of stars has 
formed in an object, it can affect the subsequent star formation pro- 
cess by photodissociating molecular hydrogen in star forming clouds 
(Omukai & Nishi 1999; Nishi & Tashiro 2000; Glover & Brand 2001; 
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Oh & Haiman 2002). In particular, Nishi & Tashiro (2000) show that 
if the molecular cloud has a metallicity smaller than about 
a single O star can seriously deplete the H2 content so that subsequent 
star formation is almost quenched. Thus, it seems plausible that stars 
do not form efficiently before the metallicity becomes larger than about 
10~^ solar. On the other hand, Susa & Umemura (2006) performed 
hydrodynamic simulations coupled to radiative transfer of ionizing and 
dissociating photons finding that, if a star forming cloud exceeds the 
threshold density of ~ 10^ cm~^ and is at a distance > 30 pc from a 
120 Mq Pop III star, it can survive negative feedback as the H2 shell 
formed in front of the ionization front is able to shield the cloud. This 
suggests that in typical primordial star forming conditions secondary 
star formation can actually take place. 

Machacek, Bryan & Abel (2001) find that the fraction of gas avail- 
able for star formation in Pop III objects of mass AI exposed to a flux 
with intensity Jlw in the Lyman-Werner band is ^ 0.06 ln(M/M(;i), 
where the mass threshold, M^^, is given by: 

(Mth/MQ) = 1.25 X 10^ + 8.7 x 10^ ( "^^^ A . (16) 

VIO ^^erg ^cm ^Hz / 

An AMR simulation by O'Shca & Norman (2007b) confirms that the 
formation of small mass, primordial objects in the presence of a Lyman- 
Werner flux is delayed (although not completely suppressed because 
enough II2 is retained in the core), and, as a result, the virial mass of 
the halo at the time of the collapse is increased compared to when 
no dissociating flux is present (see Fig. 23). The same problem has been 
analyzed by Susa Sz Umemura (2004a) by means of three-dimensional 
SPH calculations, where radiative transfer is solved by a direct method 
and the non-equilibrium chemistry of primordial gas is included. They 
find that star formation is suppressed appreciably by UVB, but baryons 
at high-density peaks are self-shielded, eventually forming some amount 
of stars. Similarly, a 1-D radiation-hydrodynamical simulation by Ahn 
&; Shapiro (2007) suggests that minihalos down to total masses of 
10^ M0 are able to self-shield against dissociating radiation from a 
nearby massive Pop III star and form dense cold clouds at their center, 
also due to shock-induced molecule formation. Their results though 
suffer from unphysical geometric effects due to the 1-D configuration 
and are not reproduced by the 2-D simulations by Whalen et al. (2007), 
who find that shock-induced molecule formation plays no role. Never- 
theless, they confirm that a LW flux does not exert a strong feedback 
on collapse of minihalos. 

The negative feedback described above could be counterbalanced by 
the positive feedback of H2 re-formation, e.g. in front of H 11 regions. 
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Figure 23. Mean halo quantities for several simulations with the same cosmic re- 
alization but a range of Lyman- Werner molecular hydrogen photodissociating flux 
backgrounds (O'Shea & Norman 2007). Panel (a): halo collapse redshift vs. 3lw- 
Panel (b): halo virial mass vs. collapse rdshift. Panel (c): halo virial mass vs. 3lw- 
Panel (d): virial temperature vs. 3lw- The J21 = "control" results are shown as 
an open square (it is at log 3lw = —24.5 in the panels which are a function of 
Jlw). In the bottom panel, the dashed lino corresponds to the fitting function for 
threshold mass by Machacek, Bryan & Abel (2001). 



inside relic H 11 regions, once star formation is suppressed in a halo 
and ionized gas starts to recombine (Ricotti, Gnedin & Shull 2001, 
2002; Mashchenko, Couchman & Sills 2006; Abel, Wise & Bryan 2007; 
Johnson, Grcif & Bromm 2007), in cooling gas behind shocks produced 
during the ejection of gas from these objects (Ferrara 1998). Ricotti, 
Gnedin & Shull (2002) study the formation/dissociation of H2 by means 
of hydrodynamics simulations coupled with a radiative transfer scheme 
following the propagation of ionizing and dissociating radiation. Their 
main result is that the above positive feedback is usually able to coun- 
terbalance the effect of H2 dissociation. As a consequence, the formation 
of small mass galaxies is not suppressed, as found by other authors. 
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Figure 24- Mass weighted radial profiles around the position of a star with mass 
= 100 M0. Panel A: electron number fraction (solid) and H2 mass fraction 
(dashed) for 5 different times, 0, 0.1, 1, 2.7, and 8 Myr after the star is born. The 
star dies at 2.7 Myr. Panel B: temperature. Panel C: gas density. Panel D: radial 
velocity. (Abel, Wise & Bryan 2007). 

Cen (2005) has also proposed that supersonic colhsions of cold atomic 
gas clouds, which were originally minihalos and accreted into larger 
halos, can shock heat and compress the gas, which then cools and 
forms H2 molecules. This would enhance the star formation process. 
Thus, a second burst of star formation might take place also in the 
small objects where it has been suppressed by H2 dissociation. H2 
production could also be promoted by an X-ray background, which 
would increase the fractional ionization of protogalactic gas. Such a 
positive feedback, though, is not able to balance UV photodissociation 
in protogalaxies with T^r < 2000 K (Haiman, Abel & Rees 2000; Glover 
& Brand 2003). Kuhlen & Madau (2005) have simulated with the code 
ENZO the effect of the X-ray radiation produced by the first quasar 
on the surrounding medium. They find that the net effect of X-ray is 
to reduce the gas clumping with the consequence of suppressing gas 
infall at overdensities 5 < 2000, but at higher overdensities molecular 
cooling is increased, although this is not enough to overcome heating 
in the proximity of the quasar. Similar arguments apply to cosmic rays 
(e.g. Shchekinov & Vasiliev 2004; Jasche, Ciardi & Enfilin 2007; Stacy 
& Bromm 2007). 

An example of positive feedback can be seen in Figure 24. 
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Photoheating filtering 



Cosmic reionization might have a strong impact on subsequent galaxy 
formation, particularly affecting low-mass objects. In fact, the heating 
associated with photoionization causes an increase in the temperature 
of the IGM gas which will suppress the formation of galaxies with 
masses below the Jeans mass. As Gnedin (2000b) has pointed out, one 
can expect that the effect of reionization depends on the reionization 
history, and thus is not universal at a given redshift. More precisely, one 
should introduce a "filtering" scale, kp, (or, equivalently, filtering mass 
Mp) over which the baryonic perturbations are smoothed as compared 
to the dark matter, yielding the approximate relation = S^rn^''^ . 
The filtering mass as a function of time is related to the Jeans mass 



Note that at a given moment in time the two scales can be very different. 
Also, in contrast to the Jeans mass, the filtering mass depends on the 
full thermal history of the gas instead of the instantaneous value of the 
sound speed, so it accounts for the finite time required for pressure to 
influence the gas distribution in the expanding universe. The filtering 
mass increases from roughly 10^ Mq at ^; « 10 to about 10^ at red- 
shift z ~ 6, thus efficiently suppressing the formation of objects below 
that mass threshold. Of course such result is somewhat dependent on 
the assumed reionization history. 

An analogous effect is found inside individual H ii regions around 
the first luminous sources. Once an ionizing source turns off, its sur- 
rounding H II region Compton cools and recombines. Nonetheless, the 
"fossil" H II regions left behind remain at high adiabats, prohibiting 
gas accretion and cooling in subsequent generations of Pop III objects 
(Oh & Haiman 2003). A similar entropy fioor could be obtained by 
heating the gas with Lya photons (Ciardi h Salvaterra 2007) or x-ray 
photons (Oh h Haiman 2003), although the latter mechanism would 
also promote H2 formation which in turn would increase the cooling 
(Kuhlen & Madau 2005). The suppression of minihalos formation in 
photoionized regions is increased if the clustering of the sources is 
considered (e.g. Kramer, Haiman h Oh 2006). O'Shea et al. (2005), 
on the other hand, find that it is possible for a primordial star to form 
within an H II region (similar results are found by Nagakura &: Omukai 
2005 by means of 1-D hydrodynamical calculation which include also 
the effect of HD). By means of an Eulerian adaptive mesh refinement 
simulation they find that the enhanced electron fraction within the 
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H II region catalyzes H2 formation that leads to faster cooling in the 
subsequent star forming halos, although the accretion rate is much 
lower. A similar result is found by Abel, Wise & Bryan (2007) who, 
via a hydrodynamic simulation with radiative transfer, follow the for- 
mation of a first star and its H 11 region, and the consequent collapse 
inside it of a second star. They also argue that an entropy floor in 
relic H 11 regions is present, but not as high as discussed by Oh &; 
Haiman (2003), because the typically high density gas does not get 
very heated. This is confirmed by Mesinger, Bryan & Haiman (2006), 
who use a hydrodynamic simulation to study the collapse of minihalos 
inside H 11 regions. They find that there exist a critical intensity for the 
UV flux such that for lower values H2 formation is enhanced, while for 
higher star formation is reduced. This critical value is a function of the 
gas density. 

Haiman &; Bryan (2006) interpret the low value of the Thomson 
scattering optical depth inferred by the 3-yr (confirmed by the 5-yr) 
WMAP release as an evidence for the efficiency of radiative feedback on 
small mass object formation. In fact, the ionizing photon production 
efficiency in these objects must have been low in order to not overpro- 
duce the optical depth. Nevertheless, a comprehensive conclusion on 
the impact of radiative feedback effects is still impossible, given the 
large number of different assumptions, approximations and parameters 
used in the (sometimes specific) studies present in the literature. With 
the only aim to give a flavor of the large discrepancies, in Fig. 25 we 
compare various predictions on the mass of halos affected by radiative 
feedback (in some of the forms discussed above). There we show the 
results of Ciardi et al. (2000) (long dashed-dotted line), Kitayama et 
al. (2000) (dotted), Kitayama et al. (2001) (short dashed and long 
dashed), Dijkstra et al. (2004) (solid) and Yoshida et al. (2003) (short 
dashed-dotted) ; obviously, the comparison is far from exhaustive. Here 
we simply summarize the main features addressed in those papers, while 
for the details we refer the reader to them. Ciardi et al. (2000) have 
derived the minimum mass a halo must have to self-shield against an 
incident UV flux (in addition to the background, the contribution from 
the direct flux from nearby objects is included). They have studied 
the non-equilibrium multifrequency radiative transfer of the incident 
spectrum in the approximation of a homogeneous gas layer composed 
of H, H^, H+, He, He^, He^+, H2, H2^ and free electrons. Kitayama 
et al. (2000) have studied the collapse of primordial objects in the 
presence of an UV background radiation, by means of a 1-D spherical 
hydrodynamical calculation combined with a treatment of the radiative 
transfer of photons. An analogous, but more reflned calculation, has 
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Figure 25. Mass of halos affected by radiative feedback in the formulation of Ciardi 
et al. (2000) (long dashed-dotted line), Kitayama et al. (2000) (dotted), Kitayama et 
al. (2001) (short dashed and long dashed), Dijkstra et al. (2004) (solid) and Yoshida 
et al. (2003a) (short dashed-dotted). The straight dotted hne corresponds to the 
formation redshift of 3-(t peaks in a ACDM model. For details refer to the text and 
the above papers. 



been performed by Kitayama et al. (2001), for different assumptions 
for the incident flux. The curve of Dijkstra et al. (2004) is the mass 
corresponding to the case in which the presence of an external ionizing 
flux reduces the gas infall of a factor of 2. The calculations are per- 
formed by means of a 1-D spherically symmetric hydrodynamic code. 
Finally, Yoshida et al. (2003a) derive the minimum mass of halos that 
host gas clouds in their N-body/SPH simulations, which include the 
non-equilibrium evolution of the relevant species. 

In summary, the presence of a UV flux delays the collapse and cooling 
of small mass halos and the amount of cold gas available primarily 
depends on the intensity of the flux and the mass of the halo. Some 
controversy though still exists on the efficiency of such feedback. 
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4.2. Mechanical Feedback 

In addition to the radiative feedback, the mechanical feedback associ- 
ated with mass and energy deposition from the first stars, can deeply 
affect the subsequent SF process in several ways. 

Blowout and Blowaway 

Depending on the mass and the dark matter content of galaxies, SN 
and multi-SN events might induce partial {blowout) or total (blowaway) 
gas removal from the galaxy itself, regulating the SF (e.g. Mac Low & 
Ferrara 1999; Nishi & Susa 1999; Ferrara & Tolstoy 2000; Springel & 
Hernquist 2003; Bromm, Yoshida & Hernquist 2003; Silk 2003; Wada & 
Venkatesan 2003; Greif et al. 2007; Whalen et al. 2008). It is found that 
only objects with masses M £ 5 x 10^ Mq can experience a complete 
blowaway, while larger objects have lower mass ejection. Thus, quench- 
ing SF in galaxies by ejecting large fractions of their gas is very difficult. 
These results have been substantiated and confirmed by Mori, Ferrara 
&; Madau (2002) who have performed a three-dimensional hydrody- 
namic simulation, using a nested grid method to follow the evolution of 
explosive multi-SN events in an object of total mass of 10^ M© at z = 9, 
representing a 2-a fluctuation of the density power spectrum. They find 
that, depending on the stellar distribution, less than 30% of the avail- 
able SN energy is converted into kinetic energy of the escaping material, 
the remainder being radiated away. It appears that mechanical feedback 
is less efficient than expected from simple energetic arguments, as off- 
nuclear SN explosions drive inward-propagating shocks that tend to 
collect and pile up cold gas in the central regions of the host halo (see 
Fig. 26). Thus, also relatively low-mass galaxies at early epochs may be 
able to retain a considerable fraction of their gas and continue forming 
stars. A variant of this study, using the SPH method and considering 
a single supernova explosion in Pop III objects at higher rcdshift, has 
been carried out by Bromm, Yoshida & Hernquist (2003). These au- 
thors find that the more energetic PISN explosions (M^ = 250 Mq; 
E ~ 10^^ erg) lead to a blowaway (this has been later confirmed by 
Greif et al. 2007 for a PISN with mass = 250 M© and E w lO^^ gj.g)^ 
while the lower explosion energy corresponding to a single Type II 
SN (E ^ 10^^ erg) leaves much of the halo intact. Nevertheless, the 
extent to which such feedback can operate depends sensitively on the 
initial configuration of the medium surrounding the explosion site and 
the ionization front expansion around a massive progenitor star can 
significantly aid gas evacuation by SN blastwaves (Kitayama &; Yoshida 
2005), as can be seen in Fig. 27. 
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Figure 26. Snapshots of the simulated logarithmic gas number density due to SN 
explosions (at five different elapsed times) in a halo with mass M = llfh-^MQ at 
redshift 2 = 9 (Mori, Madau & Ferrara 2002). The three panels in each row show 
the spatial density distribution in the X — Y plane on the nested grids. The left, 
middle, and right panels in each row correspond to the spatial resolution levels. 
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Figure 27. Eventual fate of a halo given the indicated explosion energy. The first 
letter refers to the final state of the halo prior to the explosion (E: photoevaporated; 
P: partly ionized, defined as the I-front not reaching the virial radius; N: neutral, 
or a failed H ii region). The second letter indicates outcome of the explosion (D: 
destroyed, or F: fallback). See Whalen et al. (2008) for details. 

Following the ejection of gas, H2 formation is favored in cooling gas 
behind the shocks, possibly balancing H2 dissociation, as discussed in 
the previous Section (Ferrara 1998). This process should be accessible 
to observations, as the re-formed H2 emits mid-IR roto- vibrational lines 
(as for example the restframe 2.2 //m) detectable by JSWT (Ciardi &; 
Ferrara 2001). 

In some cases, SN explosions could promote star formation. In fact, 
following an explosion, the interstellar gas is swept up and a dense 
shell is formed immediately behind the radiative shock front at ra- 
dius Rsh- When the dense shell gets dynamically overstable, it breaks 
into fragments, some of which will contract and eventually form stars 
(see e.g. Vishniac 1983; McCray & Kafatos 1987; Nakano 1998; Tsu- 
jimoto, Shigeyama Sz Yoshii 1999; Reynoso Sz Mangum 2001; Mackey, 
Bromm & Hcrnquist 2003; Salvaterra, Ferrara & Schneider 2003). Shell 
expansion tends to suppress (stabilize) the growth of density pertur- 



review_08.tex; 28/04/2008; 12:20; p. 69 



70 



bations owing to stretching, hence counteracting the self-gravity pull. 
The instantaneous maximum growth rate is: 



Rsh 




1/2 

(18) 



where Cg^gh = {kTgh/ fimH)^^"^ is the sound speed in the shell and po the 
ambient gas density. Instability occurs only if a; > 0, or 

Rsh 1 -irGpoRsh tcross 

where tcross ~ Rsh/cs,sh is the crossing time in the shell and t// is the 
free-fall time. Furthermore, if the fragment mass is not much larger than 
the Jeans mass, the oo > condition translates into Rgh/Rsk < ^/tff- 
Instability occurs only if Rsh/ Rsh < tc-oss/t'jf, where i?,^/, is the shell 
radius, tcross is the crossing time in the shell and tjf is the free-fall 
time. Hence, large shell velocity-to-radius ratios inhibit the formation 
of gravitationally unstable fragments. Under some conditions (see e.g. 
Salvaterra, Ferrara & Schneider 2003) the above instability criterion is 
satisfied and (typically low-mass) stars are formed as a result of shell 
fragmentation. 

Impinging Shocks 

The formation of a galaxy can be inhibited also by the effects of 
shocks from neighboring objects, causing heating, evaporation and/or 
stripping of the baryonic matter (Scannapieco, Ferrara & Broadhurst 
2000). In the former scenario (mechanical evaporation), the gas in a 
forming galaxy is heated above its virial temperature by the shock. The 
thermal pressure of the gas then overcomes the dark matter potential 
and the gas expands out of the halo, preventing galaxy formation. 
Only if the cooling time of the collapsing cloud is shorter than its 
sound crossing time, will the gas cool before it expands out of the 
gravitational well and continue to collapse. In the latter scenario, the 
gas may be stripped from a collapsing perturbation by a shock from a 
nearby source. In this case, the momentum of the shock is sufficient to 
carry the gas with it, emptying the halo of its baryons and preventing 
a galaxy from forming. 

In practice, the short cooling times for most dwarf-scale collapsing 
objects suggest that the baryonic stripping scenario is almost always 
dominant. This mechanism has the largest impact in forming dwarfs in 
the ^ 10^ Mq range, which is sufficiently large to resist photoevapora- 
tion by UV radiation, but too small to avoid being swept up by nearby 
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Figure 28. Predicted evolution of Pop 11 (solid lines) and Pop III (dashed) cosmic 
star formation rates, and mass-averaged metallicity (dot-dot-dashed) from Torna- 
tore, Ferrara & Schneider (2007). The results of three simulation runs at different 
resolution are shown for each quantity. As a reference, low-redshift measurements 
(points) taken from Hopkins (2004) are reported. 

outflows. Sigward, Ferrara & Scannapieco (2005) studied this problem 
by means of numerical simulations, finding that, if the forming galaxy 
is virialized, the impinging shock has a negligible effect on its evolution, 
while, if it is in the turnaround stage, up to 70% of its gas content can 
be stripped away. 

4.3. Chemical Feedback 

The concept of chemical feedback is relatively recent, having been first 
explored by Bromm et al. (2001) (following the early study by Yoshii & 
Sabano 1980) and Schneider et al. (2002), and subsequently discussed 
by Schneider et al. (2003), Mackey, Bromm & Hernquist (2003) and 
others. 
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According to the scenario outlined in Sec. 3, the first stars forming 
out of gas of primordial composition might be very massive, with masses 
« 10^ - 10^ Mq. The ashes of these first supernova explosions pollute 
with metals the gas out of which subsequent generations of low-mass 
Pop II/I stars form, driving a transition from a top-heavy IMF to a 
"Salpeter-like" IMF when locally the metallicity approaches the critical 
value Zar = W'^^^Zq (Schneider et al. 2002, 2003). The uncertainty 
in the value of the critical metallicity is due to the role of dust cooling, 
which has been only recently quantitatively assessed by Schneider et al. 
(2006b). In brief, metals depleted onto dust grains enable fragmentation 
to solar or sub-solar mass scales already at metallicities Z = 10~^ Zq; 
on the contrary, in the absence of dust grains, even at Z = O.OIZq 
metals diffused in the gas-phase lead to fragment mass scales which are 
still of the order of 100 solar masses. 

Thus, the cosmic relevance of Pop III stars and the transition to a 
Pop II/I star formation epoch depends on the efficiency of metal enrich- 
ment from the first stellar explosions, the so-called chemical feedback, 
which is strictly linked to the number of Pop III stars that explode 
as PISN, the metal ejection efficiency, transport and mixing in the 
IGM. Depending on the strength of the chemical feedback Pop III star 
formation could last for a very short time and their contribution to 
the ionizing photon production could be negligible (see e.g. Ricotti Sz 
Ostriker 2004a for a discussion). It is very likely though that the tran- 
sition occurred rather smoothly because the cosmic metal distribution 
is observed to be highly inhomogeneous: even at moderate redshifts, 
2; 3, the clustering properties of C iv and Si iv QSO absorption 
systems are consistent with a metal filling factor < 10%, showing that 
metal enrichment is incomplete and inhomogeneous (see Sec. 5). 

As a consequence, the use of the critical metallicity as a global 
criterion is somewhat misleading because chemical feedback is a local 
process, with regions close to star formation sites rapidly becoming 
metal-polluted and overshooting Zcr, and others remaining essentially 
metal-free. Thus, Pop III and Pop II star formation modes could have 
been coeval, and detectable signatures from Pop III stars could be 
found well after the volume-averaged metallicity has become larger than 
critical. 

In a seminal paper Scannapieco, Schneider Sz Ferrara (2003) stud- 
ied, using an analytical model of inhomogeneous structure formation, 
the separate evolution of Pop Ill/Pop II stars as a function of star 
formation and wind efficiencies. Their main finding was that, essen- 
tially independent of the free parameters of the model. Pop III stars 
continue to contribute appreciably to the star formation rate density 
at much lower redshift. These findings have been now confirmed by 
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additional numerical work (Tornatore, Ferrara k, Schneider 2007) and 
are summarized in Fig. 28, where the relative contribution of Pop III 
and Pop II to the cosmic star formation is shown along with the mean 
IGM metallicity evolution. Due to inefficient heavy element transport 
by outflows and slow "genetic" transmission during hierarchical growth, 
large fluctuations around the average metallicity arise; as a result Pop 
III star formation continues down to z = 2.5, but at a low peak rate 
of 10~^ M0yr~^Mpc~^ occurring at z 6 (about 10~^ of the Pop II 
one). This finding has important implications for the development of 
efficient strategies for the detection of Pop III stars in primeval galaxies, 
as discussed previously. A similar result has been found by Wyithe &; 
Cen (2007). 



5. Cosmic Reionization and IGM Metal Enrichment 

Although cosmic reionization has received great attention in the last 
decade due to the availability of high-z absorption line studies and 
progresses in the CMB experiments, the nature of the ionizing sources, 
the reionization history and its effect on structure formation remain 
unclear and highly debated. This is mostly due to uncertainties in the 
modeling of several physical issues: properties of first stars and quasars, 
ionizing photon production and radiative transfer, just to mention a 
few. Whatever their nature, only a fraction of the emitted ionizing 
photons is able to escape from their production site and reach the IGM. 
Thus, the knowledge of the escape fraction is crucial for any theoretical 
modeling of the reionization process. The same stars that produce the 
ionizing radiation axe expected to pollute the IGM with heavy elements. 
In the following we will discuss the present understanding of the above 
processes. 

5.1. Escape Fraction 

The first attempt to derive a value of the escape fraction, /e^cj for a 
Milky Way-type galaxy dates to Dove & Shull (1994), who assumed 
a smoothly varying H I galactic distribution and OB associations lo- 
cated in the Galactic plane. By integrating the fraction of Lyman 
continuum (Lyc) photons that escapes the disk over the adopted lumi- 
nosity function of OB associations, they estimated /g^c ~ 14%. Dove, 
Shull & Ferrara (2000) later improved the calculation by solving the 
time-dependent transfer problem of stellar radiation through evolving 
superbubbles. Their main result is that the shells of the expanding su- 
perbubbles quickly trap ionizing photons, so that most of the radiation 
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escapes only after the superbubbles have broken out of the Galaxy, 
when however the ionizing power of the central stellar chister has al- 
ready largely faded away. This results in an escape fraction roughly 
a factor of 2 lower than the one obtained by Dove & Shull (1994), 
although the exact value depends on the star formation history. An 
even more refined approach is the one by Fujita et al. (2003), in which 
the effects of repeated supernova explosions, including the formation 
and evolution of superbubbles, are modeled using a hydrodynamical 
simulation. They confirm that the shells can trap ionizing radiation 
very effectively until the bubbles start to accelerate, causing the shells 
to fragment. The values found for fesc ai^e roughly consistent with the 
Dove, Shull Sz Ferrara (2000) estimates. All the above studies model 
the idealized case of a smoothly varying H i distribution and sources 
positioned on the Galactic plane. A more realistic case would give a 
higher value of the escape fraction. The effects of gas density inho- 
mogeneities have been studied by Giardi, Bianchi &; Ferrara (2002) 
for a Milky Way-type galaxy. To this aim, a comparison between a 
smooth Gaussian distribution and an inhomogeneous fractal one has 
been made, including realistic assumptions for the position and emis- 
sion properties of the ionizing stellar sources based on the available 
data, and the 3-D radiative transfer of ionizing photons. While a fractal 
interstellar medium results in an escape fraction roughly constant over a 
wide range of ionization rates, in a Gaussian distribution f^sc decreases 
with decreasing rates. Based on a model that provides a good fit to the 
observed size distribution of H i holes in nearby galaxies, Clarke & Gey 
(2002) derive a simple relationship between the star formation rate and 
the interstellar porosity of a galaxy. This gives a critical star formation 
rate for which the porosity is of the order of 1. The authors expect high 
escape fractions only in galaxies whose star formation rates exceed the 
critical value, as for example Lyman Break Galaxies (LBGs). 

Additional theoretical works (Ricotti & Shull 2000; Wood & Loeb 
2000; Fujita et al. 2003) have extended the analysis to include high- 
redshift galaxies, for which the escape of Lyc photons can be modified 
by geometrical effects (i.e. disk vs. spheroidal systems) and by the 
higher mean galactic Interstellar Medium (ISM) density. Ricotti &: Shull 
(2000) use models of smoothly distributed gas in hydrostatic equilib- 
rium and estimate fesc by integrating the escaping radiation of each 
single OB association over the luminosity function of the OB associa- 
tions. They perform a parametric study to understand the dependence 
of fesc on redshift, mass, star formation efficiency, stellar density dis- 
tribution and OB association luminosity function. Generally speaking, 
the escape fraction decreases as the object virialization redshift or mass 
become larger or the rate of Lyc photons emitted smaller. Their analysis 
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Figure 29. Left panel: Escape fraction as a function of the ionization rate for a 
clumpy (circles) and a Gaussian (asterisks) density distribution (Ciardi, Bianchi & 
Ferrara 2000). Central panel: Redshift evolution of the escape fraction for different 
masses of the host dark matter halo (Ricotti & ShuU 2000). Right panel: Evolution 
of the escape fraction as a function of the mass of the first, massive star (Yoshida 
et al. 2007). 



is confined to objects with M ~ 10 Mq, but a similar result is 
found by Wood & Loeb (2000) for objects with M « 10^-^^ j^q. They 
improved on the previous calculation by following the propagation of 
ionizing photons with a 3-D radiation transfer code and considering 
also a case of a two-phase medium. They find that the escape fraction 
is enhanced in a clumpy ISM. Although the derived escape fraction is 
lower than 1% at z ~ 10 for stellar sources, for miniquasars, which have 
higher luminosities and harder spectra, fesc can be > 30% at the same 
redshift. The higher escape fraction expected for quasars suggests that, 
if sufficiently abundant, they might have played a crucial role in the 
reionization of the IGM (see next Section). A similar trend with the 
mass and the virialization redshift of the galaxies is found for single 
metal- free very massive stars (Kitayama et al. 2004), by means of 
1-D hydrodynamic code coupled with radiative transfer. In contrast 
with the above works, Fujita et al. (2003), in the previously mentioned 
paper, find f^sc > 20% at z > 5. This reversed trend would be due 
to the fact that ionizing photons escape more easily from high-redshift 
disks because the gas distribution is highly stratified under the strong 
gravitational potentials of the compact, high-redshift halos, and this 
allows quick acceleration of the superbubbles. The same trend is found 
by Razoumov & Sommer-Larsen (2006) who find that fesc increases 
from 1-2% at z = 2.39 to 6-10% at z = 3.6 by postprocessing high- 
resolution simulations of galaxies with radiative transfer. A value of 
few % for galaxies with M > 10^^ M© and SFR 1 - 5 yr"^ 
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over the redshift range 3 — 9 is found by Gnedin, Kravtsov &; Chen 

(2008) by means of AMR simulations with on-thc-fly radiative transfer. 
Escape fractions larger than 70% as high as 100% (depending on the 
mass of the star and the hosting halo) are instead expected from very 
massive stars at high redshift (Whalen, Abel k. Norman 2004; Alvarez, 
Bromm &: Shapiro 2006; Yoshida ct al. 2007). Thus, the majority of 
theoretical models converge on values fesc <15% (the number though 
depends e.g. on mass of the host object, redshift, gas distribution), but 
a clear consensus has not yet been reached, although the first, relatively- 
or very-massive stars seem to have f^sc >70% and as high as 100%. 
For an example of the dependence of the escape fraction see Fig. 29. 

5.2. Hydrogen Reionization 

It is commonly believed that the reionization process proceeds through 
different stages. According to the terminology introduced by Gnedin 
(2000a), during the initial, pre-overlap stage, the hydrogen throughout 
the universe is neutral except for isolated H ll regions due to individ- 
ual sources. Following is the overlap stage, during which individual 
H II regions overlap and reionize the low-density, diffuse intergalactic 
gas. Some neutral hydrogen remains in dense clumps which are then 
slowly reionized during the post-overlap stage. The epoch of complete 
reionization, Zianj is somewhat dependent on the definition of "com- 
plete". A possible definition could be the redshift at which the mean 
free path to the ionizing radiation is of the order of the Hubble radius, 
or the time derivative of the mean free path has a peak. A more prac- 
tical definition, associated with a quantity more easily measurable, is 
that Zjon corresponds to the epoch when the volume weighted neutral 
hydrogen fraction becomes < 10^'^ Increasing attention has been ded- 
icated to the theoretical modeling of the reionization process, adopting 
semi-analytical (from the pioneering Shapiro & Giroux 1987; Fukugita 
&; Kawasaki 1994; Miralda-Escudc & Rccs 1994 and Tegmark, Silk 
&L Blanchard 1994 to the more recent Haiman & Loeb 1997; Madau, 
Haardt & Rees 1999; Valageas & Silk 1999; Cojazzi et al. 2000; Chiu 
k Ostriker 2000; Miralda-Escude, Haehnelt & Rees 2000; Wyithe k 
Loeb 2003a; Cen 2003a; Liu ct al. 2004; Furlanetto, Zaldarriaga k 
Hernquist 2004a; Furlanetto k Oh 2005; Benson et al. 2006; Choudhury 
k Ferrara 2006; Cohn k Chang 2007) perturbative or semi-numerical 
(e.g. Zhang, Hui k Haiman 2007; Mesinger k Furlanetto 2007) and 
numerical approaches (e.g. Gnedin k Ostriker 1997; Gnedin 2000a; 
Ciardi et al. 2000; Razoumov et al. 2002; Ciardi, Stoehr k White 2003; 
Ciardi, Ferrara k White 2003; Ricotti k Ostriker 2004a; Sokasian et 
al. 2003; Iliev et al. 2006a; Trac k Cen 2007; Croft k Altay 2008). Two 
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main ingredients are required for a proper modeling of the reionization 
process: (i) a reliable model of galaxy formation and (ii) an accurate 
treatment of the radiative transfer of ionizing photons. 

Despite many applications of hydrodynamical simulations to the 
structure formation process, much of our current understanding comes 
from semi- analytical models, based on simplified physical assumptions. 
These models have the advantages of allowing a larger dynamic range 
and of being fast enough to explore a vast range of parameters. Their 
main disadvantage is that no information on the spatial distribution 
of structures is provided. For this reason a new method has recently 
been developed that combines N-body simulations of the dark matter 
component, with semi-analytical models predicting galactic properties 
"a posteriori" (e.g. Kauffmann ct al. 1999). The limitation of the above 
methods is related to the resolution of the N-body simulations, which 
turns out to be critical in studies of the reionization process. In fact, the 
resolution must be high enough to follow the formation and evolution 
of the objects responsible for producing the bulk of the ionization 
radiation. At the same time, a large simulation volume is required 
to have a region with "representative" properties and to avoid biases 
due to cosmic variance. It should be noted, however, that although 
the contribution of small mass objects to the ionizing photon produc- 
tion might be negligible due to feedback effects, they are nevertheless 
sinks of ionizing radiation (e.g. Barkana &: Loeb 2002 and references 
therein). Thus, any numerical simulation which is not able to resolve 
these objects is likely to underestimate the gas clumping factor, and 
consequently the recombination rate, resulting in an earlier reionization 
epoch. These minihalos introduce a substantial cumulative opacity to 
ionizing radiation, unless they are photoevaporated by UV radiation 
quickly enough (Barkana &: Loeb 2002; Shapiro, Iliev & Raga 2004; 
see previous Section). Ciardi et al. (2006) demonstrated by means of 
numerical simulations of cosmic reionization including minihalos, that, 
depending on the details of the minihalo formation process, their effect 
can delay complete reionization by Az ~ — 2. 

At the present state, while the evolution of dark matter structures is 
well understood, the treatment of the physical processes that govern the 
formation and evolution of the luminous objects (e.g. heating/cooling 
of the gas, star formation, feedback) is more uncertain. In particular, no 
study treats the entire range of the feedback effects that can influence 
the galaxy formation process (see previous Section). A first attempt 
to treat self-consistently a number of feedback effects ranging from the 
mechanical energy injection to the H2 photodissociating radiation pro- 
duced by massive stars has been done by Ciardi et al. (2000), followed 
by the more recent Ricotti, Gnedin & Shull (2002) and Yoshida et al. 
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(2003a), but the complexity of the network of processes makes their 
implementation not trivial. 

The next challenge for this kind of study is to develop accurate and 
fast radiative transfer schemes that can be then implemented in cosmo- 
logical simulations. The full solution of the seven dimension radiative 
transfer equation (three spatial coordinates, two angles, frequency and 
time) is still well beyond our computational capabilities and, although 
in some specific cases it is possible to reduce its dimensionality, for 
the reionization process no spatial symmetry can be invoked. Thus, 
an increasing effort has been devoted to the development of radiative 
transfer codes based on a variety of approaches and approximations 
(e.g. Umemura, Nakamoto Sz Susa 1999; Razoumov & Scott 1999; Abel, 
Norman & Madau 1999; Gnedin 2000a; Ciardi et al. 2001; Gnedin & 
Abel 2001; Sokasian, Abel & Hcrnquist 2001; Ccn 2002; MascUi, Fcrrara 
& Ciardi 2003; Whalen & Norman 2006; Mellema et al. 2006a; Qiu et 
al. 2006; Rijkhorst et al. 2006; Susa 2006; Pawlik & Schaye 2008; for a 
comparison between different codes see Iliev et al. 2006b). The first at- 
tempt to include the radiative transfer in simulations of the reionization 
process has been done by Gnedin Sz Ostriker (1997). Ideally, radiative 
transfer codes should be eventually coupled with a hydrodynamic sim- 
ulation. Although such an implementation has already been attempted 
by Gnedin (2000a), the Local Optical Depth approximation adopted for 
radiative transfer is less accurate for optical depths > 1 and therefore 
may not be suitable to describe in detail the entire reionization history. 
A similar implementation, with the Optically Thin Variable Edding- 
ton Tensor approximation for radiative transfer, is used by Ricotti & 
Ostriker (2004a). Although in some respect these two works improve 
on previous calculations (they couple radiative transfer with hydrody- 
namic simulations and include some feedback effects), their conclusions 
might be affected by possible uncertainties introduced by the small 
cosmic volumes (maximum box size of 4 Mpc) considered (see e.g. 
Barkana &; Loeb 2004 for a discussion of box dimension in simulations of 
cosmic reionization). Iliev et al. (2006a) have in fact clearly shown that 
only boxes larger than 20-30/i~^ Mpc can properly describe the global 
reionization process, while smaller boxes exhibit a large scatter. For this 
reason, boxes with such dimensions are needed (see Trac & Cen 2007 
for a self-consistent calculation on large boxes). The largest simulations 
run (applied to the study of high-z QSOs spectra) have 1280/t~''" Mpc 
side (Kohler, Gnedin & Hamilton 2007), but the resolution is very poor 
and results from smaller simulations are used to include the effect of 
scales below the resolution limit. 

An illustrative example of models of the reionization process is 
shown in Fig. 30. The Figure refers to the simulations by Ciardi, Stoehr 
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Figure 30. Roioiiizatioii maps showing the rcdshift evolution of the number density 
of neutral hydrogen in the field (L = 20/i~^ Mpc comoving; left panels) and in a 
protocluster environment (L = 10/i~^ Mpc comoving; right panels). From top to 
bottom the redshift of the simulations is 15.5, 12 and 9 (Ciardi, Stoehr & White 
2003; Ciardi, Ferrara & White 2003). 



& White (2003), employing a combination of high resolution N-body 
simulations, a semi-analytical model of galaxy formation and the radia- 
tive transfer code CRASH (Ciardi et al. 2001; Maselli, Ferrara & Ciardi 
2003). The maps are indicative of how the reionization process is af- 
fected by the environment: although in a protocluster environment the 
ionizing photon production per unit mass is higher at high redshift 



review_08.tex; 28/04/2008; 12:20; p. 79 



80 



than in a field region, as the high-density gas, more common in a 
protoclustcr, is more difficult to ionize and recombines much faster, 
filaments of neutral gas are still present after the field region is almost 
completely ionized. 

The reionization process is affected not only by the environment, 
but also by the scale considered. In fact, while on scales of order of 
few Mpc reionization might be complete, on larger scales it could still 
be quite inhomogeneous due to the large scatter in the abundance of 
star-forming galaxies at early epochs (Barkana & Loeb 2004). A useful 
mathematical tool to study such complex topology is constituted by 
the Minkowski functionals (Gleser et al. 2006) of isodensity surfaces 
which can be used to discriminate among different reionization his- 
tories. In general, we expect that the properties of the H ii regions 
depend on many poorly constrained quantities. McQuinn ct al. (2007) 
have investigated this issue by means of radiative transfer simulations 
on large scales, including small scale physics via analytic prescriptions. 
They find that the morphology of H II regions is similar for fixed val- 
ues of the mean ionization fraction, with other parameters playing a 
secondary role (the next most important dependence is on the na- 
ture of the ionizing sources, with larger and more spherical regions 
for rarer sources). Clearly topology of the H ii regions, clustering of 
sources and the dumpiness of the IGM are key issues to understand 
reionization in detail. These have been investigated by Furlanetto k. 
Oh (2005). Their main result is that the radiation background at any 
point evolves through a series of discrete jumps until it saturates at 
the mean free path, when a region becomes recombination-limited; 
afterwards the evolution proceeds very slowly. Towards the end of 
reionization fiuctuations in the mean free path within large contiguous 
H II regions dominates over the clustering of the sources. However, Liu 
et al. (2006) noted that if the UVB evolves rapidly the following facts 
can be explained without invoking a fluctuating UVB at z ^ Q: rapid 
increase and large scatter in the Gunn-Peterson optical depth, long-tail 
distribution of the transmitted flux, long dark gaps in spectra. 

The nature of the sources responsible for the IGM reionization is still 
the subject of a lively debate, but most theoretical models adopt stellar 
type sources, including a variety of spectra and IMF. These have to be 
constrained on the basis of the available observations: (i) the value of 
the IGM temperature from the Lya forest at z ~ 2 — 4 (e.g. Bryan 
& Machacek 2000; Ricotti, Gnedin k Shull 2000; Schaye ct al. 2000b; 
McDonald et al. 2001; Meiksin, Bryan k Machacek 2001; Zaldarriaga, 
Hui k Tegmark 2001); (ii) the abundance of neutral hydrogen at 2; > 6 
from the spectra of high-redshift QSOs or galaxies (e.g. Fan et al. 2000, 
2001, 2003, 2006) and (Hi) the measurements of the Thomson scattering 



review_08.tex; 28/04/2008; 12:20; p. 80 



81 



optical depth (e.g. Kogut et al. 2003; Spergel et al. 2003; Page et al. 
2007; Nolta et al. 2008). A simple scenario in which the IGM reionized 
early {zion > 10, as suggested by the last condition) and remains so 
thereafter, is ruled out as the IGM would reach an asymptotic thermal 
state too cold compared to observations (Theuns, Schaye Sz Haehnelt 
2000; Hui & Haiman 2003; Wyithe & Loeb 2003b). 

The most attractive scenario requires an enhanced ionizing photon 
emission at high redshift (compared to the low redshift), obtained either 
through a higher ionizing photon production (metal-free stars and/or 
top-heavy IMF), star formation efficiency or escape fraction. A suitable 
combination of these parameters allows an early enough reionization 
and matches the measured Thomson scattering optical depth (e.g. Cen 
2003b; Ciardi, Fcrrara & White 2003; Ciardi, Stoehr & White 2003; 
Haiman & Holder 2003; Onkcn & Miralda-Escude 2004; Sokasian et 
al. 2003, 2004; Wyithe & Loeb 2003b; Iliev et al. 2006a; but see also 
Ricotti & Ostriker 2004a). Following feedback and metal enrichment, 
it is plausible that the ionizing photon production drops, resulting 
in a partial/total recombination of the IGM, followed by a second 
reionization at z 6, produced by more standard sources (Wyithe &; 
Loeb 2003a; Cen 2003a; but see also Oh & Haiman 2003 who propose 
that the ionizing photon production is regulated by gas entropy rather 
than by a transition from Pop III to Pop II stars) or more simply 
resulting in a temporary slight decrease of the ionizing flux (Tumlinson, 
Vcnkatcsan & Shull 2004). Such reionization history might also help in 
releasing the tight constraint posed by the IGM temperature at z 3 
on the reionization redshift. The contribution to cosmic reionization 
from very massive mctal-free stars is still very uncertain. In fact, to 
produce a number of ionizing photons sufficient to substantially ionize 
the IGM, there would be an over production of metals (e.g. Ricotti Sz 
Ostriker 2004a; Rozas, Miralda-Escude &: Salvador-Sol'e 2006) and the 
transition from Pop III to Pop II stars would occur early on (Fang & 
Cen 2004). It has been proposed that while very massive metal- free 
stars would form in metal- free halos with Tyir < 10^ K, larger metal- 
free halos woTild host metal- free stars with masses < 100 Mq. In this 
case, it is found (Greif & Bromm 2006) that, due to the impact of 
radiative feedback, the contribution of very massive Pop III stars to 
reionization (and metal enrichment, see Section 5.4) is not substantial 
as their formation terminates at relatively high redshift. The formation 
of lower mass Pop HI stars would stop at lower redshift due to metal en- 
richment and feedback, and their contribution to reionization would be 
more important. This picture is supported by Schneider et al. (2006a). 
Using a semi-analytic model to follow the formation and evolution of 
dark matter halos, coupled with a self-consistent treatment of chemical 
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Figure 31. The best-fitting model Clioudiiury & Ferrara (2006). Tlie panels show 
the redshift evolution of: (a) the volume-averaged electron and HI fraction. The 
arrows show an observational lower limit from QSO absorption lines at z = 6 and 
upper limit from Lya emitters at z — 6.5; (b) the cosmic star formation history, 
with the contribution of Pop III and Pop II stars, (c) the number of source counts 
above a given redshift, with the observational upper limit from NICMOS HUDF; 
(d) the electron scattering optical depth, with observational constraints from 3-yr 
WMAP data; (e) Lya effective optical depth; (f ) Ly/3 effective optical depth; (g) the 
evolution of Lyman-limit systems; (h) photoionization rates for neutral hydrogen, 
with estimates from numerical simulations; (i) temperature evolution of the mean 
density IGM. For tha data sources, see the original paper. 



and radiative feedback, they find that the constraints imposed by the 
NICMOS observations (Bouwens et al. 2005) can be fulfilled together with 
those imposed by WMAP only if Pop III stars have a slightly top-heavy 
IMF, but masses < 100 Mq. A similar result is found by Daigne et 
al. (2004) and Choudhury &; Ferrara (2006). The last paper presents a 
comprehensive model of reionization which satisfies simultaneously all 
the available constraints set by different observations (see caption of 
Fig. 31 for a complete list). Although many details might be modified 
by future observations, a solid conclusion seems to be that reionization 
has been a gradual and long-lasting event which was started by Pop 
III stars at z ~ 15, it was 90% complete at z ~ 9 and completed only 
with the help of Pop II stars and quasars at z ~ 7. 

Following the first year release of WMAP data, that implied a high 
value for the Thomson scattering optical depth, several alternative 
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processes/contributions to enhance the high redshift ionizing photon 
emission have been proposed. For example, a component of a non-scalc- 
free isocurvature power spectrum in addition to the more standard 
scale-free adiabatic power spectrum, would induce an earlier structure 
formation and yield an efficient source of early star formation, down to 
z ~ 10, without violating the existing constraints from the Lyman-a 
forest observations (Sugiyama, Zaroubi &; Silk 2004). The same effect 
could be obtained with non-Gaussian density fluctuations (Chen et al. 
2003) . If magnetic fields generated in the early universe had a strength 
B > 0.5 nG (consistent with BBN and CMB constraints that set an 
upper limit to the comoving amplitude of ~ lOnG at 1 Mpc scale), 
they would as well induce additional density perturbation on small 
scales (Tashiro & Sugiyama 2006a). A contribution could also come 
from globular clusters, for which a high escape fraction is expected (Ri- 
cotti 2002). An enhancement of the ionizing photon production caused 
by an increased/earlier formation of low mass galaxies, though, would 
crucially depend on the efficiency of feedback effects and their ability to 
partially/completely suppress star formation in primordial objects (see 
previous Section). Moreover, e.g. Ricotti, Gnedin h Shull (2002) have 
shown that even for the case of weak feedback, small galaxies can only 
partially reionize the IGM. More specifically, reionization begins earlier 
with the inclusion of small galaxies (with the consequence of raising the 
value of the Thomson scattering optical depth), but their influence is 
limited by Jean-mass filtering in the ionized regions and the time of 
overlap is dictated by the efficiency of the higher mass halos (Iliev et 
al. 2007a). In more exotic scenarios, loops formed from a cosmic string 
network would induce structure formation at earlier times (Glum &; 
Vilenkin 2006). 

In addition to stellar-type sources, a contribution to the UV pho- 
ton budget could also come from other sources. In particular, high- 
z quasars, if present, can produce a substantial amount of ionizing 
photons (e.g. Madau, Haardt & Rees 1999; Valageas & Silk 1999; 
Miralda-Escude, Haehnelt & Rees 2000; Wyithe & Loeb 2003a; Ripa- 
monti, Mapelli, Zaroubi 2008), thanks to their higher luminosity and 
escape fraction. Madau et al. (2004), for example, have shown that 
miniquasars powered by intermediate-mass black holes (the remnants 
of the first generation of massive stars) could produce a significant 
amount of ionizing photons at z > 15, under a number of plausible 
assumptions for the amount of gas accreted onto the black holes and 
their emission spectrum. On the other hand, as the recombination 
rate and hence the required emissivity to reionize the IGM, increase 
both with redshift and with the luminosity of the sources (Miralda- 
Escude, Haehnelt & Rees 2000), this might reduce the effectiveness 
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of reionization by luminous quasars. Another crucial question about 
QSOs is their relatively short lifetime, which might effectively limit 
the growth of their ionized regions. In addition, Salvaterra, Haardt & 
Ferrara (2005) find that, unless primordial quasars are extremely X-ray 
quite (to avoid exceeding the observed soft X-ray background), their 
contribution to reionization is only secondary. 

Also more energetic photons, e.g. X-ray photons, can contribute 
to the IGM reionization (Oh 2001; Venkatesan, Giroux & Shull 2001; 
Madau et al. 2004; Ricotti & Ostriker 2004b). This component could 
come from early quasars, but also from thermal emission from the hot 
supernova remnants, or inverse Compton scattering of soft photons by 
relativistic electrons accelerated by supernova explosions (Oh 2001). As 
X-rays have much larger mean free paths, compared to photons from 
stellar spectra (the mean free path for a 1 keV X-ray photon is 10^ times 
larger than for a 13.6 eV photon), and their escape fraction is ~ 1, they 
can permeate the IGM relatively uniformly. But only if redshifted X- 
rays are taken into account, their contribution to reionization can be 
substantial (e.g. Ricotti & Ostriker 2004b). 

Another possible contribution to reionization at high redshift could 
come from decaying particles and neutrinos. While massive active neu- 
trinos (Sciama 1990; Scott, Rees & Sciama 1991) have been excluded by 
current cosmological data (see Spergel et al. 2003 for the latest results), 
a decaying sterile neutrino does not violate existing astrophysical lim- 
its on the cosmic microwave and gamma ray backgrounds and could 
partially ionize the IGM. While Hansen &: Haiman (2004) find that 
a decaying sterile neutrino could account for the observed Thomson 
scattering optical depth, Mapelli & Ferrara (2005) conclude that it 
must have played only a minor role in the reionization process. But 
there is no lack of other decaying particle-physics candidates; e.g. cryp- 
tons, R-parity violating gravitinos, moduli dark matter, superheavy 
dark matter particles, axinos and quintessinos. Chen & Kamionkowski 
(2004) have studied these particles, deriving the channels in which most 
of the decay energy ionizes and heats the IGM gas and those in which 
most of the energy is instead carried away (e.g. photons with energies 
100 keV < -E < 1 TeV). They find that decaying particles can indeed 
produce an optical depth consistent with the WMAP results, but they 
produce new fluctuations in the CMB temperature-polarization power 
spectra. The amplitude of such fluctuations generally violates current 
constraints for decay lifetimes that are less than the age of the universe, 
while it is usually consistent with the data if the lifetime is longer (but 
see also Kasuya & Kawasaki 2004, who flnd that also particles in a short 
lifetime are consistent with the WMAP data. In a subsequent version 
of the paper they favor longer lifetimes to be consistent the WMAP3 
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results). Avelino & Barbosa (2004) have shown that a contribution to 
reionization from the decay products of a scahng cosmic defect network 
could help to reconcile a high optical depth with a low redshift of com- 
plete reionization. Finally, it has been proposed (Sethi &; Subramanian 
2005) that the primordial magnetic field energy could dissipate into 
the IGM by ambipolar diffusion and by generating decaying magneto 
hydrodynamics turbulence. These processes can modify the thermal 
and ionization history of the IGM and contribute to the Thomson 
scattering optical depth. 

Neither X-rays nor decaying particles alone produce a fully ionized 
IGM, but the IGM may have been warm and weakly ionized prior to 
full reionization by stars and quasars. This scenario would also alleviate 
constraints on structure formation models with low small-scale power, 
such as those with a running or tilted scalar index, or warm dark matter 
models, which, alone, would not provide enough ionizing photons (e.g. 
Haiman & Holder 2003; Somerville, Bullock & Livio 2003; Yoshida et 
al. 2003b; Yoshida et al. 2003c). 

5.3. Helium Reionization 

The stellar sources that ionize hydrogen can as well produce singly 
ionized helium, through photons with energies hv > 24.6 eV. Moreover, 
as these two species have comparable recombination rates, their reion- 
ization history is very similar. On the other hand. He II reionization 
requires harder photons (hi^ > 54.4 eV), cither from massive Pop III 
stars (Venkatesan, Tumlinson &: Shull 2003), from quasars (Madau, 
Haardt & Rees 1999; Miralda-Escude, Haehenelt & Rees 2000; Wyithe 
& Loeb 2003a; Madau et al. 2004; Ricotti & Ostriker 2004b) or from 
thermal emission from shock-heated gas in collapsed cosmic structures 
(Miniati et al. 2004). This last contribution would be comparable to the 
QSO one at z ~ 3 and dominating at z > 4, in terms of He ii ionizing 
photons. Since H l and He I do not absorb a significant fraction of pho- 
tons with energy higher than 54.4 eV, the problem of He ii reionization 
can be decoupled from that of the other two species. 

The He II Lya absorption is generally much stronger than the H I Lya 
absorption, by a factor ?7=N(He ii )/N(H i )r~^ 4rHcii/THi, where the 
second equality is valid for optically thin lines. The larger strength 
of He II arises because it is harder to photoionize than H i (although 
the number of helium atoms is smaller than hydrogen atoms by a 
factor of 10), owing to lower fluxes and cross sections at its ion- 
izing threshold. In addition. He iii recombines 5.5 times faster than 
H II and He II (this number becomes ~ 5.9 when including the increase 
in the electron density due to the ionization of helium). This suggests. 
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supported by semi-analytical (Madau, Haardt & Rees 1999; Miralda- 
Escude, Haehenelt &: Rees 2000; Venkatesan, Tumlinson & Shull 2003; 
Wyithe & Loeb 2003a; Furlanctto & Oh 2008) and numerical (Sokasian, 
Abel & Hernquist 2001; Gleser et al. 2005; Pasclios et al. 2008) calcula- 
tions, that helium reionization occurred at z < 5, i.e. after completion of 
hydrogen reionization. The known population of quasars is sufficient to 
completely reionize He ii before 2; ~ 3 (e.g. Miralda-Escude, Haehenelt 
& Rees 2000), although this is subject to the assumed gas clumping and 
a contribution from high-z quasars (or other sources) could be needed 
(e.g. Madau, Haardt & Rees 1999). On the other hand, one might 
worry that the primordial contribution to He II ionizing photons may be 
too high and induce an almost simultaneous H i and He 11 reionization. 
Assuming a standard power-law spectrum for miniquasars, it is shown 
that only a small delay between the complete overlapping of H II and 
He III regions is expected. Thus if in addition to miniquasars Pop III 
stars and thermal emission from shock-heated gas also contributed to 
He II reionization in the early universe, one might expect a double 
reionization, a first time by the above sources, followed by recombi- 
nation as Pop n stars with softer spectra took over, and then by a 
second overlap phase at 2; « 3, driven by the known population of 
quasars (e.g. Venkatesan, Tumlinson & Shull 2003; Wyithe & Loeb 
2003a) . Ricotti & Ostriker (2004b) find that reionization of He 11 from 
primordial black holes can be almost complete at z ^ 15. In their 
model, the full reionization at z ~ 3 is due to a hardening of the 
background radiation rather than QSO activity. The He iii ionizing flux 
should fluctuate substantially because, in a random region of space, it is 
dominated by only a few active sources. This can produce optical depth 
fluctuations in the He II forest on scales comparable to the mean free 
path, « 7000 km s^^. These fluctuations may also be due, however, to 
the random fluctuations in the number of underdense voids in the IGM 
crossed by the line of sight. A strong proximity effect is also expected, 
and it is observed in the spectra of quasars for which the He 11 Lya 
forest has been observed. 

Associated with IGM reionization is also an increment of its temper- 
ature (e.g. Miralda-Escude & Rees 1994; Theuns et al. 2002a; Ricotti & 
Ostriker 2004b). Because its cooling time is long, the low-density IGM 
retains some memory of when and how it was reionized. The post- 
ionization temperature is generally higher for harder spectra of the 
ionizing radiation, but it depends also on the type of source. While the 
IGM can be heated up to T ~ 10, 000 K during hydrogen reionization, 
it can reach T !=a 20, 000 K during helium reionization (see also Fig. 40). 
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5.4. IGM Metal Enrichment 

The same stars that ionize the IGM are responsible for its metal en- 
richment. Differently from their present counterparts, Pop III stars do 
not suffer strong stellar winds and thus the main contribution to the 
metal pollution comes from the later stages of their evolution. Wind- 
driven mass loss is believed to be metallicity-dcpendent, and a scaling 
law oc Z^/^ has been suggested for hot stars (Kudritzki 2000; Nugis 
Sz Lamers 2000). This scaling breaks down at Z < 10"^ Zq, where 
the power-law becomes steeper (Kudritzki 2002). Marigo, Chiosi & 
Kudritzki (2003) have studied the mass loss by stellar winds for stars 
with initial masses in the range 120 — 1000 M©, including the two 
main mass loss driving processes, i.e. radiative line acceleration (e.g. 
Kudritzki 2002; Krticka & Kubat 2006) and stellar rotation (Macdcr 
&; Meynet 2000). According to their calculations, wind mass-loss for 
stars with metallicity of 10~^ Zq is significant only for very massive 
objects (750 - 1000 Mq). For a 1000 M© non-rotating star, 17.33% 
(0.03%) of the initial mass is ejected via mass loss as helium (oxygen) 
during the stellar lifetime. A small amount of carbon and nitrogen is 
ejected as well. Several important effects are introduced if rotation 
is taken into consideration. First, chemical mixing is more efficient 
for lower metallicity for a given initial mass and velocity. This favors 
the evolution into the red supergiant stage, when mass loss is higher, 
and it enhances the metallicity of the surface of the star, boosting the 
radiatively driven stellar winds. Second, although metal-poor stars lose 
less mass to winds compared to their more enriched counterparts, they 
also lose less angular momentum, so that they have a larger chance of 
reaching the break-up limit (when the outer layers become unbound 
and are ejected) during the main sequence phase. Meynet, Ekstrom 
&; Maeder (2006), including the above effects of rotation, find that a 
60 M0 star, with Z = 10~^ and initial velocity of 800 km/s can lose 22% 
of its mass. At Z = though the mass loss is of the order of few percent 
(Ekstrom, Meynet &: Maeder 2006). Finally, whereas pulsation-driven 
mass loss is important, if not dominant, at high metallicity, Baraffe, 
Heger & Woosley (2001) have shown that it is negligible at very-low 
metallicity. If indeed mass loss in very massive stars is dominated by 
optically thick, continuum-driven outbursts or explosions as suggested 
by Smith &; Owocki (2006), nitrogen would be expelled, but no other 
heavy metal because these stars are not in advanced O, Si and subse- 
quent burning stages. We can conclude then that, as metal pollution 
from stellar winds is negligible for very-low metallicity stars, the main 
contribution to the IGM metal enrichment com,es from those stars that 
end up their lives as PISN and core-collapse SNe (see Fig. 16). 
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Figure 32. Mean mctallicity vs. ovcrdonsity for three different simulations at z = 3 
and for a simulation at 2 = (sec details in Aguirre et al. 2001c). The top axis 
is log Nhi of an absorber corresponding to the bottom-axis overdensity at 2 = 3. 
The vertically striped box outlines the approximate current constraints from the 
Lya forest at ^; « 3, while the horizontally striped box shows an extension of these 
constraints to lower densities. 



As metal factc^ries (stars) are almost always associated with high 
density regions, in the absence of any efficient diffusion/transport mech- 
anism, one would then expect a very strong density/metallicity corre- 
lation. Although such positive correlation is seen, it is not as strong 
as expefcted on the basis of the above naive assumption. However, 
a relation is predicted by simulations, as is clearly seen in Fig. 32 
(Aguirre et al. 2001c). In any case, the observed large-scale clustering 
of metal absorbers encodes valuable information about the masses of 
the objects from which they were ejected. Likewise, as the maximal 
extent of each enriched region is directly dependent on the velocity at 
which the metals were dispersed, measurements of the small clustering 
of these absorbers is likely to constrain the energetics of their sources. 

Different mechanisms have been suggested to remove metals from 
the galaxies in which they are produced, ranging from dynamical en- 
counters between galaxies, to ram-pressure stripping, supernova-driven 
winds or radiation-pressure driven dust cffliix. While dynamical re- 
moval undoubtedly occurs at some level, it is not clear if it can account 
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for the observed metallicity in the z 3 IGM (Aguirre et al. 2001c; but 

sec Gncdin 1998, who finds merging events as a dominant mechanism 
for transporting heavy elements into the IGM). The same argument 
apphes to the ram-pressure stripping, which has a minor impact on the 
overall metallicity and filling factor as it occurs only in the densest and 
most polluted regions of space (Scannapieco, Ferrara & Madau 2002). 

Metal ejection by galactic winds is one of the most popular mecha- 
nisms for the IGM pollution (e.g. Gen &; Ostriker 1999; Aguirre et al. 
2001b; Madau, Ferrara & Rees 2001; Scannapieco, Ferrara & Madau 
2002; Qian & Wasserburg 2005; Scannapieco 2005; Fieri, Martel & 
Grenon 2006; Tornatore, Ferrara &: Schneider 2007). Madau, Ferrara 
Sz Rees (2001) find that pre-galactic outflows from the same primordial 
halos that reionize the IGM, could also polhitc it to a mean metal- 
licity of Z ^ 10~'^ Zq. The advantage of metals produced in these 
low-mass halos is that they can more easily escape from their shallow 
potential wells than those at lower redshift. In addition, the enriched 
gas had to travel much shorter distances between neighboring halos at 
these early times, and it might therefore have been easier to obtain 
non- negligible cosmic metal filling factors (see also Scannapieco 2005). 
Moreover, the velocities associated with these pre-galactic outflows are 
small enough to leave the thermal and structural properties of the IGM 
unperturbed, in contrast with stronger galaxy winds (with velocities 
Z 200 — 300 km s"^), which, although capable of enriching the IGM 
to the mean level observed, are likely to overly disturb it (e.g. Gen & 
Ostriker 1999; Aguirre et al. 2001b). Rayleigh- Taylor instabilities at the 
interfaces between the dense shell which contains the swept-up material 
and the hot, metal-enriched low-density bubble may contribute to the 
mixing and diffusion of heavy elements. The volume filling factor of 
the ejecta can be up to 30% or higher, depending on the star forma- 
tion efhciency, with the majority of the enrichment occurring relatively 
early, 5 ^ z £ 14 (Madau, Ferrara & Rees 2001; Scannapieco, Ferrara 
& Madau 2002; Thacker, Scannapieco & Davis 2002), and a possibly 
large contribution from very massive stars (Qian & Wasserburg 2005). 
Such early enrichment model can explain quite well the observed red- 
shift evolution of elemental abundances derived from QSO absorption 
line experiments, as seen from Fig. 33, where model predictions are 
compared with available data. Nevertheless, enrichment by outflowing 
galaxies is likely to have been incomplete and inhomogeneous. Op- 
penheimer & Dave (2006), via hydrodynamic simulations, find that 
momentum-driven winds from galaxies at z < 5 are able to reproduce 
the measured G iv over the range 2; ~ 2 — 5. Thus it seems that a very 
early generation of stars is not required to explain these observations. 
The worry remains that the hydrodynamical motions induced by the 
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Figure 33. Metallicity evolution for low-density (12 < log Nhi < 15 ) IGM. Points 
are the data of Songaila (2001); dotted, solid, and dashed lines are the predicted 
values for star formation efficiencies 0.5, 0.1, and 0.01, respectively (see Scannapieco, 
Ferrara & Madau 2002 for details). 

outflow may spoil the nice agreement of Lya absorption line statistics 
derived from studies in which these phenomena were not included. This 
problem has been convincingly washed out by the results of Bertone 
&: White (2006), who, by combining N-body simulations and semi- 
analytical descriptions of the winds, find that winds do not sensibly 
alter the statistical properties of the Lya forest. 

But if indeed, as discussed in Sec. 3, the first stars were in the range 
for PISN and were produced, in the absence of feedback effects, in halos 
of masses of M ~ 10^ — 10^ M©, the more energetic explosions would 
expel ^ 90% of the stellar metals into a region of ~ 1 kpc (Bromm, 
Yoshida & Hernquist 2003). Due to this burst-like initial SF episode, a 
large fraction of the universe could have been endowed with a metallic- 
ity floor Z Zq already at z ^ 15. Similarly, Matteucci &: Calura 
(2005) find that ^ 110 PISN are enough to pollute a cubic comoving 
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Mpc to such metallicity floor. The authors flnd that, although it is easy 
for PISN to enrich the gas at high rcdshift, a contribution from Pop II 
stars is needed to reproduce the metalhcity patterns observed in Lya 
forest and in damped Lya systems. 

An alternative mechanism for the enrichment of the IGM is based 
on the ejection of dust grains, within which heavy elements are locked, 
by radiation pressure, along the lines suggested by Ferrara et al. (1991). 
Aguirre et al. (2001a) and Bianchi &; Ferrara (2005), have investigated 
the IGM pollution caused by such a mechanism. An attractive feature of 
this model is that, unlike winds, enrichment by dust would not impact 
the thermal/structural properties of the IGM, as no shock waves are 
involved in the transport. According to this study, dust pollution can 
account for the observed level of C and Si enrichment of the z ~ 3 
IGM, but not for all the chemical species observed in clusters at low 
redshift. Thus, the authors suggest a possible hybrid scenario in which 
winds expel gas and dust into galaxy halos but radiation pressure dis- 
tributes the dust uniformly through the IGM. Of course, it is still to 
be demonstrated that grains are efficiently sputtered in the IGM and 
therefore can release their metals. 

As discussed in Sec. 3, the presence of dust, as well as metals, 
can deeply affect the primordial star formation process. Thus, it is 
important to determine the dust production in the early universe. 
Recently, this issue as been discussed by Nozawa et al. (2003) and 
Schneider, Ferrara & Salvaterra (2004), which investigated the forma- 
tion of dust grains in the ejecta of PISN and core-collapse SNe. The 
fraction of mass locked into dust grains is 2 — 5% of the progenitor 
mass for SN II and 15 — 30% for PISN; however, PISN dust depletion 
factors (fraction of produced metals locked into dust grains) are smaller 
(< 40%) than SN II ones. These conclusions depend very weakly on 
the thermodynamics of the ejecta, which instead affect considerably 
the grain formation epoch, composition and size distribution (Nozawa 
et al. 2003). Thus, if the first generation of stars were very massive, a 
large amount of dust grains might be produced in the early universe. 

An additional problem is that our cosmic inventory of the cosmic 
metals associated with the observed star formation history is lacking a 
consistent fraction of heavy elements. This discrepancy is often referred 
to as the missing metals problem (see Pettini 2006 for the most recent 
review on the subject). One possibility is that missing metals escaped 
detection because they reside in the (hot) halos of the parent galaxies, 
as suggested by Ferrara, Scannapieco & Bergeron (2005) who found 
that up to 95% of the produced metals could be hidden in such a 
hot gas phase. Consistently with this hypothesis, Scannapieco et al. 
(2006) find that the correlation function of C iv absorbers favor a metal 
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enrichment mode in which metals are confined within bubbles with a 

typical radius of 2 Mpc around galaxies of mass ~ 10^2 M0. Note that 
this is not in contradiction with the fact that these metals could have 
been injected by early progenitors of these galaxies, as predicted by the 
pre-enrichment scenario. 

5.5. Key Observations 

5.5.1. Escape Fraction 

Observationally, a wide range of values for the escape fraction has been 
deduced, but it appears that - broadly speaking - relatively low values 
of fesc are favored. For example, most of the detections of starburst 
galaxies with HUT^^ (Hopkins Ultraviolet Telescope) and FUSE^^ (Far 
Ultraviolet Spectroscopic Explorer) (Leitherer et al. 1995; Hurwitz, 
Jelinsky 8z Dixon 1997; Heckman et al. 2001) are consistent with f^sc < 
10%, although objects have been observed with /e^c < 57% (Hurwitz, 
Jelinsky & Dixon 1997). Tumlinson et al. (1999), based on Ha images, 
derive a value of the escape fraction lower than 2%. However, absorption 
from undetected interstellar components could allow the true escape 
fractions to exceed these upper limits. Values of f^sc < 10% arc found 
from the detection of flux beyond the Lyman limit in a composite 
spectrum of 29 LBGs at z = 3.4 (Steidel, Pettini & Adelberger 2001; 
Haehnelt et al. 2001). An f^sc ^ 20% would be compatible also with 
the observed UVB (Bianchi, Cristiani &: Kim 2001). Fernandez-Soto, 
Lanzetta &; Chen (2003) found that, in a sample of 27 spectroscopically 
identified galaxies of redshift 1.9 < z < 3.5 in the Hubble Deep Field, 
on average no more than 4% of the ionizing photons escape galaxies 
and even when all systematic effects are included, the data could not 
accommodate any escape fraction larger than 15%. A high value of 
fesc has been measured only by Bland-Hawthorn & Maloney (1999; see 
Erratum 2001), who used optical line emission data for the Magellanic 
Stream to derive fesc ~ 45%. For a compilation of the existing measure- 
ments see Inoue, Iwata & Deharveng (2006). In conclusion, observations 
of the escape fraction, with few exceptions, give fesc < 15%. 

5.5.2. Gunn-Peterson Test 

The discovery of quasars at z > 5.8 (e.g. Fan et al. 2000, 2001, 2003, 
2004; Goto 2006) arc finally allowing quantitative studies of the high- 
redshift IGM and its reionization history. In particular, the detection of 
a Gunn-Peterson trough (Gunn Sz Peterson 1965) in the Keck^'^ (Becker 

http:/ /praxis. pha.jhu.cdu/hut. html 
http://fuse.pha.jhu.edu 
^° http://www2.keck.hawaii.edu 
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Figure 34- Evolution of the effective Lya optical depth. Points correspond to data. 
The dashed line shows the best-fit at z < 5.5 (Fan et al. 2006), the solid line is 
calculated from the lognormal distribution of Lya optical depths. For details see 
Becker et al. (2007). 

et al. 2001) and VLT^^ (Very Large Telescope) (Pentericci et al. 2002) 
spectra of the Sloan Digital Sky Survey quasars SDSS J1030-0524 at 
z = 6.28, and SDSS J1148+5251 at z = 6.37 (White et al. 2003) 
might indicate that the universe is approaching the reionization epoch 
at z ~ 6, consistent also with the dramatic increase in the neutral 
hydrogen optical depth at z > 5.2 observed by Djorgovski et al. (2001). 
However, while the theoretical optical depth distribution commonly 
used to measure the H i ionization rate gives a very poor fit to the 
observations, more accurate models of the density distribution (i.e. 
an inverted temperature-density relation, see Sec. 2) provide a better 
agreement and suggest that the measured dramatic increase might be 
explained other than with a reionization at z ~ 6 (Becker, Ranch & 
Sargent 2007). Analytical fits to the measured effective optical depth 
are shown in Fig. 34. 

These observations suggest that the IGM is in the post-overlap 
stage (Barkana 2002). It should be noted, though, that an analysis of 
the same high redshift quasars which takes into account the observed 
dimension of the ionized region surrounding them, implies a neutral 



http:/ /www. eso.org/projects/vlt 
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fraction of > 0.1, much higher than previous lower hmits (Wyithe 8z 
Loeb 2004; Mesinger & Haiman 2004; Wyithe, Loeb k Carilli 2005; but 
see Oh & Furlanetto 2005 for a discordant view). The main uncertainty 
in the above estimates is the quasars hfetime. Also analyzing the QSOs 
and modeling their spectra separately, Mesinger Sz Haiman (2007) find 
a lower limit of 0.033. Maselli et al. (2007), using a combination of 
SPH simulations and radiative transfer and analyzing mock spectra 
along lines of sight through the simulated QSO environment, show that 
the observed dimension of the ionized region is typically 30% smaller 
than the physical one. This induces an overestimate of the neutral 
fraction, which they find to be < 0.06. This value is consistent with that 
found by Yu Sz Lu (2005) using a semi-analytic model to describe the 
QSO's environment (including a clumping factor and the contribution 
to ionization from stars). A similar result has been found by Bolton & 
Haehnelt (2007a) for a variety of quasars lifetime and luminosity and 
including the treatment of He in their radiative transfer simulation. The 
bottom line is that, according to these studies, the size of the observed 
ionized regions cannot put stringent limits on IGM ionized fraction. 
Stronger constraints can be obtained considering the relative size of 
the ionized regions in the Lya and Ly,(^ (Bolton & Haehnelt 2007b). 
An additional complication follows from the location of QSOs in biased 
regions. In fact, as the extent of the H ii regions produced by nearby 
galaxies can be substantial, their contribution should be included when 
using the H II regions from QSO to determine the abundance of H I in 
the IGM (e.g. Alvarez & Abel 2007; Lidz et al. 2007). 

By itself, the detection of a Gunn-Peterson trough does not uniquely 
establish the fact that the source is located beyond the reionization 
epoch, since the lack of any observed flux could be equally caused by 
ionized regions with some residual neutral fraction, individual damped 
Lya absorbers or line blanketing from lower column density Lya forest 
absorbers. On the other hand, for a source located at a redshift beyond 
but close to reionization, the Gunn-Peterson trough splits into disjoint 
Lyman-a, -/? and possibly higher Lyman series troughs, with some 
transmitted flux in between them, that can be detected for sufficiently 
bright sources and used to infer the reionization redshift (Haiman & 
Loeb 1999). Alternatively, metal absorption lines, rather than H ones, 
have been proposed to probe the pre-reionization epoch (Oh 2002). To 
this goal, one should select absorption-line probes which are still unsat- 
urated when the IGM is predominantly neutral (this is possible for ions 
that are less abundant than H i or have smaller oscillator strengths), 
have ionization potentials similar to that of H, and that should lie 
redward of the Lya wavelength (to avoid confusion with lower-rcdshift 
Lya forest). This suggested experiment remains speculative as it is not 



review_08.tex; 28/04/2008; 12:20; p. 94 



95 



clear how it would be possible to have a highly neutral, but metal 
polluted IGM. 

A recently suggested technique to study reionization is through 
the statistical analysis of QSO absorption spectra. In particular, the 
dark gap width distribution, together with the complementary leaks 
distribution (Fan et al. 2006; Gallerani, Choudhury & Ferrara 2006; 
Gallerani et al. 2007; Feng et al. 2008; Liu, Bi & Fang 2008) could help 
in discriminating between a late and an early reionization history, as 
for the former >30% of the line of sights to QSOs in the redshift range 
5.7 < z < 6.3 arc expected to have dark gaps with widths > 50A, 
while none is expected for the latter (Gallerani, Choudhury &; Ferrara 
2006). In order to discriminate between different reionization histories 
though, a sample at higher redshift is required (Gallerani et al. 2007). 

5.5.3. Lya Emitters 

Lya emission has been crucial in determining the redshift of distant 
galaxies (e.g. Hu et al. 2002; Kodaira et al. 2003; Rhoads et al. 2003; Hu 
et al. 2004; Pello et al. 2004; Kurk et al. 2004; Nagao et al. 2004; Stern 
et al. 2005; Wang, Malhotra & Rhoads 2005; Wilhs & Courbin 2005; 
Kashikawa et al. 2006), with the highest z Lya emitter being at z = 6.96 
(lye et al. 2006; Ota et al. 2007; Stark et al. 2007b have found two 
possible gravitationally lensed candidates at z « 10) and it represents 
an invaluable tool to probe the high redshift universe and, possibly, the 
sources responsible for the IGM reionization. As Lya is a resonant tran- 
sition with a large cross section, even a small hydrogen neutral fraction 
in the intervening absorber can cause large opacity to Lya photons 
from high-z sources. In particular, only the blue side of the Lya line 
would be scattered away, while the red side would be observed, with a 
profile that depends on the ionization state of the IGM surrounding the 
galaxy or the presence of galactic winds (Miralda-Escude 1998; Santos 
2004). As a result of the absorption, the Lya photons are scattered 
until they redshift out of resonance, resulting in a compact halo of 
Lya light surrounding the source of the Lya emission (Loeb &; Rybicki 
1999). Although the brightness of such diffuse Lya component is very 
weak, this radiation is highly polarized and thus in principle detectable 
(Rybicki & Loeb 1999). If, on the contrary, the source is surrounded 
by a large H ii region, the Lya photons can avoid absorption, provided 
they redshift out of resonance before they reach the boundary of the 
H II region. Madau & Rees (2000) and Cen & Haiman (2000) (see also 
Cen 2003c) show that a sufficiently bright quasar placed at z > 6 can 
produce a Stromgren sphere large enough to render its Lya emission 
line detectable, e.g. by JWST or Keck. Haiman (2002) has revisited this 
problem for a fainter source, showing that even for a galaxy as the one 
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discovered by Hu et al. (2002) a significant fraction of the emission line 
can remain observable. Moreover, if one considers the contribution to 
the formation of the H II region by neighboring galaxies, the prospect 
for observing a high- 2; galaxy are even better (Furlanetto, Hernquist &: 
Zaldarriaga 2004) , also prior to complete reionization (Wyithe k. Loeb 
2005). By means of numerical simulations of cosmological reionization, 
Gnedin & Prada (2004) show that it is possible to find galaxies at 
2; !=a 9 (and possibly at 2; > 10) that are barely affected by the dumping 
wing of the Gunn-Peterson absorption. Later on, Mesinger, Haiman 
& Cen (2004) have shown that it should be possible to statistically 
extract relevant parameters, including the mean neutral fraction in 
the IGM and the radius of the local cosmological Stromgren region, 
from the flux distribution in the observed spectra of distant sources. 
Malhotra & Rhoads (2006) have estimated a lower limit for the volume 
average ionization fraction of 20-50% by assigning to each observed 
Lya emitter a minimum H 11 region. The evolution in the luminosity 
function of Lya emitters observed by Kashikawa et al. (2006) from 
z = 5.7 to z = 6.5 has been interpreted as due to a rapid evolution in 
the IGM ionization state. Dijkstra, Wyithe & Haiman (2007), though, 
argue that the observations are consistent with a mild or no evolution 
in the IGM transmission and can be explained with the evolution in the 
mass function of dark matter halos hosting Lya emitters. Following an 
approach similar to that of Malhotra & Rhoads (2006) they estimate a 
volume averaged ionization fraction > 80%. 

5.5.4. H I 21 cm Line Emission/ Absorption 

It has long been known (e.g. Field 1959) that neutral hydrogen in the 

IGM and gravitationally collapsed systems may be directly detectable 
in emission or absorption against the CMB at the frequency corre- 
sponding to the redshifted H i 21 cm line (associated with the spin-flip 
transition from the triplet to the singlet ground state). Madau, Mciksin, 
& Rees (1997) flrst showed that 21 cm tomography could provide 
a direct probe of the era of cosmological reionization and reheating. 
In general, 21 cm spectral features will display angular structure as 
well as structure in redshift space due to inhomogeneities in the gas 
density field, hydrogen ionized fraction, and spin temperature. Several 
different signatures have been investigated in the recent literature: 
(i) fluctuations in the 21 cm line emission induced by the "cosmic 
web" (Tozzi et al. 2000), by the neutral hydrogen surviving reioniza- 
tion (Ciardi &: Madau 2003; Furlanetto, Sokasian &l Hernquist 2004; 
Zaldarriaga, Furlanetto k. Hernquist 2004; Furlanetto, Zaldarriaga k. 
Hernquist 2004b; He et al. 2004; Mellema et al. 2006b; Santos et al. 
2007) , and by minihalos with virial temperatures below lO'' K (Iliev et 
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Figure 35. Expected rms brightness temperature fluctuations as a function of beam 
size for a model with reionization at z « 7.5 (left panel) and at 2 ~ 12.5 (right 
panel). A fixed bandwidth At' = 1 MHz has been assumed. Each curve corresponds 
to a different emission redshift or, equivalently, observed frequency (see Ciardi & 
Madau 2003). 



al. 2002, 2003; Furlanetto & Oh 2006; Shapiro et al. 2006); (ii) a global 
feature ( "reionization step" ) in the continuum spectrum of the radio sky 
that may mark the abrupt overlapping phase of individual intergalac- 
tic H II regions (Shaver et al. 1999); (in) and the 21 cm narrow lines 
generated in absorption against very high-redshift radio sources by the 
neutral IGM (Carihi, Gnedin, & Owen 2002; Furlanetto 2006) and by 
intervening minihalos and protogalactic disks (Furlanetto & Loeb 2002; 
Furlanetto 2006). While an absorption signal would be preferable, being 
of a higher intensity, it relies on the existence of powerful radio sources 
at high redshift, which have not yet been found. A possibility would be 
the radio afterglow of GRBs or hypernovae (which could in principle be 
visible up to 2; ~ 20 — 30) , but the estimated absorption lines are very 
difficult to detect also with the next generation of radio telescopes (loka 
&: Meszaros 2004). On the other hand, a signal in emission, although 
weaker than the one in absorption, would always be present, as long 
as the IGM is not completely ionized and its temperature is above the 
one of the CMB. 

Ciardi & Madau (2003) were the first to use large scale numerical 
simulations of hydrogen reionization by stellar sources to investigate 
the 21 cm signal expected from the diffuse, neutral IGM. They found 
that the predicted brightness temperature fluctuations on arcmin scales 
have typical values in the range 5 — 20 mK on scales below 5 arcmin. 
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with the maximum corresponding to the epoch when roughly 50% of 
the IGM is ionized, in agreement with subsequent studies (see Fig. 35). 
Ilicv et al. (2002, 2003) find that the signal expected from minihalos is 
comparable to the signal from the IGM if feedback effects are assumed 
to be negligible. Shapiro et al. (2006) use numerical simulations to 
calculate the relative importance of IGM and minihalos, finding that 
the latter dominate the emission at 2; < 18, while they are negligible at 
higher redshifts. In their analysis though they neglect the effect of Lya 
pumping and feedback, which can be crucial once the first sources of 
radiation start to shine. The inclusion of such mechanisms shows that 
the contribution from minihalos is almost always lower than that from 
the diffuse IGM (Furlanetto & Oh 2006). 

More specific applications of the 21 cm signal allow to track the phys- 
ical properties of individual H II regions as they form and grow around 
reionization sources, in particular Pop III stars and QSOs (e.g. Wyithe, 
Loeb & Barnes 2005; Zaroubi & Silk 2005; Chen & Miralda-Escude 
2006; Cen 2006; Rhook k Haehnelt 2006; Liu et al. 2007; Thomas 
&; Zaroubi 2008). As the expected angular resolution is sensitive to 
scales that correspond to the transition between ionized and neutral 
regions for hard ionizing photons sources, additional information on 
their nature could come form observations of single ionized regions 
(Zaroubi & Silk 2005), also in the form of spectral dips in frequency 
space along single lines of sight (Kohler et al. 2005). This technique 
though would be compromised once ionized regions start to overlap. In 
addition, the 21 cm radiation emitted from the partially neutral IGM 
outside H 11 regions, if detected, could provide a novel way to probe the 
growth of BHs (Rhook & Haehnelt 2006). 

A key physical mechanism to ensure a detectable 21 cm signal is the 
decoupling of the CMB temperature from the spin temperature, as the 
latter regulates the 21 cm line absorption/emission. It has been shown 
(Chuzhoy & Shapiro 2006) that the exact value of the spin temperature 
depends on physical processes that usually are not included in the 
calculations, the most important being the backreaction of resonant 
scattering on the pumping radiation and the scattering by resonant 
photons other than Lya. In addition, it is usually assumed that the 
spin temperature is decoupled from the CMB temperature due to scat- 
tering with LycK photons (the Wouthuysen-Field effect; Wouthuysen 
1952; Field 1959) but its strength can be overestimated if a proper 
treatment of the Lyman series cascades is not included (Hirata 2006; 
Pritchard &: Furlanetto 2006). Moreover, before a Lya background is 
established, the bias in the galaxy distribution can induce fluctuations 
in the Lya, and thus 21 cm, flux (Chen & Miralda-Escude 2006). These 
fluctuations could be used to probe the distribution of the very first 
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galaxies (Barkana & Loeb 2005b; Wyithe & Loeb 2007). Additional 
fluctuations arise if the natural broadening of the line is considered. 
In this case in fact the optical depth in the wings causes scattering to 
happen nearer to the source, with a distribution function that depends 
on the optical depth and the frequency of the photon (Chuzhoy k, Zheng 
2007). These fluctuations though should affect the 21 cm signal only 
during a narrow range of redshift near the beginning of reionization, as 
a later times there are so many UV photons that these corrections are 
minor. 

Kuhlen, Madau & Montgomery (2006) use hydrodynamic simula- 
tions of structure formation to study the effect of X-rays on 21 cm 
line. In a simulation without sources they find that, while in regions 
with 5 < 1 decoupling is not effective, for higher density contrasts 
H-H collisions are efficient at decoupling and most of the gas is shock 
heated, producing a signal in emission. When the effect of mini-quasars 
is included, the X-ray radiation preheats the IGM to a few thousand 
kelvins and the increased electron fraction boosts both the H-H and the 
H-e collisions (see also Nusser 2005a for the efficiency of H-e collisions in 
the decoupling during the very early stages of reionization). Thus also 
low density regions in the IGM can be seen in emission. It should be 
noted though that, although X-ray are usually considered to uniformly 
penetrate the IGM, this is not true at all times and might affect the 
21 cm line signature (Pritchard k, Furlanetto 2007). 

Also the mechanism that can heat the gas above the CMB temper- 
ature is crucial for the observation of the line in emission. In addition 
to Lya and X-ray photons (e.g. Madau, Meiksin & Rees 1997; Chen 
& Miralda-Escude 2004; Pelupessy, Di Matteo & Ciardi 2007; Ciardi 
& Salvaterra 2007; Ripamonti, Mapelli & Zaroubi 2008), also UHE- 
CRs and decaying dark matter particles can heat the gas (Shchekinov 
k. Vasiliev 2007). But it's been shown that Lya scattering is not an 
efficient mechanism to heat the IGM (Chen & Miralda-Escude 2004; 
Chuzhoy & Shapiro 2006; Furlanetto k, Pritchard 2006), while X-rays 
are better, although they can uniformly heat the IGM only sometime af- 
ter the first structure formation (Pritchard k Furlanetto 2007; Zaroubi 
et al. 2007). 

Another signal detectable in the 21 cm radiation is its polarization. 
Although intrinsically unpolarized, in the presence of a magnetic field, 
the 21 cm line would show a left- and right-handed polarized component 
due to the Zeeman effect (Cooray k Furlanetto 2005). In addition, scat- 
tering between electrons produced during reionization and the 21 cm 
quadrupole would induce polarization (Babich k Loeb 2005). Such 
observations would probe the magnetic field on scale of H II regions 
around bright quasars or the intergalactic magnetic fields. An alter- 
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Figure 36. Contours of the angular power spectrum (in mK) of the the 21 cm fluctu- 
ations (shaded area) and the angular power spectra due to intensity fluctuations of 
radio galaxies (solid contours) and radio halos (dot-dashed contours) as a function of 
angle (/) and frequency. The lines in the right panel show the contours for the total 
foreground signals. In the left panel, foreground sources above a flux limit Sc ~ 0.1 
mjy have been removed (see Ciardi & Madau 2003; Di Matteo, Ciardi & Miniati 
2004). 



native method to probe such magnetic fields would be through the 
imprint they leave on the brightness temperature fluctuations (Tashiro 
& Sugiyama 2006b). 

The signatures in the brightness temperature discussed above can 
possibly be detected by the planned facilities SKA, LOFAR, 21cmA/PAST 
and MWA (see e.g. Valdes et al. 2006; Mellema et al. 2006b for the ob- 
servability of its fluctuations). The same telescopes should also be able 
to distinguish between different reionization models and discriminate, 
e.g., a history with discrete H ii regions from one with partial uniform 
ionization or from a double reionization (Furlanetto, Zaldarriaga & 
Hernquist 2004b). In principle, if peculiar velocities are associated with 
H I , they could induce non linear distortions in the redshift space and 
reduce the power spectrum of the brightness temperature. Wang & 
Hu (2006) though find that the next generation of radio telescopes are 
unlikely to be affected by such distortions. 

However, these experiments are extremely challenging due to fore- 
ground contamination from unresolved extragalactic radio sources (e.g. 
Di Matteo et al. 2002), free-free emission from the same reionizing halos 
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(e.g. Oh & Mack 2003), synchrotron emission from cluster radio halos 
and rehcs (Di Matteo, Ciardi & Miniati 2004) and the Galactic free- 
free and synchrotron emission (Shaver et al. 1999). As the foregrounds 
are slowly varying as a function of frequency, it appears that the best 
strategy for measuring the 21 cm signal would be to search for frequency 
rather than angular fluctuations (Di Matteo et al. 2002; Gnedin & 
Shaver 2004; Zaldarriaga, Purlanetto & Hernquist 2004; Santos, Gooray 
& Knox 2005; Wang et al. 2006; Gleser, Nusser & Benson 2008). Di 
Matteo, Ciardi & Miniati (2004) have used cosmological simulations 
to predict the impact of extragalactic foreground produced by both 
extended (e.g. cluster radio halos and relics) and point-like (e.g. radio 
galaxies and free-free emission from ISM) sources in a self-consistent 
way. They find that the contribution to the angular fluctuations at 
scales 9^1 arcmin is dominated by the spatial clustering of bright 
sources. Hence, efficient removal of such sources may be sufficient to 
allow the detection of angular fluctuations in the 21 cm emission free 
of extragalactic foregrounds (see Fig. 36). 

5.5.5. CMB Footprints 

In addition to the temperature anisotropies discussed in Sec. 2, reion- 
ization introduces also a polarization signal in the CMB spectrum. 
It is well known that the primordial polarization is generated at the 
recombination epoch through Thomson scattering of the quadrupole 
of the temperature anisotropy. The same mechanism operates during 
the reionization epoch and again distorts the shape of the polarization 
power spectrum. The WMAP telescope during the first year of operation 
detected a correlation between CMB temperature and polarization on 
large angular scales (Kogut et al. 2003; Spergel et al. 2003), whose 
amplitude is related to the total optical depth of CMB photons to 
Thomson scattering, giving Tg ~ 0.16 ± 0.04 (the uncertainty quoted 
for this number depends on the analysis technique employed). This 
high value strongly favored an early reionization; e.g. 0.12 < Tg < 0.20 
translates into a reionization redshift of 13 < z < 19 for a model 
of instantaneous reionization. The limits shift at lower redshift for a 
more realistic reionization, their precise value depending on the adopted 
history. In the light of the first three years of the WMAP experiment, the 
above results have been revised to Tg = 0.09 ± 0.03 (see Fig. 37; Page 
et al. 2006; Spergel et al. 2006). Such value has been only marginally 
updated by the five years WMAP results, having now converged to = 
0.087 ± 0.017 (Dunkley et al. 2008); note, however, the considerably 
smaller errorbar. The updated Thomson scattering optical depth can 
be better interpreted by a reionization redshift z ~ 10 (Spergel et al. 
2006; Dunkley et al. 2008), substantially lower than the first year value. 
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Figure 37. Comparison between the likelihood curves for r in the WMAP-lyr and 
WMAP-3yr data. Different curves refer to different analysis (see Page et al. 2006 for 
details) . 



An interesting analysis by Hansen et al. (2004) shows that the estimate 
of Te varies substantially when inferred from power spectra computed 
separately on the northern and the southern hemispheres, and that 
the full sky estimate by the WMAP team could in large part originate in 
structures associated with the southern hemispheres. In addition, the 
interpretation of the Thomson scattering optical depth could change 
if magnetic fields were present during the epoch of reionization. In 
this case, the observed anisotropies at large scales could be partially 
explained with the existence of magnetic fields (Gopal & Sethi 2005). 

Another valuable tool to study the reionization process is to detect 
the imprint it should leave on the CMB on small angular scales. When 
the IGM gets reionized, Thomson scattering off the newly formed elec- 
trons would, on one hand, blend photons initially coming from different 
directions thus damping the primary anisotropies, and on the other 
hand, create small-scale secondary anisotropies. Several groups have 
studied the CMB temperature anisotropies produced by an inhomoge- 
neous reionization (Knox, Scoccimarro & Dodelson 1998; Bruscoli et 
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al. 2000; Benson et al. 2001; Gnedin & JafFe 2001; Santos et al. 2003; 

McQuinn et al. 2006; Salvaterra et al. 2005; Zahn et al. 2005; Iliev et 
al. 2006c). In particular, calculations of anisotropies from simulations 
of patchy reionization give peak values of the power spectrum between 
«i 10-1^ (e.g. Gnedin & JafFe 2000) and few 10-^^ (e.g. Salvaterra 
et al. 2005; Iliev et al. 2006c), which might be probed by future inter- 
ferometers like ALMA, ACT^^ or CQ^^. The detection of such anisotropies 
would be an invaluable tool to discriminate between different sources of 
ionizing photons (e.g. McQuinn et al. 2006), especially in combination 
with measurements of 21 cm line emission (Salvaterra et al. 2005; Al- 
varez et al. 2006a; Holder, Iliev &: Mellema 2007; Adshead &: Furlanetto 
2007; Tashiro et al. 2008; but see also Cooray 2004), taking advantage 
of the correlation between the two signals. 

Future missions, like PLANCK'^'*, should also be able to distinguish 
between different reionization models (e.g. single, double reionization) 
based on the temperature-polarization power spectrum (e.g. Holder et 
al. 2003; Hu & Holder 2003; Kaplinghat et al. 2003; Naselsky & Chiang 
2004; Zhang, Pen k Trac 2004; Colombo et al. 2005). The overall differ- 
ences in the temperature (polarization) angular power spectra between 
prompt and extended reionization models with equal optical depths 
are less than 1% (10%) (Bruscoli, Fcrrara & Scannapieco 2002). The 
same observations should also allow to discriminate a contribution to 
reionization from decaying particles (e.g. Naselsky &z Chiang 2004). 

Analogously to free electrons, heavy elements leave an imprint on 
the CMB due to scattering in their fine-structure lines, partially erasing 
original temperature anisotropies and generating new fluctuations. The 
main difference with the Thomson scattering is that the latter gives 
equal contribution over the whole CMB spectrum, while the line scat- 
tering gives different contributions to different spectral regions. Since 
every line affects only a given range of redshift, having observations at 
different wavelengths can provide an upper limit to abundances of dif- 
ferent species at various epochs (Basu, Hernandez-Monteagudo &; Sun- 
yaev 2004; Hernandez-Monteagudo, Verde Sz Jimenez 2006; Hernandez- 
Monteagudo et al. 2007, 2008). 

It should be noted that the time and duration of recombination 
directly influence the characteristics of CMB anisotropies. Several au- 
thors (e.g. Leung et al. 2004; Dubrovich Sz Grachev 2005; Chluba &; 
Sunyaev 2006) have found that the detailed treatment of processes 
such as the physics of the H+ and Hc^ populations, the softening 
of the matter equation of state due to the transition from ionized 

http:/ /www.hep.upenn.odu/angclica/act/act.html 

http://brown.nord.nw.m/CG/CG.htm. 

http://astro.estec.esa.nl/SA-general/Projects/Planck 
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Figure 38. Predicted number of GRBs with observed flux greater than a limiting 
flux F, at different observed wavelengths (see Ciardi & Loeb 2000 for more details). 



to neutral matter, the induced emission of the softer photons which 
change the two-photon decay rate of the 2s level of H, affect the value 
of the anisotropies at large multipoles at a percent level. Thus, in this 
epoch of precision cosmology, also those processes that can affect the 
anisotropies at a percentage level, should be included in the analysis. 

5.5.6. Gamma Ray Bursts as Cosmic Lighthouses 
The observation of the high-z universe, beyond the epoch of reioniza- 
tion, could be feasible directly if bright enough sources are present. If 
indeed GRBs are related to the SFR of the host galaxy and occur 
at high redshift (the highest GRB has been detected by SWIFT at 
z = 6.29, Watson et al. 2006), their afterglow emission can be used 
to probe the ionization history of the IGM (Totani et al. 2006). In 
contrast to other sources, such as galaxies or quasars, the observed 
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infrared flux from a GRB afterglow at a flxed observed age is only 
a weak function of its redshift (e.g. Ciardi & Loeb 2000). Thus, the 
afterglow of these sources could be detected by the next generation of 
telescopes, e.g. JWST, and used to probe the SF and ionization history 
of the universe (see Fig. 38; e.g. Ciardi &: Loeb 2000; Blain k, Natarajan 
2000; Lamb & Reichart 2000; Porciani & Madau 2001; Bromm & Loeb 
2002; Choudhury & Srianand 2002; Gallerani et al. 2008; McQuinn 
et al. 2008), as well as the characteristics of the circumburst medium 
(e.g. Frail et al. 2006). A proof that very high-redshift GRBs could be 
observed is GRB 0803 19B, that would have been easily detected in K 
band with meter-class telescopes up to z « 17 (Bloom et al. 2008). 
In particular, it has been noted (Inoue, Yamazaki & Nakamura 2004) 
that as the detection of a GRB afterglow through a NIR filter i beyond 
the Lya drop-out redshift, ^Ly^out' probes the ionization state of the 
IGM around ^Ly^out' •^i ^ ^ filters would be suitable to 

probe the ionization state at 2; ?a 5 — 8, 8 — 11, 11 — 15 and 15 — 20, 
respectively. Thus, the null detection of the GRB afterglows at z « 11 
in J-band would imply a high neutral fraction at z « 8 — 11. The 
Subaru^^ telescope has sufficient resolution to do the photometry. To 
better constrain the redshift, spectroscopy is required and it will be 
only feasible with SWIFT or future missions as ASTRO-F or JWST. 

Free electrons along the light path cause distortions in the spec- 
trum and light curve of an early radio afterglow. This distortion can 
be quantified through the Dispersion Measure (DM). At z > 6 the 
DM due to the ionized IGM probably dominates the contribution from 
the Galactic plasma and the plasma in the host galaxy. Thus, DM of 
GRBs at various redshifts and directions would help in determining the 
reionization history and possibly the topology of the H II regions. Both 
SKA and LOFAR should be able to detect bright afterglows and measure 
the DM towards high redshift GRBs (loka 2003; Inoue 2004). 

5.5.7. Helium. ReionizatAon 

The Gunn-Peterson trough technique to infer the abundance of neu- 
tral hydrogen in the IGM can be extended to the determination of 
He II abundance. The first evidence of such a trough has been found 
in the spectrum of the quasar Q0302-003 (z = 3.286, Jakobsen et al. 
1994), followed by more determinations in the range 2.4 < z < 3.5 
(e.g. Davidsen, Kriss Sz Zheng 1996; Hogan, Anderson & Rugers 1997; 
Anderson et al. 1999; Smette et al. 2002; Jakobsen et al. 2003; Zheng et 
al. 2004). In particular, high resolution spectra of the quasar HE 2347- 
4342 at z = 2.885 (and more recently of QSO 1157+3143, Reimers et 

^® http://subarutelescope.org/ 
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Figure 39. Effective optical depth as a function of redshift as measured from a 
sample of 1061 QSO spectra from the SDSS database (Bernardi et al. 2003). 

al. 2005a) have shown a largely fluctuating ry, in the range 0.1 — 1000, 
with an average of 80 (e.g. Reimers et al. 1997; Kriss et al. 2001; 
Smette et al. 2002; Shull et al. 2004). This implies that the metagalactic 
radiation field exhibits inhomogeneities at the scale of ~ 1 Mpc and 
has been interpreted as evidence for a patchy He ii ionization in the 
IGM. However, it is unclear if the fluctuations could be caused by 
fluctuations of the IGM density or of the ionizing background flux 
(Miralda-Escude, Haehnelt & Rees 2000; Kriss et al. 2001; Bolton 
et al. 2006), a wide range in the QSO spectral indexes (Teller et al. 
2002), finite QSO lifetimes (Reimers et al. 2005b) or the filtering of 
QSO radiation by radiative transfer effects (Shull et al. 2004; Maselli 
& Ferrara 2005). In fact, also QSOs that are believed to be in the 
post-reionization phase (e.g. HS 1700+6416, that, together with HE 
2347-4342, is the only line of sight with He ii absorption resolved in the 
Lya forest structure) show a large variation in -q (Fechner et al. 2006). 
Particular care should be taken when estimating the r/ value. In fact, 
if the same line width is assumed for He ii and H i absorption features 
without including thermal broadening, as done in most analysis, this 
leads to a spurious correlation of the inferred r] value with H i column 
density (Fechner & Reimers 2007). 

An He ii reionization at z ~ 3 was initially favored by a claim of an 
abrupt change at that redshift in the Si iv /C iv ratio in a sample of 13 
QSO spectra, indicating that, while at z ^ 3 the ionization balance is 
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Figure 40- Temperature at the mean density (left panel) and slope of the effective 
equation of state (right panel) as a function of redshift. The dashed lines are from 
a simulation designed to fit the data. For more details see Schaye et al. (2000b). 



consistent with a pure power-law ionizing spectrum, at higher redshift 
the spectrum must be softer. An analogous result is obtained by Heap et 
al. (2000), who analyzed the Gunn-Peterson He ii absorption trough at 
z = 3.05, finding that the spectrum is much softer than the one deduced 
at z = 2.87. However, a more recent VLT/UVES study (Kim, Cristiani 
& D'Odorico 2002) using 7 QSOs finds no strong discontinuity for the 
quantity Si iv /C iv around 2; « 3 and suggests that it might not be a 
good indicator of He 11 reionization. Also the analysis of 19 high-quality 
quasar absorption spectra by Aguirre et al. (2004) seems to indicate 
no evolution in the Si IV /C IV ratio. Finally, Boksenberg, Sargent & 
Ranch (2003) find no jump at 2; ~ 3 in the median column density ratio 
Si IV /C IV , and also other ionic ratios vary continuously with redshift. 
On the other hand, the presence of a significant amount of O VI at 
z > 3 suggests either a considerable volume of He iii bubbles embedded 
in the more general region in which the ionizing flux is heavily broken, 
or the addition of collisional ionization to the simple photoionization 
models (Songaila 1998). 

As usually r]=N(H.e 11 )/N(H i )> 1 and rneii = (?7/4)thi, He 11 can 
be tracked into much lower density regions than H i , making the 
He II absorption a good probe of low-density regions or hot, coUisionally 
ionized regions in the IGM. Shull et al. (2004) find that r] is system- 
atically larger in voids (thi < 0.05), compared to filaments, possibly 
implying a softer ionizing radiation field. 

There have also been suggestions (Ricotti, Gnedin & Shull 2000; 
Schaye et al. 2000b; Theuns, Schaye & Haehnelt 2000) for a relatively 
sudden increase in the line widths between z = 3.5 and 3.0, which 
could be associated with entropy injection resulting from the reioniza- 
tion of He II . In particular, Theuns et al. (2002b) use hydrodynamical 
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simulations to show that the observed temperature increase leads to a 
relatively sudden decrease in the hydrogen Lya effective optical depth. 
Bernardi et al. (2003) find clear evidence for such a feature in the 
evolution of the effective optical depth determined from a sample of 
about 1100 quasars obtained from the SDSS (Fig. 39). However, the 
existence of the peak in the IGM temperature evolution is still quite 
uncertain, as the data are also consistent with no feature at 2; = 3. Prom 
Fig. 40 it is clear that the temperature peaks at z ^ 3 and that the gas 
becomes close to isothermal (adiabatic index 7 1), but the present 
constraints arc not sufficient to distinguish between a sharp rise (as 
indicated by the dashed line) and a more gradual increase. In addition, 
several different reionization histories and radiative transfer effects can 
fit the observed temperature data equally well, as is discussed in detail 
by Miniati ct al. (2004). 

The IGM temperature might be affected also by the energy de- 
posited by dark matter particles if these are either decaying or an- 
nihilating (the most interesting candidates being sterile neutrinos with 
masses of a few keV, light dark matter with mass ss 100 MeV and 
neutralinos) . If so, the IGM temperature might deviate from adia- 
batic evolution before the heat input from luminous sources comes into 
play (Mapelli, Ferrara & Pierpaoli 2006; Ripamonti, Mapelli &; Perrara 
2007a; 2007b; Valdes et al. 2007). 

5.5.8. Metal Enrichment 

Observationally, our best view of the IGM is still provided by the 
Lya forest in the spectra of distant QSOs and the detection of metals 
in the forest ranks as one of the most significant discoveries initially 
made possible by the Keck telescopes (Cowic ct al. 1995; Buries & 
Tytler 1996) and now also exploiting the power of VLT/UVES. From an 
analysis of these data, Hellsten et al. (1997) and Ranch, Haehnelt, &; 
Steinnietz (1997) concluded that the measured column density ratios 
N{C IV )/A^(H I ) imply a typical [C/H] ~ -2.5 at z ~ 3, with a one 
order of magnitude dispersion in the metallicity of different clouds 
about this mean value. These measurements, however, still refer to 
overdense regions of the universe, traced by Lya clouds with column 
densities in excess of log iV(H i ) = 14.5 . 

The situation is far less clear-cut when we turn to regions with 
log Ar(H I )< 14.5 and density increasingly closer to the mean (a good 

fitting formula is given by Pb/Pb — 0.8 [N-^ j / 10^'^ cm^'^]^ ''''), where 
the observations are very challenging even with a 10-m telescope. Two 
early studies addressed this problem, with conflicting results. Lu et 
al. (1998) applied a stacking technique to nearly 300 C iv regions in 
QSO spectra, but still found no composite signal. They interpreted 
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this non-detection as evidence for a highly non-uniform degree of metal 
enrichment at z ~ 3, with the voids having [C/H] < —3.5 . On the other 
hand, Cowic & Songaila (1998) used a pixel-to-pixel optical depth tech- 
nique to conclude that the average C iv /H i ratio remains essentially 
constant over the full range of neutral hydrogen column densities tested, 
down to log N{Il I ) ~ 13.5 . If this is indeed the case, the ejection and 
transport of metals away from galaxies must have been much more 
efficient than envisaged. Ellison et al. (1999, 2000) re-examined the 
problem and found that both approaches suffer from limitations which 
had not been properly taken into accoTint in previous analysis. 

To date, the number of direct metal line detections associated with 
log A^(H I )< 14.5 lines is small. For this reason, searches of warm 
photoionizcd O VI have been used. Schaye et al. (2000a) reported the 
first (statistical) detection of O VI in the low-density IGM at high 
redshift (2 < z < 4.5). By performing a pixel- by-pixel search for 
O VI absorption in eight quasar spectra, they detected O vi in the form 
of a positive correlation between the H I Lyct optical depth and the 
optical depth in the corresponding O vi pixel, down to tjj j « 0.1. 
However, the interpretation of the O vi data is not straightforward be- 
cause the ionization mechanism is often ambiguous. Simcoc et al. (2002) 
found evidence that the O Vl corresponding to strong H I absorbers 
(log A^(H I )> 15) closely resemble hot, colhsionally ionized gas found 
near regions of significant overdensity {6 > 300) with a tiny filling factor 
of about 10~^; furthermore, they seem to be clustered on velocity scales 
of 100 — 300 km and show weaker clustering out to 750 km s~^. 
The O VI corresponding to low H i absorbers (log N(H i ) < 15) on the 
other hand, is more consistent with warm, photo-ionized gas (Carswell 
et al. 2002; Bergeron et al. 2002; Simcoe et al. 2004). The presence of 
O VI has been detected down to rjj j « 0.2 also by Aracil et al. (2004), 
but associated only with gas located within !=a 400 km s~^ from strong 
H I lines, i.e. possibly associated with galactic winds. As a caveat on 
these studies one should keep in mind that such difficult experiments 
could be affected by a number of errors. For example, the apparent 
O VI absorptions could be due to spurious coincidences with H i 
absorptions at other redshifts and thus the O VI optical depth would 
be independent on that of the H i (Fieri &; Haehnelt 2004). 

Considerable progress has been made concerning the distribution 
of predominantly photoionized species, as C iv and Si iv . Schaye et 
al. (2003) have measured the distribution of carbon in the IGM as a 
function of redshift z and overdensity 6, obtaining the following fit to 
its median abundance: 

[C/H] = -3.471HI + Omt'olliz - 3) + Omt',]!i{log6 - 0.5), (20) 
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which is vaHd in the range 1.8 < 2; < 4.1 and —0.5 < log 5 < 1.8 (the 
quoted errors are statistical; systematic errors due to the uncertainty in 
the spectral shape of the UV background radiation can be substantially 
greater, see Schaye et al. 2003). Two points are worth noticing about 
this formula. The first is that even in underdense regions {6 = —0.5) 
carbon is detected with an abundance [C/H]?» —4. Secondly, it appears 
that very little evolution has occurred in such a long cosmic time span 
(2.16 Gyr). The slow evolution with column density and redshift is 
observed up to z = 5.7 (Ryan- Weber, Pettini & Madau 2006). 

The former evidence implies that metal transport must have oc- 
curred efficiently and at early times to simultaneously account for 
the cool temperatures and widespread distribution of elements like C 
and Si. Thus, observations of metals can possibly allow testing of the 
mechanisms that distributed them or establish when/where they were 
produced. For example, if metals are carried by galactic winds, they 
should be seen in absorption against bright background sources, such 
as quasars or GRBs, in narrow lines, with characteristic equivalent 
widths 0.5 A< < 5 A (Furlanetto k Loeb 2003). 

The "no evolution" trend seems to extend to even higher redshifts, 
as already pointed out by Songaila (2001) (and further supported by 
Pettini et al. 2003) , who concluded as well that the C IV column density 
distribution function is consistent with being invariant throughout the 
redshift range z = 1.5 — 5.5. This range has been extended to z !=a 6 by 
Simcoe (2006). This is quite a puzzling result as, even if the observed 
elements were produced by a high-z pre-enrichment episode, large vari- 
ations of the ionizing radiation field would be expected anyway during 
such extended time intervals. Alternatively, the C iv systems may be 
associated with outflows from massive star-forming galaxies at later 
times, while the truly intergalactic metals may reside in regions of the 
LycK forest of density lower than those probed up to now. Songaila 
(2006) finds that, in the redshift range z = 2 — 3.5, of a total sample of 
53 absorbers, only half of those with NC iv > 2 x 10^^ cm~^ could be 
associated with galactic outflows, while all the rest lie in the IGM. 

Only recently, since a larger number of observations has been avail- 
able, long-duration GRBs (which are thought to be associated with 
the death of massive metal-poor stars) have been used to probe the 
metallicity of the host galaxies. The available sample is still to small to 
draw conclusive results, but it looks like they can probe more metal-rich 
galaxies compared to QSOs absorption lines. Although these studies 
are still in their infancy, they look extremely promising in probing 
metal-enrichment (Savaglio 2006). 
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6. Conclusions 

The investigation of the first cosmic structures has started to attract a 
huge collective effort, as the distant universe is progressively disclosed 
by the most advanced observing facilities. Compared to just a few 
years ago, the progresses made in the understanding of several of the 
fundamental aspects in the field arc indeed impressive. It is then not 
inappropriate to claim that a qualitatively clear scenario of these early 
epochs, whose current status is broadly summarized in this review, 
is becoming rapidly established. However, as for any new frontier, we 
do not yet fully know its extent and the territories beyond it. The 
uncertainties in the basic cosmological constants are pinned down to 
the point that "precision cosmology" has become a fashionable way to 
present the data. Nevertheless, the shape of the primordial fluctuation 
spectrum is still uncertain and needs to be constrained, along with 
the degree of non-gaussianity of the density field. The nature of dark 
matter is likely to remain an extremely challenging problem, which 
may have intimate connections with the highly debated halo density 
profiles. Both problems affect the way in which galaxies build up and 
their observable properties; hence their importance extends beyond the 
realm of purely theoretical speculations. 

Aside from these general issues, the formation of the first structures 
exhibits its complexity in the large number of different aspects, ranging 
from galaxy formation to stellar evolution. In principle, the advantage 
of having to deal with the simplest native environment of the first stars 
(no heavy elements, dust, dynamically significant magnetic fields, UV 
radiation background) is rapidly dissolved within the lack of similar 
counterparts in the local universe that can be used as working labo- 
ratories. Although the formation of small mass stars in a primordial 
environment is possible, we tend to believe that the first stars were 
very massive, but as yet, no ultimate study has shown the formation of 
such stars starting from realistic conditions, including all the relevant 
physics (turbulence, radiative transfer, rotation) in detail. For this rea- 
son, we have essentially no clue regarding what the mass distribution, 
i.e. the IMF, of such stars could be; even indirect observations are 
of limited help due to the poor number statistics. Other constraints 
have been proposed, but they remain very weak. Also the evolution 
of massive stars is under work: for example, the combined effects of 
metallicity, rotation, and mass loss still need to be thoroughly studied 
and understood. 

These uncertainties must be added to the one related to the physical 
mechanism driving the characteristic stellar mass to lower (and more 
common) values with time. Metallicity increase is probably the most ac- 
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cepted explanation, but one can suspect that other agents may drive the 
transition: turbulence in the protostcllar clumps, strengthening of the 
UV radiation field, or heating due to cosmic rays. A sudden change in 
the properties of the first stars would noticeably affect the evolution of 
reionization history, as the ionizing power would drop dramatically with 
increasing (decreasing) stellar metallicity (mass). An equally uncertain 
quantity is the fraction of emitted ionizing photons able to escape from 
galaxies. Theoretical approaches to this problem have not yet clarified 
the dependence of such crucial parameter on redshift, galactic mass, 
and internal structure, and it is likely that the next major advance will 
come from a dedicated observational campaign. 

Despite the wealth of literature on the formation and evolution of 
low mass galaxies and on the variety of feedback effects which may 
hinder their formation, these remain controversial topics in the study 
of the first structures. Also, no general consensus has been reached 
on their contribution to the global star formation rate, IGM metal 
enrichment and reionization. For this reason, a stronger effort will be 
required in the future to clarify these issues. 

The widespread interest stimulated by reionization has produced, as 
a side effect, a rapid growth of numerical radiative transfer techniques. 
At the same time additional physical ingredients, like the plethora 
of feedback effects which are equally important as radiative transfer 
in shaping the reionization history, should be properly included in 
simulations. On the observational side, a major cornerstone to map 
reionization history will be provided by 21 cm HI line measurements, 
which hold the promise to open a paved road toward the distant "gray 
ages" . 

Similar arguments should also apply to the study of the intermediate- 
low redshift IGM, i.e. the Lya forest. Despite its simplicity, which 
allowed early successes in matching the observed properties in the 
framework of hierarchical cosmological models, few issues remain un- 
solved. A long standing puzzle involves the observed temperature dis- 
tribution, which seems to indicate an IGM warmer than predicted by 
numerical simulations. Also, an alleged temperature jump at z 3 
has not yet been either fully interpreted theoretically or confirmed 
experimentally. For this reason, it is crucial to devote more attention 
to the study of He ii reionization. In addition, the IGM temperature 
deserves further investigation as it is often used to set a maximum 
redshift for hydrogen reionization. 

Finally, the importance of winds from galaxies, both at intermediate 
and high redshifts, has only recently started to be appreciated, both 
theoretically and observationally. Do all galaxies go through a wind 
phase? What fraction of their mass did they lose in such events? What 
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fraction of the IGM has been inside galaxies, and possibly polluted by 
nucleosynthetic products of massive stars, as a function of rcdshift? 
What is the cosmic filling factor of metals and do we understand its 
evolution? The answers we currently possess on these and many other 
questions concerning the cosmic chemical evolution are very partial. 
The hope is that the pace we have kept so far will bring us close to the 
final answers during our life time. 
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